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Certainly new aviation projects take many 

forms, but manned, unmanned, research or 
commercial, they all depend on efficient 
hydraulics. 


The drawing board stage, or before, is the time 

to consult Dowty Equipment of Canada Ltd. . . 
a self contained organization dedicated to the 
highest standard of quality. 
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Mail coupon for Shell credit card 
—honoured by an ever-expanding network of airports 
throughout Canada and the U.S. 


Shell is North America’s 
largest supplier of commer- 
cial jeé fuel and aviation 
gasoline. Send for your Shell 
credit card—and charge 
famous Shell products at 
any of hundreds of airports 
across the continent. 

You can be sure that Shell 
quality will never vary no 
matter where you buy Shell 
products. 

HELL has pioneered many of the 

most important advances in 
aviation and aerospace fuels and 
lubricants. 100-octane gasoline 
is one example. Another is the new 
ashless dispersant oils now widely 
used by piston engine operators. 
Yet another example: the rocket fuel 
used to put U.S. satellites in orbit. 

The same kind of research that 

produced these products has gone 
into all Shell fuels and lubricants in 
use today. Here is a brief rundown 
of these outstanding aviation 
products: 
Shell aviation fuels. Shell gasolines 
areavailablein grades 80/87, 100/130, 
and 115/145. Also available: Aero- 
shell* turbine fuels 650 and JP-4. 
Aeroshell Oils. Aeroshell Oil W is 
Shell’s unique non-ash, dispersant 
oil with multi-viscosity characteris- 
tics. Use grade 80 for small engines, 
grade 100 for large engines. 

Shell also makes outstanding 
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straight mineral oils in grades 65, 
80, 100, and 120. 

Aeroshell Fluids. Shell offers you 
nine Aeroshell lubricating fluids— 
everything from squirt-can oil to a 
low-volatility instrument oil that is 
effective.to —65 degrees F. 
Aeroshell Greases. The nine Aero- 
shell greases include a high tempera- 
ture grease, a graphite grease, a 
grease with an extreme pressure 
additive, and a special grease for 
helicopters. 

Aeroshell Compounds. These include 
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two anti-corrosion compounds, a 
de-icer, an anti-frost compound, and 
an anti-seize compound containing 
graphite. 

You can charge these and many 
other Shell aviation products with a 
Shell credit card. 

You can also charge a number of 
services, according to the conditions 
specified on Shell’s credit card. 
These services include minor emer- 
gency repairs; transient hangar and 
tie-down fees ; landing fees; and tires, 
batteries and accessories. 


--——— Mail coupon today for Shell credit card ----- 
Shell Oil Company of Canada, Limited 


Gentlemen: | would like to have a Shell credit card. 
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Versatile New Engine 
by CANADIAN PRATT & WHITNEY 


Canada’s first turboprop is a radical departure 
in engine design. The PT6 is an axial-centri- 
fugal flow, moderate pressure ratio turbine 
engine with a free turbine drive for either 
turboprop or turboshaft installations. This 500 
h.p. power plant is ideally suited for single or 


multi-engine fixed wing aircraft, helicopters or 
high speed marine installations. From a func- 
tional point of view, it weighs just 250 pounds, 
yet combines rugged service and dependable 
performance with the practical economies of 
low fuel consumption and ease of maintenance. 


Conadion [eats Whithey Udircraft- 


COMPANY, LIMITED Longueuil, Montreal, P.Q. 
A subsidiary of United Aircraft Corporation, East Hartford, Conn. 
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Professor T. S. Lowe’s first balloon ascent 
over Ottawa on the 17th July, 1858. This 
ascent is believed to have been the first 


man-carrying flight in Canada. 














EDITORIAL 


ACHIEVEMENTS IN COOPERATIVE RESEARCH 
AND DEVELOPMENTS IN NATO 


I AM particularly pleased to have been asked to con- 
tribute to your Journal in a guest editorial capacity. 
The reason for this is not only that I have another 
opportunity to develop a topic that I have been ad- 
vocating for many years but because, judging from 
the response Canada has made through her participa- 
tion in cooperative international endeavour, I have 
readers who are already largely convinced. I have 
been the Chairman of the Advisory Group for Aero- 
nautical Research and Development and, after some 
ten years of experience in seeing: put into practice, in 
a significant way, those terms and expressions every- 
one recognizes as being desirable aims, I am more 
than ever convinced that international cooperative 
activity at the most creative level is correct in prin- 
ciple and, much more significant, possible in practice. 
Ten years ago we were told that the nations would 
not be able to cooperate effectively, because reasons 
of national security, national proprietary interests, na- 
tional pride, creative individuality etc would prevent 
any real achievement. In spite of these counsels we 
persisted in the formation of our Advisory Group, 
with its permanent secretariat in Paris and its volun- 
tarily subscribed panels. From the beginning, at which 
time we had 3 panels drawing together some 80 
talented NATO nationals in the study of three broad 
scientific disciplines, we have built up, by 1961, 7 
panels involving over 250 of among the foremost 
scientists. of the NATO nations. The thought and 
effort of these people is put voluntarily at the dis- 
posal of NATO in meetings normally held twice a 
year; the work of each panel is monitored by an 
executive officer, who is a member of the permanent 
secretariat located in Paris. Dr. John J. Green, with 
his early enthusiastic support of AGARD, was a great 
help in this development. In 1955 we were able to 
hold the Fifth General Assembly in Ottawa, and were 
fully supported by your Government. Examples of 
achievements to which Canada, as one of the NATO 
nations, can point with pride are various and occur 
at all scientific levels and in every scientific discipline: 
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— Work on the NATO Light Strike Fighter 
— The Training Center for Experimental Aerodyna- 
mics, Brussels 
— Design criteria for aircraft undercarriages; rough- 
ness of NATO runways 
— Missile Range in Greece 
—A large number of publications which provide a 
review of certain research endeavours and often 
give member countries access to information other- 
wise only available through costly and time con- 
suming research programs. Examples: 
Flight Test Manual 
Manual on Aeroelasticity 
Materials Handbook 
Multi-Lingual Dictionary 
— Symposia on Operational Research 
— Work in the area of documentation. That the ac- 
celeration of scientific activity has resulted in an 
overwhelming increase in the quantity of technical 
literature is apparent to everyone — what is not 
exactly clear is the response that should be made 
to the problem, in the international sense. 

The list of accomplishments is long — and grati- 
fying. It is gratifying to me personally but this is of 
small importance. It is more important that it be a 
source of gratification to the nations, such as yours, 
which made it possible. Make no mistake, the achieve- 
ments did not come about through some form of in- 
ternational prestidigitation: no magic formulae have 
produced these results. In each case dedicated men of 
many nations have given of their time and effort to 
find a solution to a scientific situation, exactly as has 
always been done in national practice; having found 
a possible approach, they have applied the normal 
development routines to achieve the progress we now 
see. The only magical element that can be found is 
that, in many cases, the result has the status of an 
agreement among nations, or a course of action which 
is acceptable to these nations. When one realizes that 
it took years to agree on a standard atmosphere in 
the West, one begins to appreciate the significance of 
what has been accomplished. 
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To the individuals who have worked on your be- 
half within our framework we all owe many thanks 
— not simply because they have displayed a great deal 
of energy within our organization but because they 
have helped to show the way. In most panels, there 
are Canadians directly involved. I cannot name all of 
the excellent participants from Canada. I shall name 
only a few. Dr. MacPhail of the National Research 
Council, who is Chairman of the Fluid Dynamics 
Panel; Mr. Wood of the National Aeronautical Estab- 
lishment, Chairman of our Flight Mechanics Panel. 
To Mr. F. R. Thurston of the National Aeronautical 
Establishment, we all owe much for the application 
of his considerable energy and talent to the Structures 
and Materials Panel, for five years as a member and 
then as Chairman. The same panel also has an excellent 
executive officer from Canada, Mr. Gerald H. Cooper. 
The work of these men indicates that in Canada, as 
in the other nations which are contained in our struc- 
ture, significant progress in cooperation on an inter- 
national scale is inevitable, providing that (a) there 
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W: are indebted to Mr. K. M. Molson, Curator of 
the National Aviation Museum, for our frontis- 
piece this month. 


The complete story of this balloon ascent has 
not yet been assembled but it seems that Professor 
T. S. Lowe was an American, invited to Ottawa to 
take part in the celebrations marking the laying of 
the first trans-Atlantic cable. (A balloon ascent was 
justifiably regarded as a telling feature of any cele- 
brations in those days, whether it was appropriate 
or not.) Though his balloon had been built in the 
States, apparently Professor Lowe had not made an 
ascent before coming to Ottawa, so that the flight in 
Ottawa on the 17th July, 1858, was his first. It is re- 
ported that he made other ascents in the area that 
summer. 

It is believed to have been the first man-carrying 
flight in Canada. 

Later Professor Lowe became possibly the fore- 
most US balloonist and was placed in charge of the 
balloon unit of the Union forces in the American Civil 
War. 


Those who know Ottawa can get their bearings 


exists a strong structure of organization within which 
their combined talents can find application, and (b) 
the value of their work within this structure is so fully 
appreciated nationally that they are supported in their 
participation. 

In closing, I should like to sound a note of cautious 
optimism. In my almost sixty years of intimate in- 
volvement in science, I have never lost sight of the 
possibility for good of international actions — we are 
further toward this goal now than I would have 
thought possible even fifteen years ago. But we have 
a tremendous distance yet to go. All that can be said 
is that a start has been made; a possibility has been 
demonstrated; a series of tests, statistically not suffi- 
cient, has confirmed an hypothesis. The future is up 
to the workers of the nations, which, like yours, show 
a will. 


Dr. THEODORE VON KARMAN, 
Chairman, Advisory Group 
for Aeronautical Research 
and Development, NATO 
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Prof. Lowe’s first balloon ascent over Ottawa on 
the 17th July, 1858. 


from the Rideau Canal Locks. The picture was ap- 
parently “taken” from the site of the present Parlia- 
ment Buildings, looking east. Though appearances 
may be misleading — and probably are in this case — 
Professor Lowe seems to have taken off, appropriately 
enough, from the neighbourhood of Rockcliffe. 
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MAN-POWERED FLIGHT 
AND TODAYt 


by Helmut Haessler 


Two of the greatest pioneering achievements in 
the realm of man-powered flight were the Haessler- 
Villinger and Bossi-Bonomi aircraft of the mid-30s. In 
the December 1960 issue of the Journal we published 
a paper by Mr. Bossi reporting on his work and we 
have now received this paper from Mr. Haessler 
giving details of his aircraft and some stimulating ob- 
servations about future possibilities. — Sec. 


SUMMARY 


This paper reviews the flights of the Haessler-Villinger 
man-powered aircraft in the years 1935-37 and gives the 
specification of this aircraft. It considers the influence of 
weight and span on pene, and wing weights are calcu- 
lated for a range of spans, areas and loads, on the basis of 
the weight of the Haessler-Villinger aircraft. It concludes with 
a proposal for a man-powered aircraft capable of restricted 
flights of 1 hour, with and without ground effect, by a trained 
racing cyclist, or unrestricted flights with glide ratio improved 
to 1:100, by an average amateur. 


INTRODUCTION 
eo the early aeronautical pioneers tried to fly, 
they thought the use of man-power would be 
the quickest and easiest way to do it. For centuries 
they were unable to succeed by this means and the 
first actual flights were made with gliders, without 
any power at all. When power was used, they found 
that much more was needed than a human being could 
develop. Aircraft powerplants, which in the early days 
developed about 10 hp, now have more than a 
thousand times this power output, making it possible 
to fly beyond the limits of the atmosphere into outer 
space. 


It seems that with the beginning of the space age, 
which coincides with the completion of the conquest 
of the air, interest in man-powered flight is becoming 
livelier than ever. There can be no practical reason 
for this, unless man considers the conquest of the air 
to be incomplete until he can fly with his own muscles. 
It is not surprising that this could not be done at the 
beginning of the development of aircraft, since a 
horizontal flight with that wee bit of man-power de- 
mands extraordinary refinements in the airplane shape 
and structure. 


I have been asked to report on the development 
and flying of the Haessler-Villinger man-powered air- 
craft 26 years ago and I would like to add some 
thoughts about the present situation, 





+Received on the 17th January, 1961. 
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IN 1935-37 





Figure 1 


FIRST HORIZONTAL MAN-POWERED FLIGHTS 
First flight 

On the 29th August, 1935, at the airport of Frank- 
furt/M in Germany, the Haessler-Villinger man- 
powered aircraft was carried out of the hangar to the 
airfield to attempt its first man-powered flight; some 
previous test flights had been made without power on 
the airport at Dessau. The airplane looked like an 
ordinary sailplane of 45 ft span, except for the pro- 
peller mounted on a pylon over the fuselage (Figures 1 
and 2). The technical specification is given in Table 1. 





Figure 
Man-powered aircraft — Haessler-Villinger 
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TABLE 1 
' SPECIFICATION OF THE HAESSLER-VILLINGER AIRCRAFT 


Wing area | 9.65 sq m | 104 sq ft 
Span | 13.5 m | 44.3 ft 
Wing loading | 11.5kg/sqm _ | 2.29 Ib/sq ft 
Aspect ratio 18.8 18.8 


Wing section Gottingen 535 


Rib distance | 20 cm 7.9 in 
Chord | 76cm | 30.0 in 
Sinking speed | 0.52 m/sec 20.5 in/sec or 102.5 

| ft/min 
Fuselage cross-section | 0.33 sq m | 3.56 sq ft 
Prop dia. | 1.5m | 4.9 ft 
Prop rpm 600 rpm | 600 rpm 
Prop eff. | 0.82 0.82 
Glide ratio | 1:24 | 1:24 
Empty weight | 34 kg | 75 Ib 
Energy acc. 10 kg 22 Ib 

| 143 Ib 


Pilot 65 kg 








A prize had been instituted by the Polytechnische 
Gesellschaft of Frankfurt for the first flight by man- 
power. There were regulations covering the distance 
to be flown and the storage of the energy for takeoff. 
In spite of the fact that I could not agree to all points 
of the regulations (the weight of the energy accumu- 
Jator of 20 Ib had to be carried in the aircraft and 
the flight path had to go around two marks 1700 ft 
apart) the aircraft was designed to meet the regula- 
tions as far as possible. 

In preparation for the takeoff the airplane was 
anchored to the ground by a pair of anchor-pins at 
the centre of the fuselage, and the pilot stretched 
several rubber cables, one after another, from the 
fuselage to a second anchor at a certain distance in 
front of the aircraft. Then he entered the cockpit and 
released the anchor at the fuselage. The aircraft ac- 
celerated and climbed to a height of about 10 ft. At 
this point the pilot put his power into two pedals 
connected to the propeller and continued to fly hori- 
zontally for a distance of about 400 ft. When his peak 
power output was over he relaxed and finished his 
flight in a normal glide. 

The total distance covered from takeoff point to 
standstill was 780 ft and was checked by a committee 
consisting of members of the Physikalische Gesell- 
schaft and other officials. After consideration of all 
circumstances, this committee pronounced this flight 
as the first flight over a horizontal airfield accom- 
plished by the unaided power of the pilot. 

Later calculations and tests showed that on this 
flight at least 420 ft were covered in horizontal flight 
under the power of the pilot applied to the propeller, 
and not more than 360 ft resulted from the stored 
energy. Since the storing ability of the accumulator 
changed due to hysteresis in the rubber, and friction 
factors between the skid and ground changed, the 
distance of 360 ft is an optimal figure. Test flights 
without using the propeller covered distances between 
280 and 360 ft, comprising the acceleration path on 
the ground, the climb path until the stored energy 
was used up, and the glide path without power 
(Figure 3). 

When the news of this first flight went out, people 
said that such a short flight with the aircraft launched 
by a rubber bungee, could have been done with any 
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good glider without man-power. When it was pointed 
out that a glide of such a distance could not be made 
from a launching height of 10 ft, unless the aircraft 
had a glide ratio of at least 1:62, it was argued that 
the launch had been done at excess speed, so that 
the glide distance could be attributed to kinetic 
energy. A calculation of this condition shows that 
launching at excess speed with a given amount of 
stored energy is less efficient and results in shorter 
flight distances than those actually flown. Besides that, 
the aircraft in question was not capable of storing the 
necessary energy in excess speed, according to its 
lift/drag diagram and its flight envelope, since it would 
have fallen apart before that point could be reached’. 

This illustrates the situation at that time, when a 
clear definition of man-powered flight had not been 
established. The committee should have said that for 
the first time a horizontal flight had been achieved 
under the direct power of the pilot, for a distance of 
at least 420 ft; the total flight of 780 ft included a 
takeoff and climb for a distance of 120 ft, using the 
pilot’s energy (stored in an accumulator), and a glide 
of not more than 240 ft. We will see later that 
elongation of the glide path by continuing to apply 
cruising power was not feasible. 


Subsequent flights 


After several improvements to the aircraft (in par- 
ticular, controls) and changing the pilot for a stronger 
one, further flight tests were made in August 1936 
at the airport at Frankfurt and a total distance of 
1300 ft was covered. The horizontal flight path was 
at least 940 ft compared with 420 ft of the year before. 

A few months later on the 21st November, 1936, 
the aircraft was flown at Hamburg before members 
of the Lilienthal Gesellschaft. A short movie strip 
from a newsreel which shows one of these flights is 
still in existence and shows a horizontal flight path of 
more than 1060 ft, of a total distance of 1420 ft. From 
Figure 4, taken from this strip, it can be seen that, 
after recovering the normal height of about 10 ft, 
the pilot maintains this height for the length of the 
man-powered flight. 


In the following year we demonstrated the aircraft 
at three air shows at Schleusingen, Schwarza and 
Leipzig, each time covering distances from 1400 to 
2000 ft. We would not have agreed to undertake these 
scheduled demonstrations if we had not been confident 
of achieving this performance regularly even in poor 
weather conditions. At Schleusingen, the flights had 
to be done in the rain. 


}-————— WOR! ZOWTAL FLIGHT WITH MANPOWER — —— 
| FIRST FLIGHT 420 FT 
LONGEST FLIGHT 2000 FT 











DISTANCE FLOWN WITH STORED POWER — 
280 - 360 FT 


TOTAL DISTANCE ————- -—_——___- ____ 
FIRST FLIGHT 780 FT 
LONGEST FLIGHT 2360 FT 


Figure 3 
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Figure 4 


Taken from a movie strip made at Hamburg in November 1936 


On the 4th July, 1937, at the airport at Meiningen, 
a distance of 2350 ft was attained and measured by 
airport officials. This included a horizontal flight of 
2000 ft flown in 1 min 12 sec at a height of 10 ft. The 
horizontal part of this flight was five times longer than 
that of the first flight in 1935. These results were re- 
ported in the newspapers, but they are apparently not 
known in aviation circles today. 


A weight of 20 lb, equal to the weight of the 
rubber bungee and the winding equipment, was car- 
ried at all times in the nose of the fuselage and could 
not be removed without seriously affecting the posi- 
tion of the centre of gravity — especially not after the 
alteration of the tail surfaces in the second year, 
which resulted in higher tail weight. In all cases the 
stored takeoff energy was the same; the accumulator 
system was not changed. 


These results were the best we could expect with 
this aircraft and with this pilot. We had no intention 
of looking for another pilot, although this could have 
led to flights of 2 minutes with this aircraft. In the 
end we could not afford the experiments any longer 
and gave the aircraft to the Luftfahrt Museum in 
Berlin, in February 1938. 


Between July 1935 and February 1938 a total of 
about 120 flights were made with this aircraft. About 
80 flights were made with man-power; the others were 
done for testing changes in trim and controls and 
training two pilots to the peculiarities of man-powered 
flight. Of the 80 man-powered flights at least 15 were 
of .a duration of 1 minute or more. 


Other flights 


There were other man-powered flights before and 
after ours. In 1927 an airplane with flapping wings 
was designed by Mr. Lippisch and flown at the 
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Wasserkuppe Soaring Centre. Several of its flights 
were stretched glides up to double the normal dis- 
tance, on a not horizontal terrain. Takeoff was ac- 
complished with rubber bungee and teams. 


In Italy a man-powered aircraft was built by 
Bossi-Bonomi in 1936. Takeoffs were made with 
teams’ and the results were stretched glides. The 
longest flight, with assistance by team to a height of 
25 ft, was % of a mile. In some cases, takeoffs were 
made under the power of the pilot alone, through a 
driven wheel and the propellers. The flying distance 
after this type of takeoff is reported by Mr. Bossi to 
have been 300 ft. 


Today in 1961, 26 years after our first man- 
powered flight, another prize is being offered, by Mr. 
Kremer, and the competition will be conducted by 
the Royal Aero Club of the U.K.* The regulations do 
not allow the use of stored energy and specify a flight 
in a figure of eight around two marks half a mile 
apart, which will require a flying time of between 3 
and 4 minutes. It seems that the regulation banning 
stored energy was made because of the recording 
difficulties which could arise from its use. 


PERFORMANCE LEVELS IN MAN-POWERED FLIGHT 


Now it would be appropriate to establish defini- 
tions of the different performance levels in man- 
powered flight. 


When we consider a glide elongated by applying 
man-power, the lower limit is the normal glide with 
zero power, the upper limit horizontal flight. The 
range between these limits can only be considered as 
an intermediate state, which I will call “the man- 
powered glide”. The aircraft of Lippisch belongs to 
this category and most of the flights with the Bossi- 
Bonomi aircraft were done in this regime. 
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Figure 5 
Scheme for calculation of takeoff height 


Only the horizontal part of a flight with man- 
power can be considered as “man-powered flight.” 
If stored or external energy is used for takeoff, the 
amount of this energy has to be checked and the part 
of the flight attributable to it has to be deducted from 
the whole. Another way would be to record the 
height of the flight-path by photography and to time 
only that part of the flight which is above release level. 
This leads to practical difficulties when the release 
level is not at the eye level of the recording person. 
The ideal case is the takeoff under the direct power 
of the pilot, since here the whole flight is above take- 
off level, if the airfield is level. 


I want here to inject a proposal for defining assisted 
takeoff for any particular aircraft by establishing, for 
that aircraft, a specific limited amount of stored 
energy used for takeoff and the best use, in terms 
of flying time, that can be made of that energy. 


The aircraft should be launched on a level airfield 
so that the height at power shut-off is near zero — 
let us set an arbitrary maximum of 8 ft on this release 
height, sufficient to give the pilot a margin for the 
correction of slight deviations from the horizontal 
flight path. 

Release height H (Figure 5), defined as the dis- 
tance between the lowest point of the aircraft and 
the ground, is to be calculated from the total distance 
A flown with the stored energy, the glide ratio (tan €) 
of the aircraft under conditions of ground effect, and 
the horizontal distance B from start to power shut-off. 
The best flying time from start to finish (the dis- 
tance A) must be established by a series of test flights, 
all done with the same amount of stored energy. 


The time so established must be deducted from 
the times of subsequent flights using man-power and 
launched with the same amount of stored energy. The 
result in each case would be the time of the horizontal 
part of the flight, achieved under pilot-power applied 
to the propeller. 

Once the figure for the duration of that part of 
the flight attributable to stored energy has been estab- 
lished, no distances would have to be scaled and a 
stop watch would be sufficient for recording the man- 
powered element. Under this arrangement the energy 
storing equipment would never have to be carried in 
the aircraft; in fact all considerations affecting the 
takeoff would be standardized _and established as a 
datum, from which the “man-powered flight” can be 
measured. 

Figure 6 shows the performance characteristics of 
a professional racing cyclist. We know that a man 
has a certain power level which he can maintain for 
a long time, up to several hours. By working at this 
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cruising power level, at which the lungs are capable 
of maintaining the oxygen supply required by the 
body, no energy reserves are used up. Man has a 
certain capability of storing oxygen in his blood, 
which enables him to do a certain amount of work 
without adding new oxygen by breathing. This 
reserve of energy is of the order of 1100 m kg. 
This energy can be used in a very short time, e.g. 
10 seconds, or it can be spread over hours. If the total 
power output is distributed over more than 5 minutes, 
the expenditure of reserve energy becomes very little 
compared with the cruising power. We will consider 
the power up to 5 minutes as the short-time power 
and over 5 minutes as the long-time power. 


Performance Level 1 


The horizontal flights 25 years ago were short-time 
flights and done at heights of less than the wing span, 
that is within ground effect. This will be considered 
as Performance Level 1 (Figure 7). The Haessler- 
Villinger aircraft is typical of this group with assisted 
takeoff and Bossi-Bonomi with unassisted takeoff. 


Performance Level 2 


If a man-powered aircraft is to be built today, its 
level of performance should be an improvement on 
that of 25 years ago. Horizontal flights of 1 minute 
duration have already been made. If a new aircraft 
could fly for say 2 minutes, it would be an improve- 
ment but it would not be a significant advance in the 
art of man-powered flight, and would scarcely be a 
sufficient reward for the effort entailed. 


To achieve real progress the power necessary to 
maintain horizontal flight must be equal to or less than 
the power output of.a man for 5 minutes or more. 
We will see later that it seems to be possible today to 
fly for several hours at low heights, with the help of 
ground effect, using stored energy for takeoff in the 
manner described above. That will be Performance 
Level 2. To be realistic, the regulations of any com- 
petition today should conform to this level. 
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Figure 6 
Typical man-power envelope 
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Performance Level 3 

If the drag of the aircraft can be reduced to the 
point where the power necessary to maintain hori- 
zontal flight at a height equal to or greater than the 
span of the aircraft, ic. out of ground effect, is no 
greater than the cruising power output of a man, we 
would come to Performance Level 3. This perform- 
ance looks feasible today for a limited number of 
pilots. If it can be accomplished, the takeoff problem 
will be automatically solved, since the power necessary 
for horizontal flight at ground level will be lower 
than that for cruising out of ground effect; the pilot 
will be able to take off without using up his reserves 
and climb to a height equal to the span, where he can 
continue to fly for several hours. 


Some hold the view that only flight with unassisted 
takeoff should be considered as man-powered flight. 
However, if that part of a flight which is directly 
attributable to stored energy or other assistance is 
checked in a reliable way and deducted from the 
whole distance flown, there can be no reason why 
the remainder should not: be considered as a man- 
powered flight. It seems to be unreasonable to limit 
a man-powered aircraft to short-time flights by a pre- 
mature requirement that it should take off under direct 
pilot-power, when man-powered flights lasting for a 
respectable time could be made with assisted takeoff 
as described. During a takeoff run, the pilot would use 
a great part of his reserve energy, limiting his time of 
actual flight, whereas, with assisted takeoff, he would 
fly on cruising power alone. Takeoff under direct 
pilot-power will come into its own as soon as it is 
possible to fly at cruising power without ground 
effect. 


Performance Level 4 

If it were feasible to fulfill strength and weight 
requirements for a combination of soaring and man- 
powered flight, then we would get Performance 
Level 4 for unrestricted flights. To perform at this 
level a glide ratio of 1:75 with man-power should be 
considered as a minimum. 


MAN-POWER 

I made my first tests on available man-power, with 
cyclists in 1933. The cyclist to be tested had to pull 
a second cyclist of the same weight and in the same 
attitude and position, giving the same conditions for 
both (Figure 8). The second cyclist did no work and 
had only to watch a scale indicating the tension in the 
connecting cable. A rubber damper was included in 
the connection, to absorb shocks caused by roughness 
of the road. A third person with a stop watch checked 
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the time taken to cover a measured distance. Since the 
power used to pull the second cyclist was equal to the 
power to move the first, the total power output was 
twice the amount shown by the scale and stop watch. 
The whole setup was simple and worked very satis- 
factorily. 


Finding a pilot 

I tested first myself, then some friends, then every 
available member of my gliding club, and the results 
were discouraging for use in man-powered flight. 
Then I went to a camp of professional racing cyclists; 
and that made ends meet. I found a fellow who was 
willing to work with me, who developed a power 
output of 1.3 hp for 30 seconds. I made arrangements 
with my gliding club that this man should be taught 
to fly as quickly as possible. When I started designing 
the aircraft for him, I laid it out in such a way that 
I would get as wide a power margin as possible, to 
allow for still unknown losses. I realize now how near 
I was, theoretically, to achieving flights up to 5 
minutes with this design and this cyclist. 

The first difficulty came when I filed an entry for 
the competition and the committee insisted that the 
accumulator for takeoff power should be carried in 
the aircraft, this meant an additional weight of 20 lb. 
I had to decide whether to drop the project altogether 
or to add the weight, which would make the chances 
of complete success very marginal. Success meant not 
the fulfillment of the requirements of the competition 
— this could not be done — but rather the achievement 
of the first horizontal flight by man-power; time and 
distance were of secondary importance. Here my 
friend, Franz Villinger, came to my assistance and we 
decided to continue together. All the work was done 
at night and on week-ends. 


When the aircraft began to take shape, I found 
out that our pilot was not making sufficient progress 
with his flying lessons; as a professional cyclist, he was 
more interested in cycling than in flying and, because 
of his profession, he was out of town most of the time. 
Finally I had to look for another pilot with gliding 
experience. Since I could not find the right man in 
my own club, I went to other clubs and eventually 
settled on a Mr. Duennebeil. His power output was 
far too low for our purpose (60% of the professional 
cyclist) but he was the best man available. It was he 
who made the first man-powered flight at Frankfurt 
in 1935. The fact that we were able to do this even 
with this low level of power output proves that the 
design was based on reliable. data, though in his 
Mitchell Memorial Lecture at Southampton* Mr. B. S. 
Shenstone is reported as saying “The Haessler-Vil- 
linger design . . . was based on quite false data on 
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Figure 8 
Power test with cyclists 
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Figure 9 


the power a man could produce”, without any sup- 
porting figures. After the first flight we started 
testing again, in search of a better man. Since the 
glider pilots did not look promising and the profes- 
sional cyclists were too hard to get, we looked for 
amateur racing cyclists in different towns and selected 
a Mr. Hofmann as our next pilot. His power output 
was about halfway between those of Duennebeil and 
our first professional. Hofmann was unemployed at 
this time and we sent him to a flying school at our 
expense; then we looked for a job for him in Dessau, 
where Villinger and I were living. 

Altogether we tested about 50 cyclists. In the first 
year, we ran our tests from 20 seconds to 40 seconds, 
but in the next year we extended them up to 5 minutes. 
Our discovery of the sharp power drop in the first 
few minutes was discussed with Mr. Ursinus and he 
arranged to carry out some further tests in the 
“Muskelflug Institute” which he had founded in 1935. 

At the time of our first man-powered flights, 
doubts were expressed about our pilot seating ar- 
rangement and many considered the rowing move- 
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Figure 10 
Available man-power 
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ment to be more efficient than cycling. Mr. Ursinus 
investigated the different types of possible movements 
and seating arrangements and concluded that the row- 
ing movement was less efficient than cycling and that 
our seating arrangement was superior both to the 
normal bicycle arrangement and to a prone position 
of the pilot. In addition, the position of the pilot in 
our aircraft was unique in not requiring the pilot to 
use his hands to balance himself, so that no reaction 
loads due to balancing were conducted into the con- 
trol system (Figure 9). 

Mr. Ursinus also found out that a combination of 
cycling with the legs and hand-cranking gave a signi- 
ficantly higher power output for periods of less than 
2 minutes. 
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Figure 11 
Available man-power 


Summing up my own results, I found that, ac- 
cording to his state of training, every person has a 
certain cruising power level where oxygen exchange is 
balanced and no reserves are used up. For a pro- 
fessional cyclist, the energy reserves of the body, 
above cruising power, are of the order of 1100 m kg 
or 8000 ft lb and his cruising power is about 0.46 hp; 
a trained amateur cyclist has a cruising power of 
0.37 hp and an average cyclist about 0.28 hp. This 
refers to an average pilot weighing 65 kg. I was the 
first to publish these figures, in 1937°. 


Cruising power and reserve energy 

The energy reserves above cruising power level 
can be used at different short-time power levels 
(Figure 6). The envelope of these short-time power 
levels is a hyperbola. When a cyclist has maintained a 
certain level of power output until it reaches this en- 
velope, he is exhausted and has to relax; he cannot 
continue following the envelope line; his power out- 
put drops sharply and then, as his oxygen level is 
restored, it returns gradually to cruising power level. 
In flights of less than 5 minutes this recovery would 
not be possible before the aircraft had lost height and 


Canadian Aeronautical Journal 
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20.38 0.33 0.30 0.29 0.28 

0.55 0.46 0.41 0.39 0.38 

0.79 0.60 0.52 0.48 0.47 
RECORD NATIONAL AMATEUR 


RECORD WORLD AMATEUR 1.03 
RECORD WORLD PROFESSIONAL 1.16 0.89 
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Figure 12 


landed. The recovery line is difficult to establish, but 
experience shows the trend. 

At the beginning, I did not know the nature of this 
power decay and, when I found that the power drop 
in the first two or three minutes was greater than 
expected, the tests with the two cyclists were extended 
up to 5 minutes. Beyond 5 minutes the recording of 
the tension in the tow cable became less and less 
reliable and I had to extrapolate the results to one 
hour after repeated tests to establish the reliability of 
the figures below 5 minutes. 


Mr. Nonweiler has given the power used by 
cyclists in establishing national and world records’. 
In Figure 10 I have plotted cruising power against 
energy reserves from zero to world record perform- 
ance. The results show a close relationship between 
these two factors so that it is now possible to calculate 
the power decay against time for each cruising power 
level (Figures 11 and 12). To find the performance 
of a cyclist, only one test of over 1 minute is neces- 
sary, and the power decay against time can then be 
read from the graph. 

Figure 13 shows available cruising power against 
weight of cyclist. This is based on the experience 
that persons weighing about 140 Ib have the best 
power/weight ratio. Above and below this weight this 
ratio becomes less favourable. 

All the tests referred to were done with cyclists. 
The efficiency of the chain drive of the bicycle was 
not checked and separated from the results, so that 
the actual man-power was in fact higher than in- 
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Figure 13 
Power/weight ratio for cyclists 
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dicated, by the amount of the frictional losses. This 
must be borne in mind when calculating the efficiency 
of an aircraft drive system. 


EFFICIENCY OF THE PROPULSION SYSTEM 

I want here to comment upon a remark attributed 
to Mr. Shenstone, in his Henson and Springfellow 
Lecture entitled Man-powered Flight’, that the power 
transmission in the Haessler-Villinger aircraft was 
inefficient because it incorporated a belt drive. He 
claimed that a chain drive would be more efficient, 
were it not for the twist in the chain necessitated by 
the two axes of the power drive having to be arranged 
at right angles to each other (Figure 14). 

The efficiency of belt drives depends in the first 
place on the desired lifetime of the belt. When a belt 
is to be used for 40 hours a week in commercial work, 
it should have a lifetime of several months without 








PROPOSAL NONWEILER 


Figure 14 
Propulsion systems 


maintenance. The materials used to meet this require- 
ment are limited in flexibility and give efficiency co- 
efficients below those of a straight chain drive. In a 
drive for a man-powered aircraft conditions are quite 
different; its life is counted in minutes and it is under 
constant maintenance. We tried to find a material with 
the highest flexibility and good surface friction, and 
were not interested in long life. First we tested all 
available commercial belt materials, none of which 
were satisfactory. Then we started to develop our 
own, using fabric belting of different gauges, differ- 
ent widths and different types of weaving, which we 
rubberized with a special solution to get less slip and 
better contact with the pulley. The pulleys were 
covered with a layer of rubber on a metal base. In 
this way we evolved a belt drive of an efficiency 
better than that of a straight chain drive. The only 
trouble we experienced was stretching after several 
flights, due to the fabric base being made of a natural 
fibre. We changed the belt after every six flights. 
Today materials like nylon would make this un- 
necessary. 
The efficiency of the propulsion system of our 
aircraft has been unjustly criticized and in defence I 
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Figure 15 
Efficiency of propulsion systems 


would compare our arrangement with the unconven- 
tional systems of the Bossi-Bonomi aircraft and the 
Nonweiler proposal. 


To determine the total efficiency of a propulsion 
system we must consider not only propeller efficiency 
and frictional losses in the drive system (for which 
only the part above the level of the losses in a bicycle 
are significant) but also the power necessary to keep 
the propulsion system aloft; this amounts to about 
0.005 hp/lb. We should not overlook the fact that 
the propulsion system is not a part of the structure 
of the airframe and, to support and accommodate it, 
the structural weight of the aircraft is increased by 
some 80% of the weight of the propulsion system 
itself. We will call this part the indirect weight of the 
propulsion system. This weight factor is so important 
that, if we are to fly at cruising power, the direct 
weight of the propulsion system, comprising the whole 
drive system, bearings, fittings, fillers, doublers and 
propeller should be limited to 6 lb, since the loss due 
to the combined weight, direct and indirect, should 
not exceed 10% of the power available. 

At 1 hp the total efficiency of the Haessler-Vil- 
linger drive system was 78 %, compared with 49% for 
the Bossi-Bonomi system (Figures 15 and 16). At 
cruising power (0.5 hp), we get for the Haessler- 
Villinger 75% and for the Bossi-Bonomi 17%. In the 
latter the power loss due to weight alone amounts to 
64%, compared with 7% for the Haessler-Villinger. 
This shows clearly that the influence of the weight of 
the propulsion system can easily offset any efforts to 
increase propeller efficiency. 
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Figure 16 
Total efficiency of propulsion systems 






Figure 17 


CONTROLS 

In controlling a man-powered aircraft the feet, 
which normally move the rudder, cannot be used for 
this purpose. In the Haessler-Villinger aircraft, move- 
ments of the control wheel were normal for longi- 
tudinal and lateral control. For moving the rudder, 
the wheel rotated about a vertical axis (Figure 17). 
After some training this arrangement did not cause 
any trouble to the pilot. 


To save weight, we did not use the conventional 
control surfaces for lateral and longitudinal control. 
For lateral control, the angle of incidence of the two 
halves of the wing could be changed differentially in 
opposite directions and, for longitudinal control, the 
angle of incidence was changed symmetrically. The 
stabilizer was fixed. 


This arrangement was used in the first year of our 
experiments. Lateral control was satisfactory as long as 
the aircraft did not get into extreme positions; there 
is a danger of unsymmetrical stall and I would not 
use this device again. Longitudinal control was not 
satisfactory; it caused over-controlled takeoffs, with 
too steep climbs and consequent loss of speed, and 
resulted in some hard landings and repairs. In the 
second year we changed back to a normal stabilizer 
and elevator arrangement with far better results. At 
the same time the size of the rudder was increased. 





Figure 18 
Control surfaces — 1935 
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Figure 19 
Control surfaces — 1936 


The difference in the tailplane design and rudder ar- 
rangement in 1935 and 1936 can be seen in Figures 18 
and 19. 


The use of spoilers for lateral control, as used on 
the Bossi-Bonomi aircraft, cannot be recommended, 
since the response of the aircraft to control move- 
ments is too slow. 


As a result of these experiences, I think that the 
control system of a man-powered aircraft should be 
completely conventional. Any deviation will cause de- 
lays arising from incorrect takeoffs, damage from hard 
landings, changes found necessary from flying ex- 
perience and so on. 


WEIGHT OF THE HAESSLER-VILLINGER AIRCRAFT 

Many designers assume that the natural way to 
achieve minimum weight in a man-powered aircraft 
would be to reduce the strength — which seems to be 
justifiable in this special case, since no flying would 
be done at heights greater than 50 ft. When we de- 
signed the Haessler-Villinger aircraft, we started out 
with an ultimate load factor of 2.5. A check of the 
wing deflections, however, showed that the strength 
in bending had to be increased until deflections could 
be considered normal. In this way we ended up with 
an ultimate load factor of about 6. 

We placed a restriction on diving speed to limit 
the torsion loads. Takeoff and landing loads could 
not be reduced. Our aircraft was only a quarter of 
the weight of a similar sailplane but still had almost 
full strength (Figures 20 and 21). 


The rough breakdown in Table 2 indicates how 
the weight goes down when a sailplane design is con- 
verted to suit the conditions of man-powered flight. 

It is important to stress every part of the aircraft 
carefully and keep a strict weight control during de- 
sign and actual building (Figures 22 and 23). 

Some pounds can be saved by selecting materials 
of lower specific weight and gauges below usual limits. 
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Figure 20 
Weight comparison 
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Figure 21 


Instead of using birch plywood with a minimum gauge 
of 0.8 mm (0.032 in), on many parts we used cedar 
plywood with lower specific weight and a minimum 
gauge of 0.6 mm (0.024 in). For covering the wing 
and tail surfaces, we used silk with a special type of 
impregnation. The weight of the fabric covering after 
impregnating was 160 gr/sq m, compared with the 
normal 300 gr/sq m on sailplanes. Special care has to 
be taken with the use of paints, lacquer and dope. 
No brush should be used, only spraying by a skilled 
painter. 


With reduction of cross-sections to a minimum the 
usual methods of gluing cannot be used. Special sup- 
ports have to be made for use during the gluing pro- 
cess, so that clamping pressures can be safely applied, 
and the weight of these supports and clamps must it- 
self be supported, since it cannot be carried by the 
airplane structure. 


The fuselage was covered with 0.6 mm cedar ply- 
wood to give a sufficiently stiff structure. Plywood 
gauges around the pilot’s seat and skidbox were 
normal. At first we thought to save weight with a 
fabric-covered fuselage, but calculations of deflections 
soon showed the advantages of the plywood covering 
without sacrificing weight. 


TABLE 2 
SCHEME FOR REDUCTION OF A/C WEIGHT TO MINIMUM 


Sailplane single-seater with cantilever 

wing of 13.5m span 150 kg 
1) Remove instruments 10 kg 
2) Reduce payload from 110 to 65 kg 

(Diff. payload = 41% at a rate of 


0.8 = 33%) 46 kg 
3) Convert cantilever wing to one with 
wire bracing 19 kg 
4) Reduce wing area to min weight/ 
span 10 kg 
5) Limit diving speed 5 kg 
6) Limit ability to withstand handling 
for parts with low stress 15 kg 105 kg 
45 kg 
7) Co-ordinate all factors in original 
min layout 10 kg = 


Weight of man-powered aircraft 35 kg 





Figure 22 


Aerodynamic efficiency 

In his lecture mentioned above’, Mr. Shenstone is 
reported to have criticized the Haessler-Villinger air- 
craft as aerodynamically inefficient. Since our aircraft 
with 13.5 m span had a glide ratio of 1:24, better than 
other sailplanes of the same span, and a sinking speed 
of 0.52 m/sec, which is below that of any comparable 
sailplane even today, I cannot see the grounds for 
this assertion. Probably the unusual wing design, 
with wire bracing, does not fit the normal concept 
of an aerodynamically efficient aircraft. The drag 
of the bracing wires was 3.5% of the whole drag, 
demanding 3.5% of the whole power. With 0.9 hp 
total, this amounted to 0.031 hp. This power re- 
quirement could be offset by a weight saving of 
6.2 lb, since the power requirement to keep 1 lb aloft 
with this aircraft was 0.005 hp, and the actual weight 
saving with our wing design, in comparison with 
a similar cantilever wing, was of the order of 30 lb. 
Thus we effected a saving of 17% of the available 
power and demonstrated the superiority of a braced 
wing over a cantilever wing. 

Our aspect ratio was 18.8. In modern sailplanes the 
limit for the aspect ratio of cantilever wings is about 
20. Man-powered aircraft will have to go to aspect 
ratios of 30 and over and this will not be feasible with 
any cantilever wing; we shall have to resort to other 
arrangements. 





Figure 23 
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Table 3 gives the weights for the main components 
of the Haessler-Villinger aircraft, which will be used 
as a basis for the calculations that follow. 


LIST OF SYMBOLS 

P = Required power — hp 

W .= Weight — kg or lb 

Via = Sinking speed — m/sec or ft/sec 
Efficiency of propulsion system 
Span — m or ft 

Aspect ratio 

Cy = Drag coefficient 
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Ci. = Lift coefficient 

D, = Induced drag 

A = Lift 

Q = Dynamic pressure 

S = Wing area — sq m or sq ft 


POWER, WEIGHT AND SPAN 


After taking stock, I will now try to demonstrate 
what we can expect today of a new man-powered 
aircraft. Since in man-powered flight the available 
power is a minimum, we will have-to search for an 
aircraft where the required power becomes a mini- 
mum. 


TABLE 3 
WEIGHT BREAKDOWN OF THE HAESSLER-VILLINGER AIRCRAFT 




















1) Wing 
a) Main Spar 
Spar caps 2.910 kg 
Spar web 1.115 
Uprights & fillers 0.456 4.481 kg 
b) Rear Spar 
Spar cap 0.325 
Web 0.471 0.796 
c) D-Nose skin 4.850 4.850 
d) Ete. 
Trailing edge 0.465 
Tip edge 0.100 
Qty 70 Main ribs 50 gr 3.500 
Qty 60 Skin stiffeners 15 gr 0.900 
14.3 sqm Impregnated fabric 
covering 2.280 
Fittings 0.400 
Lacquer & glue 1.238 8.883 19.010 kg 
2) Fuselage 
Frame 1 0.045 
es 2 0.150 
ite 0.250 
. 4 0.593 
eae 0.327 
: 6 0.280 
- 7 0.220 
ae 0.180 
7 0.140 
. 10 0.120 
= il 0.110 
12 0.085 2.500 
Longerons 2.520 2.520 
Skid 0.600 
Skidbox—Skin 0.900 
Skidbox—Formers , 0.160 
Seat & Floor 0.800 2.460 
Outer-skin 0.6 mm plywood 2.800 2.800 10.280 kg 
3) Empennage 4.000 
4) Controls in fuselage 0.500 : 
5) Wire bracing 0.600 5.100 kg 
EMPTY WEIGHT OF AIRCRAFT 34.380 kg 
Empty weight 34.380 kg 
Propulsion system 2.000 kg 
Energy accumulator 9.700 kg 
(Rubber bungee = 7.00 kg) 
Weight of aircraft 46.080 kg 
Pilot 65.000 kg 
TOTAL WEIGHT 111.080 kg 
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Figure 24 
Drag breakdown — Haessler-Villinger aircraft 


The power required for horizontal flight is — 








W Vain W Vain 
= - - kg/sec or - 


ft Ib/sec 1 
759 _ oo. 


The efficiency of the propulsion system, , was dis- 

cussed earlier. We will now have a closer look at the 

factors which affect the sinking speed — 
Veins = 4V W/8? VX Cp,.,/CL"* 


or 29V W/b? VA Cp,,./ Ci" 


m/sec 
ft/sec (2) 
at sea level. 


W//b’ is known as the span loading and is the main 
factor affecting the sinking speed. A in the formula 
calls for low aspect ratios. And the last element, 
Cotor/Ci'", requires the total drag of the aircraft to be 
as small as possible and the lift near the stalling value. 

Since 4 = b*/S, one could write \/W/b’ V/A as 
VW/b’ Vb°/S ie. \/W/S, indicating the primary 
importance of low wing loading. This misled several 
designers of man-powered aircraft into building wings 
of tremendous areas, which resulted in the aircraft 
becoming too heavy (1930 Zaschka aircraft — 475 sq 
ft; Bossi-Bonomi — 226 sq ft; Haessler-Villinger — 
107 sq ft). 

This sinking speed formula contains about every 
possible factor affecting the performance of the air- 
craft and the examination of each to find the design 
with minimum power requirements is a very time 
consuming business. 


Many possible ways of achieving aerodynamic im- 
provements have been discussed in other papers on 
man-powered aircraft. 1 will try to approach the 
power problem in another way. We have to decide 
which factors have the greatest influence on perform- 
ance, so that we can concentrate on them; all other 
factors must be packed into one constant and the 
upper and lower limits of this constant considered in 
the light of existing experience. 

When we take a look at Figure 24 showing the 
relative sizes of the different components of the total 
drag of the Haessler-Villinger aircraft, we see that 
the induced drag is 57% of the total. Reduction of 
the induced drag is the greatest problem in designing 
a man-powered aircraft. The profile drag is of next 
importance, but presents no design problem provided 
we confine ourselves to airfoils already used on 
modern sailplanes, the drag coefficients of which are 
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well established. All other components of the total 
drag are so small that any improvement in this area 
will have a negligible effect on the performance of 
the aircraft. The induced drag is 


r Ww? 


= Ont? = kg (3) 


i Orb? 
Since the lift equals the weight, we see the im- 
portance of weight and span. 


Figure 25 shows the span loading against sinking 
speed for a random assortment of well known gliders. 
The upper and lower limits of these values give lines 
going to zero in the shape of a parabola. For these 
lines we can now write 


0.66 < 
me _ m/sec (4) 
0.59 Ve 


This formula will give reliable results provided 
the design approximates to the conventional shape of 
gliders. Final performance calculations will have to 
be done in the usual way. 


Voink 


Using this formula, Eq. (1) giving the power 


necessary for horizontal flight can be written 
< 0.0089 wy W 1 cWVW 
0.0078 5 7 n Ob 


(For the Haessler-Villinger c = 0.0089 : for Proposal 
H4 c = 0.0078 metric). 

This formula shows that the main problem is to 
design a wing of maximum span and minimum weight. 
All other points are of secondary importance. That 
means that designing a man-powered aircraft is 
strictly a structural problem and not an aerodynamic 
one. Of course, the aerodynamic design has to be of 
a high order but it must not interfere with the weight 
problem. It is not necessary to break fresh ground 


hp = (5) 
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in aerodynamics: on the other hand, on the structural 
side, new design techniques will have to be used. 


WING WITH MINIMUM WEIGHT 

Now we will try to find the wing having minimum 
weight, for a given span and a given load. I will use 
the actual weights of the Haessler-Villinger aircraft as 
a basis and calculate the changes in weight for a range 
of wing areas from 4 to 28 sq m (40 to 300 sq ft) and 
for three different spans. In Figure 26 the weight of 
the bending material is plotted against wing area for 
spans of 46, 62, and 79 ft. For a given airfoil section, 
increase in wing area results in corresponding increase 
in wing thickness and consequently the weight of the 
bending material goes down; and vice versa. 


In the graph showing the weight of the spar web, 
we see the influence of span. Theoretically, for a 
given span and a given load, the weight of the web 
should not change with changing spar depth (i.e. wing 
thickness), but in practice changes arise from limit 
gauges of materials and the changing weight of the 
necessary stiffening. 

In the graph showing the weight of the D-nose, 
which handles the torsion loads, the influence of limit 
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Figure 27 
Wing weight against area for a single seater 


100 


weienT 





WING AREA 


Figure 28 
Wing weight against area for a two seater 


gauges and increasing stiffening requirements with in- 
creasing wing area is also apparent. 


The total weight for all other parts of the wing of 
the Haessler-Villinger aircraft, such as ribs, trailing 
edge, fabric covering etc, amounted to 0.21 lb/sq ft 
(1 kg/sq m) of wing area, and this figure is assumed 
to be unaffected by changing area. 

When we add the weight of these four structural 
groups together we can plot the weight of the 
wing against wing area for the three chosen spans 
(Figure 27). Each curve has a clearly marked mini- 
mum, which lies on a straight line going to zero. This 
gives the theoretical minimum weight of wing at each 
span. For a practical design, I would choose a point 
at a greater wing area, say, on the dotted line. The 
weight difference would be small and, since the aspect 
ratio of the theoretical minimum is very high and ac- 
cordingly gives high deflections, a wing at the dotted 
line, with lower aspect ratio, would be stiffer. The 





Figure 29 
Weight of wing, fuselage, tailsurface and power drive 
for a single seater 
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Figure 30 
Weight of wing, fuselage, tailsurface and power drive 
for a two seater 


best compromise, to locate this point for an actual 
wing, can only be determined by calculation of de- 
flections. A wing of 14 m span would certainly be 
satisfactory and the Haessler-Villinger wing had such 
dimensions. 


In Figure 28, I have repeated this process for the 
higher wing loading of a two seater aircraft. In this 
case, the theoretical minimum weights occur at 
slightly greater wing areas than in the first case, but 
not in proportion to the increased loads. The minimum 
weight of the two seater wing of 24 m span exceeds 
the minimum weight of the single seater wing by 
only 65%. 


AIRCRAFT WITH MINIMUM WEIGHT AND 
ITS PERFORMANCE 


Having found the optimum weight for a wing of 
given span, we can now draw a graph combining these 
results with the weights of the other components of 
the aircraft (Figures 29 and 30). The weights of the 
Haessler-Villinger aircraft will again serve as a basis, 
at 13.5 m span. The weight of the tail surfaces goes 
up in the same way as the wing. The propulsion sys- 
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Figure 31 
Minimum empty weight against payload for spans of 
40 ft to 100 ft (wing with bracing wires) 
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tem is not affected by the span and stays constant. 
The length of the fuselage changes with increasing 
span and this adds slightly to its weight. Adding up 
the weights of these components we can get, as the 
final result, the minimum weight of the whole aircraft 
for any given span, for the single seater and the two 
seater. For the 14 m span in each case, the two seater 
is heavier by 85%; for 30 m. span a two seater would 
be heavier by only 71% of the single seater. 

I have shown the above results in Figure 31 and 
drawn lines indicating the empty weight against pay- 
load for spans from 40 to 100 ft; these will be useful 
later for checking other cases. 
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Figure 32 


Required man-power against span 


Now we can calculate by Eq. (1) the required 
power for horizontal flight at any given weight and 
span. The results for both the single and two seater 
are shown in Figure 32. For better comparison be- 
tween single and two seater, I have added a curve 
showing the 50% power requirement of the two 
seater, which should correspond to the full power 
requirement of the single seater if performances were 
equal. This graph indicates clearly the superiority of 
the single seater, since in a two seater the power re- 
quirement for each man is 30% higher. This is again 
in conflict with Mr. Shenstone who is reported’ to 
have said that the Haessler-Villinger aircraft failed to 
achieve better results, because being a single seater 
the optimum power was not available — the optimum 
size of crew is considerably more than one. 
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REQUIRED POWER (HP) 





GROSS WEIGHT 
Figure 33 
Required power against gross weight for a one 
_. and two seater 


The point is borne out by Figure 33, which shows 
the relationship between required powér and gross 
weight. Since the required power is proportional to 
the 1.5th power of the weight, only the lightest single 
seater. structure has scope for future development. 
This does not mean that flying a two seater would 
not be possible. By using short-time power, flights 
could be achieved in a two seater up to 1.5 minutes, 
even without takeoff assistance. Under the same con- 
ditions, long-time flights of 1 hour could be made in 
a single seater, by using ground effect, as indicated in 
the graph. 

Since the level of the required power is directly 
related to the time for which this power is available, 
we are now able to plot possible flying time against 
gross weight (Figure 34). As a basis for this calcula- 
tion I have used the power of a professional racing 
cyclist weighing 65 kg, having a cruising power of 
0.46 hp and reserve energy of 1100 m kg. Although 
the Haessler-Villinger aircraft is shown on the graph 
as having a possible flying time of 2 minutes with 
ground effect — when we actually achieved only 1 
minute — it must be remembered that we did not 
actually use a pilot of this power level. If the aircraft 
could have been built without the weight penalty of 
the energy accumulator, the possible flying time with 
the professional cyclist would have gone up to 5 
minutes. On the other hand, for the Bossi-Bonomi 
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aircraft on the same power basis we get a possible 
horizontal flying time of 38 seconds with ground 
effect (assuming the same efficiency of the power 
drive). 

It would be helpful in giving a clearer picture of 
the actual quality of an aircraft if other authors would 
use the same “standard pilot” as a basis for comparison. 

The Haessler-Villinger aircraft, without the 
weight of the energy accumulator, had a performance 
similar to our calculated aircraft of 46 ft span. The 
aircraft of 62 ft span has a possible flying time of 
1.5 minutes without ground effect and could make 
long-time flights of over 1 hour with ground effect 
at a height of 10 ft. The aircraft of 79 ft span would 
be capable of a flight of 2.2 minutes without ground 
effect, and a flight of over 1 hour with ground effect 
at heights of more than 10 ft. 


INFLUENCE OF PILOT WEIGHT ON 
AIRCRAFT PERFORMANCE 

These calculations have been based on the standard 
pilot with a weight of 65 kg. Now we will check 
what influence a change in the weight of the pilot 
has on the performance of the aircraft. The empty 
weight of the aircraft differs according to the payload 
for which it is designed. We will assume that, for 
each pilot weight, an aircraft with the corresponding 
empty weight is available, giving us the line relating 
gross weight to pilot weight in Figure 35. 

In Figure 36 the required power, with and without 
ground effect, plotted against gross weight and pilot 
weight, is compared with various available powers. 
We see that the feasibility of long-time horizontal 
flights is increased by selecting pilots of lower weights. 
If it were possible to find a professional racing cyclist 
weighing 120 lb, long-time flights without ground 
effect would be possible with aircraft weights calcu- 
lated on the basis of that of the Haessler-Villinger. 

Since the cruising power of national record break- 
ing amateurs would be 10% higher, and of world 
record breaking professionals 30% higher than that of 
the average professional cyclist shown in the graph, 
we see that the chosen power level has little hope of 
attainment: the average amateur cyclist has no hope 
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Figure 34 
Performance with professional cyclist — 65 kg 
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Figure 35 
Pilot weight against aircraft weight for a span of 79 ft 


of achieving long-time horizontal flights. The power 
which was available for my tests with the pilots 
Duennebeil and Hofmann would be insufficient even 
today, with better aircraft, for flights of over 5 
minutes. Bossi was lucky in finding a pilot with the 
right power level, and with flying experience at the 
same time. 


PROPOSAL 

Having reviewed possibilities of designing a suc- 
cessful man-powered aircraft, we should now select 
a performance goal which would justify the necessary 
work and expense. Then the aircraft should be de- 
signed for this performance, leaving a good: margin 
for unforeseen losses, to minimize the risk of failure. 
Performance for restricted flights 

This goal should be a flight of 1 hour, at cruising 
power, at a height equal to the span i.e. just out of 
ground effect. Takeoff should be unassisted, since at 
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Feasibility of long-time flight 
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Figure 37 
Man-powered aircraft — Proposal Haessler H4 


ground level the ground effect helps to bridge the 
power gap between takeoff and ceiling height. 

The lowest acceptable performance would be a 
flight of 1 hour, at cruising power with the ground 
effect, at a height of 10 ft. Takeoff could be done 
with limited assistance. The aircraft would have to be 
designed in such a way that either a skid or a power- 
driven wheel could be used. Testing would begin with 
the skid and takeoff assistance, and, when the per- 
formance has increased with training and growing 
experience — which could take a month or a year — 
to the point at which the pilot is able to get near to 
the ceiling height, the skid should be replaced by the 
power-driven wheel and the pilot can take off by 
himself. 

If the weight of the aircraft cannot be reduced 
sufficiently or the pilot fails to come up to specifica- 
tion or both, then discarding the weight of the 
wheel, saving the pilot the effort of takeoff and allow- 
ing him to fly in the ground effect still makes it pos- 
sible to accomplish a 1 hour flight, with these 
limitations. 

An alternative lowest acceptable performance 
would be a takeoff using the power-driven wheel, the 
pilot using up his energy reserves in a short-time flight 
of less than 5 minutes. This could be accomplished 
at the same weight and pilot-power as the restricted 
1 hour flight with ground effect. 


Performance for unrestricted flights 

If we want to fly an aircraft without the above- 
mentioned restrictions, the structure must conform to 
normal strength requirements for a load comprising 
the average pilot, parachute and instruments. How- 
ever, the application of man-power might still in- 
fluence the performance sufficiently to justify the 
attempt. If we assume that a glide ratio of 1:75, with 
man-power, could be considered as the lower limit 
for unrestricted flight, by leaving out instruments and 
parachute and using a professional cyclist, such an 
aircraft would still be able to fly, with ground effect, 
for 5 minutes or more in restricted flight. 

It is not within the scope of this paper to speculate 
how a combination of soaring and man-power could 
improve performance. Experience will show the actual 
advantages and disadvantages of such a combination. 
Only the feasibility will be checked here. 

The increase in weight to meet normal strength 
requirements would be small, since only the weight 
of the D-nose will go up for greater torsion loads. 
Since the weight level achieved with the Haessler- 
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Villinger aircraft cannot be considered as an absolute 
minimum, further improvements may well be possible, 
and we will make the assumption that, as a limit case, 
any increase in the weight of the D-nose can be offset 
by other weight savings. In the following Proposal H4, 
I use the already calculated weights as a basis, and 
show the percentage by which the weight could go 
up in certain areas. without materially affecting the 
performance of the aircraft. 


Specifications for Proposal H4 
The dimensions of the aircraft (Figure 37) are 
given in Table 4. 


TABLE 4 
ProposaL H4 





Span | 25.00 m | 82 ft 

Wing area | 16.75 sq m | 180 sq ft 

Wing loading | 6.6-12.0 kg/sqm | 1.35-2.44 Ib/sq ft 
Aspect ratio | 37.3 37.3 

Chord | 0.67 m 26.4 in 

Glide ratio 1:44 1:44 

Fuselage cross-section | 0.35 sq m 3.77 sq ft 

Prop dia 1.70 m 67 in 

Prop eff. | 0.84 0.84 


Wing section | NACA 633-618 


i 








An aircraft with these dimensions can be built for 
different payloads at correspondingly different empty 
weights. Three possible weight cases, A, B, and C, 
are shown in Figures 38 and 39. 


Case A would be an aircraft for restricted 
flights with man-power only. Pilot weight would 
be between 132 and 141 lb and power requirements 
call for an average professional racing cyclist or 
top notch amateur. Empty weight should be be- 
tween 110 and 152 Ib (10% below and 17% above 
calculated weight). Flights could be made for 1 
hour, up to heights of 60 ft at the lower weight 
limit and at 10 ft at the upper weight limit. 

Case B is an aircraft for restricted and un- 
restricted flight. Without parachute and instru- 
ments a professional cyclist could fly this aircraft 
for 1 hour at a height of 10 ft. 

An average amateur racing cyclist, with a 
weight between 132 and 145 lb, could fly un- 
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Figure 39 
Performance of Proposal H4 (span of 82 ft) 
with cases A, B and C 


restricted with parachute and instruments and 
improve the glide ratio to 1:90 with his cruising 
power and to 1:100 with 5 minute power. The 
empty weight of 186 lb can be 17% above the 
calculated weight. 

Case C is the upper weight limit for the air- 
craft, where a glide ratio of 1:75 can still be 
achieved in unrestricted flight at the 5 minute 
power level of the average amateur. The highest 
acceptable weight of the empty aircraft could be 
238 lb, which would be 40% above calculated 
weight. With 11 Ib less the glide ratio of 1:75 
could be held with cruising power for 1 hour. 


We see that the wide range of acceptable perform- 
ances allows sufficiently wide weight tolerances to 
enable an aircraft to be built without too much diffi- 
culty. The main factor in these wide limits is the 
full strength of the aircraft. In building man-powered 
aircraft reduction of strength does not give the hoped 
for weight reductions and limits the use of the aircraft. 
If the man-powered aircraft is to have any future, 
it will probably have to be built to normal strength 
requirements. 
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FALLING PROBE ELECTRON DENSITY MEASUREMENTS? 


by Professor P. A. Forsyth* and A. Kavadas** 


Institute of Upper Atmospheric Physics, University of Saskatchewan 


Ov of the more fundamental measurements that 

can be carried out in the ionosphere is a deterinina- 
tion of the electron content. For that part of the iono- 
sphere that is accessible to ground based sounders, the 
vertical electron density profile of the undisturbed 
ionosphere is fairly well known. For many years 
regular measurements have been carried out by iono- 
spheric sounders at stations located all over the world. 
Before the advent of the research rocket there remained 
certain ambiguities in the sounder results, but these 
have been resolved by means of the rocket measure- 
ments of Seddon, Jackson arid others. The method used 
in these early rocket experiments is shown in Figure 1, 
and depends upon the differential retardation suffered 
by radio waves of different frequencies when passing 
through an ionized medium. The usual method in- 
volves the use of two harmonically related frequencies, 
one near the critical frequency of the region to be 
investigated, the other being the sixth or seventh 
harmonic of the first. These two frequencies are 
transmitted simultaneously from the rocket as-it pro- 
ceeds upward into the ionosphere. Due to the retarda- 
tion suffered by the lower frequency, the received 
signals on the ground change their relative phases. 





Figure 1 


Illustrating the principle of the propagation experiment of 
Seddon and Jackson for measuring ionospheric 
electron densities. 





tPaper read at the Canadian High Altitude Research Sympo- 
sium in Ottawa on the 21st October, 1960. 
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Figure 2 
A typical electron density profile for the undisturbed 
ionosphere. The broken line suggests a modification 
that might occur during an auroral disturbance. 


The phase comparison is made on the ground by first 
producing the appropriate harmonic of the lower of 
the two received signals and then comparing this 
harmonic with the signal received at the higher fre- 
quency. By measuring the relative phase change of 
the two signals as the rocket goes into the ionosphere, 
it is possible to deduce the total electron content of 
the ionosphere between the rocket and the ground 
station. It is the variation with height of this Seca 
density, integrated along the line of sight from the 
ground station to the rocket, that provides the distri- 
bution of electron density as a function of height in 
the ionosphere. Clearly the sensitivity of this method 
to local variations is very limited. A typical electron 
density profile for the undisturbed ionosphere is shown 
in Figure 2. 


When the ionosphere is disturbed the electron 
density can be expected to vary considerably from 
point to point, and under these conditions the propaga- 
tional type of experiment, using either a ground based 
or a rocket borne transmitter, will be subject to con- 
siderable error. This was recognized several years ago 
and various attempts have been made to overcome the 
difficulty. The problem becomes most severe in the 
case of auroral measurements, because it is believed 
that the ionization associated with an auroral display 
is very inhomogeneous and that the electron density 
varies through wide limits in distances of the order of 
a few meters. The radio evidence concerning this 
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Figure 3 
Illustrating the principle of operation of the RF impedance 
probe method of measuring the electron density. 


ionization, obtained by means of propagation experi- 
ments, is still incomplete and poorly understood. A 
rocket experiment that would measure the electron 
density in relatively small volumes of the ionosphere 
during an auroral display would be of great value. It 
would also enhance considerably the value of the 
information if it could be obtained simultaneously for 
two or three separate trajectories in the ionosphere. 
This additional information would help to distinguish 
between the temporal variations and the spatial varia- 
tions. 


In an attempt to study the various methods that 
could be used for electron density measurements in 
the disturbed or auroral ionosphere, a NASA team, 
consisting of Bourdeau, Jackson, Kane and Serbu, 
flew a specially instrumented rocket from Churchill 
in September 1959. This rocket carried a conventional 
propagation experiment and, in addition, an RF probe 
experiment. The principle of the RF probe is shown in 
Figure 3. The probe consists essentially of 2n antenna 
driven by a relatively weak source of constant RF 
energy. The impedance of this antenna changes when 
it is immersed in an ionized medium, and this change 
in impedance can be used to measure the electron 
density in the surrounding medium. This type of 
experiment gives great promise for measurements that 
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Figure 4 a 


Comparison of the results obtained in a single rocket using 
both the RF impedance probe and the propagational 
method of measurement. (After Bourdeau, Jackson, Kane 
and Serbu, International Space Science Symposium, 
Nice, 1960.) 
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must be made in inhomogeneous media such as the 
auroral atmosphere. Figure 4 shows a comparison of 
the electron density results from the two experiments 
in the same rocket (the propagation experiment and 
the impedance probe experiment), from which it may 
be seen that the propagation experiment gave a value 
which was some 1.5 times as great as that obtained 
from the impedance probe experiment near the peak 
of the trajectory. Because this flight was a test of the 
two methods, it was flown during a period of quiet 
ionospheric conditions, and so it is believed that the 
propagation experiment gave the best results. Fortuna- 
tely, the experimental team had also put into the same 
rocket a Langmuir probe and an electric meter, from 
which it was possible to deduce that the rocket ac- 
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f/f 
Figure 5 
Illustrating the variation of impedance (Z) of an antenna 
in an ionized medium, as a function of electron density 
and frequency (f). Z, is the free space antenna impedance 
and f, is the critical frequency of the medium. 


quired a negative potential in the ionosphere and that 
the RF probe measurements were, therefore, in error, 
because of an ion sheath that formed about the probe 
antenna. The effect of this positive ion sheath of a few 
centimeters thickness forming about the RF probe 
antenna is to decrease the electron density in the 
immediate vicinity of the probe and therefore to de- 
crease the effect of the ionized medium on the impe- 
dance of the probe. In spite of this sheath formation, 
the RF probe method is likely to be of great value in 
investigations of the inhomogeneities in the ionosphere. 


With this as background, we are now describing 
another electron density experiment which we hope 
to try next year. The experiment is a modification of 
the RF impedance probe, and is intended to investigate 
the structure of the ionization in an auroral display. 
The basic difference in the new experiment is that the 
RF source is not fixed in frequency. The frequency 
is made variable and caused to track the critical fre- 
quency of the surrounding medium continuously. It 
is hoped in this way to overcome the difficulties intro- 
duced by sheath formation. Figure 5 shows the varia- 
tions of impedance of an RF probe as a function of 
electron density (f.") and frequency (f). Just at the 
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Figure 6 


Illustrating the approximate trajectories of the three 
ejected packages. 


critical frequency (f/f. = 1) a large change of impe- 
dance takes place and, in particular, the radiation 
resistance of the antenna becomes infinite. The so- 
called critical frequency at which these sudden changes 
take place is determined by the electron density 
according to the expression f° = Ne’/xm, where N, e 
and m are the electronic density, charge and mass, 
respectively. The effective volume surrounding the 
probe has a radius of the order of a wavelength, and 
the measurement should not be significantly affected 
by the presence of a very thin layer surrounding the 
probe, particularly if this layer is one that is deficient 
in electrons. 


It has already been mentioned that it would be 
desirable in studying the auroral ionization to make 
the measurement at two or three points simultaneously. 
In our experiment we plan to eject simultaneously 
from the rocket three separate instrumented packages, 
each of which will carry out continuous measurements 
of the local electron density along its particular trajec- 
tory. A sketch of these trajectories is shown in Figure 
6. The ejection of a small instrumented package from 
the parent vehicle also has the very considerable 
advantage of removing the measuring device from the 
contaminating influence of the rocket body and its 
surroundings. 





Figure 7 
A sketch of the rocket nose cone showing the location of 
the packages before ejection. 
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Figure 8 


A sketch of the ejected package showing the 
internal construction. 


Figure 7 shows a sketch of the rocket nose cone 
with the three packages in position. As soon as the 
rocket is well clear of the lower atmosphere, these 
three packages are simultaneously ejected approxima- 
tely in the horizontal plane but in directions separated 
by 120° in azimuth. They then travel on their indivi- 
dual trajectories which continue to separate farther 
and farther apart until impact. Each package is, of 
course, equipped with its own telemetry transmitter so 
that the three sets of measurements can be continuous- 
ly monitored from the ground. 


Figure 8 shows a sketch of the ejected package. 
The four short antennas are for the telemetry trans- 
mitter, but also serve as a ground plane for the long 
antenna which is the sensing probe. Internally, the 
package consists only of the telemetry transmitter, a 
suitable modulator, the sensing circuit and battery 
pack. In addition, of course, facilities have to be pro- 
vided for the monitoring of the operation of the 
internal circuits up until the time of ring. 


Figure 9 shows the measuring or sensing circuit. 
The measurement is performed with the aid of a vari- 





Figure 9 
Block diagram of the probe instrumentation. 
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Figure 10 
Graph of the variation of probing frequency with bridge 
unbalance voltage. 


able frequency oscillator and an impedance bridge, 
which remains approximately balanced as long as the 
radiation resistance of the antenna maintains a value 
close to that experienced in free space. The external 
sensing antenna is in one of the arms of the impedance 
bridge and in the other is a network having the same 
impedance variation with frequency as that of the 
antenna. Thus, as long as the antenna is in free space, 
the bridge is balanced for all frequencies and the 
oscillator seeks its lowest permitted frequency. If the 
local electron density increases so that the plasma fre- 
quency (or local critical frequency) exceeds the oscil- 
lator frequency, then the sudden sharp change in the 
antenna impedance will unbalance the bridge and 


cause the oscillator frequency to increase until it 
exceeds the local plasma frequency and the bridge 
balance is again restored. This feed-back action keeps 
the oscillator frequency just above the plasma frequen- 
cy. The actual frequency controlling element is a 
varicap, a device which changes its capacity in response 
to an applied de potential. In order to monitor the 
local electron density it is only necessary to monitor 
the frequency of oscillation of the sensing oscillator. 
This is done by a separate circuit and the resulting 
signal is applied to the telemetry transmitter. In Figure 
10 is shown the variation of frequency of the oscillator 
as a function of applied voltage. The choice of a fre- 
quency range for the experiment is not easy because 
so little is known about the auroral ionization, but we 
have finally settled on that shown in Figure 10. 


CONCLUSION 


This experiment is the first endeavour in the rocket 
field undertaken by our University group. As a result, 
we are learning as we go along and there have already 
been several false starts. We have completely revised 
the package layout and shape several times, but we feel 
that we are now close to the final design layout. We 
have not yet solved all of the problems connected with 
the ejections of the packages from the rocket although 
we do have the design of an ejection package which 
we hope we will be able to fit into the rocket nose 
cone. In any case, the project is going forward quite 
rapidly and we hope that all our problems will be 
solved in time to permit the firing of the rocket early 
next spring. 
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AERODYNAMICS OF BLASTSt 


by Dr. I. I. Glass,* A.F.C.A.I. 


Institute of Aerophysics, University of Toronto 


SUMMARY 


For the past decade theoretical and experimental investiga- 
tions have been conducted at the Institute of Aerophysics on 
shock tube flows. Recently, this work was extended to include 
the analogous flows generated by spherical and cylindrical ex- 
plosions. Glass spheres and cylinders are pressurized by means 
of compressed gases or combustible mixtures and are shattered 
to generate an explosion. Similér methods can be used for 
implosions, underwater explosions and wave interactions. The 
finite mass, strength and breaking time of the glass diaphragms 
impose some limitations on certain experiments, where this 
technique is used to generate a blast wave. Nevertheless, the 
method has pene very valuable in the study of many basic 
properties of spherical and cylindrical blast phenomena that 
have been investigated, using piezo-pressure gauges and several 
schlieren and shadowgraph techniques. Some consideration is 
given to intense explosions from concentrated energy sources 
for spherical, cylindrical and planar blasts, and explosions and 
implosions generated from finite sources with the same geome- 
tries. Blast wave simulators for aerodynamic tests and the 
dynamic testing of structural components are briefly discussed. 


LIST OF SYMBOLS 
a sound speed (ft/sec, mm/y sec) 
a(y) function of y 

specific heat at constant volume 

(ft Ib/slug°R, cal/gm°C) 
C, specific heat at constant pressure (as above) 
E. total energy in spherical flows (ft lb, 
equivalent gm of TNT) 
energy per unit length i in cylindrical flows 
energy per unit area in planar flows 
E is proportional to E., E = E./a(y) 
E./p:, an energy parameter (ft’, ft’, ft 
depending on v) 
WE./p, ft, for spherical flows 
flow quantity at a given radial distance r or R 
G(r/R) a function of the reduced radial 
distance 
C,/C, specific heat ratio 
r/e non-dimensional energy reduced distance 


« B 


a 


~~” 


7? 


(S/a,) shock Mach number 
number of grid zones in shock transition 


non-dimensional distance 


meg” 





+Paper read at the Joint 1.A.S./C.A.I. Meeting in Montreal on 
the 18th October, 1960. 
*Professor of Aeronautical Engineering 
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n a power index 

v 1, 2, 3 for planar, cylindrical, spherical flows 

p pressure (lb/ft’, dynes/cm’, or atm) 

AP. peak overpressure at the shock front (atm) 

7 3.1416, also non-dimensional parameters 

q artificial viscosity [non-dimensional form 
Eq. (17)] 

R shock wave radius (in, ft) 

r radial distance (in, ft) 

r Lagrangian or initial radius (in, ft) 

To diaphragm radius (in) 

p density (slugs/ft’, gm/cm") 

S shock velocity (ft/sec, mm/y sec) 

t time (sec) 

~ a, t/e non-dimensional energy reduced time 

T temperature (°K) 

eS t 

TFT Ep an = non-dimensional time 

u et velocity (ft/sec) 

Ui free stream velocity in steady flow (ft/sec) 

V volume (ft’, cm’) 

X % (r:/e)*, a non-dimensional Lagrangian - 
variable in spherical flows 

AX _ grid spacing in numerical computation 

Subscripts 

0 a reference atmosphere 14.7 psi, 0°C 

1 ambient or atmospheric conditions 

2 conditions immediately behind the shock wave 

3 conditions immediately behind the rarefaction 
wave 

+ conditions in the diaphragm 

Pa p:/p. shock pressure ratio 


Pa p./p, diaphragm pressure ratio 
Tn T:/T; shock temperature ratio 
Ta T./T, diaphragm temperature ratio 


SS forward and backward facing shock waves 

(particles enter from right and left, 

respectively) 

RR ___ forward and backward facing rarefaction 
waves 

Cc contact surface 

Air/Air diaphragm separating interior air at high 
pressure (and temperature) from exterior 
air at low pressure; similarly for HE/Air 
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INTRODUCTION 


B LAST waves are familiar phenomena. They occur 

in nature in a relatively mild form during thunder- 
storms or they may become very intense during 
volcanic (steam) explosions. It is estimated that the 
Krakatao explosion of 1883 resulted in an energy re- 
lease equivalent to 5000 megatons of TNT’. However, 
man was also able to generate controlled or uncon- 
trolled shock waves of increasing intensity when he 
learned how to crack a whip’, discovered propellants 
and explosives, made use of coal mines and grain 
elevators and, finally, when he invented nuclear 
weapons’. 


The aerodynamics of blasts is not only of interest 
in itself, but it is also closely related to supersonic 
and hypersonic flight, and its now familiar shock wave 
patterns and sonic booms. It is of interest to note, for 
example, that a pound mass, at the earth’s escape speed 
of 7 miles/sec, possesses a kinetic energy equivalent 


to over 14 lb of TNT (Figure 1). Consequently, a _ 


ballistic re-entry at such speeds is very much like a 
powerful explosion. It is fortunate that most of this 
energy is absorbed through the bow shock in the 
denser layers of the atmosphere, by heating the air 
which is left behind in the wake of the vehicle. Only 
a very small percentage of this energy is finally trans- 
ferred to the vehicle itself, otherwise it would quickly 
vaporize. In a similar manner, the finite amount of 
energy in an explosion spends itself through shock 
wave heating of an ever increasing mass of air until 
the shock becomes vanishingly weak. 


Mechanical, chemical, electrical or nuclear blasts 
all dissipate their energy through shock waves. In the 
latter, however, about half of the energy is dissipated 
through thermal and nuclear radiation, although at 
very early times the thermal energy leakage is small’. 
Large flow velocities, temperatures and pressures are 
induced behind a shock wave and it is the pressure 
and drag forces that cause the structural damage. The 
flow quantities decay in value with time and distance, 
and when the shock wave is far enough away from 
the point of explosion it has attenuated to a sound 
wave, and the physical quantities behind the wave are 
then only vanishingly different from those in the am- 
bient atmosphere. A similar type of decay applies to 
the shock system surrounding a supersonic aircraft so 
that when it flies at a sufficiently high altitude its 
effects on the ground are negligibly small. 


The properties of both idealized intense explosions 
from concentrated constant energy sources and very 
weak explosions can be predicted analytically* * °. On 
the other hand, the analysis of actual explosions from 
finite sources must be treated by numerical methods” *. 
Such analyses have only become possible in recent 
years with the advent of high speed electronic com- 


_ puters and the increased knowledge of the thermo- 


dynamic properties of explosives and of the air itself. 


Some of the analytical solutions for intense ex- 
plosions can also be applied to steady hypersonic flight 
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Figure 1 


Variation of kinetic energy per pound mass 
(TNT equivalent) with re-entry velocity 


through the so-called blast wave theory. For example, 
the solution for the intense cylindrical blast wave 
carries over to steady hypersonic flow over axisym- 
metric blunt-nosed slender bodies, and the planar blast 
to two-dimensional flows” * "”. 


The method of generating cylindrical and spherical 
explosions, implosions and wave interactions that has 
been developed at UTIA, consists of bursting glass 
spheres or cylinders from 1 to 6 inches diameter. For 
explosion studies these glass “diaphragms” are pres- 
surized by means of compressed gases or combustible 
mixtures and exploded in a facility which consists of 
a 3 ft diameter steel sphere whose ambient atmosphere 
can also be varied. For implosion experiments the steel 
sphere is pressurized and the glass spheres and cylinders 
are evacuated and implode owing to the existing pres- 
sure differential or they can be broken mechanically. 
In the case of an underwater blast the sphere is ex- 
ploded through overpressure in a small tank filled with 
water. A few preliminary experiments of exploding 
spheres in a supersonic flow at M = 2.5 have also been 
tried®. The head-on collision of relatively weak 
spherical shock waves has been successfully investi- 


gated”. 


An explosion generated from a 2 inch diameter 
glass sphere pressurized to 400 psi is driven by a blast 
energy equivalent to about 50 milligrams of TNT. 
Although the scaling of strong explosions in a facility 
of this type is aerodynamically quite possible it would 
only apply to very small models and at relatively 
small overpressures (up to several atmospheres). Con- 
sequently, shock tubes, sector tubes and other devices 
that are driven by TNT, or other chemical explosives, 
must be used to produce high pressure pulses of long 
duration (fractions of a second), which are required 
for the proper simulation of the flow over large 
models with relatively long structural response times. 
The present facility is best suited to basic studies of 
the gasdynamics or the aerodynamics of the blast it- 
self or the flows induced over small models by 
explosions. 
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INTENSE EXPLOSIONS FROM CONCENTRATED 
ENERGY SOURCES — SPHERICAL, CYLINDRICAL 
AND PLANAR BLASTS 


The instantaneous release of a finite amount of 
energy in a concentrated form at a point, line or plane 
gives rise to an intense blast wave in a spherical, 
cylindrical or one-dimensional geometry, respectively. 
All of these blast waves possess the property that the 
profiles of any thermodynamic and dynamic flow 
quantity f (pressure, temperature, density and flow 
velocity) relative to the corresponding value im- 
mediately behind the shock front remain invariant 
with time t, when the distances 7 are expressed as 
fractions of the shock wave radius from the point of 
the explosion, r/R. That is, as the decaying Dlast 
wave engulfs an ever greater mass of gas, at a given 
fractional distance, the physical quantities possess the 
same ratio, or f(7r,t)/f(R,t) = G(r/R). This simi- 
larity is predicted on the assumption that the shock 
wave is so strong that the ambient pressure and tem- 
perature are negligibly small in comparison with those 
behind the shock wave. Consequently, the similarity 
solution applies to the intermediate region between 
the origin of the explosion, where p>0 and T>o , 
and that radius of the shock front where the shock 
wave has attenuated to an extent that the strong-shock 
assumption is invalid. 


This type of problem was considered independ- 
ently and almost simultaneously for the spherical case 
by Taylor‘, Sedov’ and von Neumann’. The cylindri- 
cal case was considered later by Lin”, who also noted 
that his solution should be applicable to the analysis 
of a steady flow over an axisymmetric, blunt-nosed 
slender body at zero angle of attack at hypersonic 
speeds (where r is equivalent to y and t is equivalent 
to x/U,)". An observer viewing a passing hypersonic 
body through a slit normal to the flight velocity will 
see very significant changes in the radial flow velocity 
as the body appears to drive the shock wave in a 
piston-like manner. Consequently, the flow in the 
(y, x)-plane appears like a cylindrical explosion in the 
(r,t)-plane. The flow in transverse planes exhibits a 
self-similarity in both cases when the shock wave is 
very intense and can be treated as having one inde- 
pendent variable. In a hypersonic axisymmetric flow 
the total drag is related to the energy per unit length 
in a cylindrical explosion. The planar, cylindrical and 
spherical blast waves were treated by Sedov’ and later 
by Sakurai”. 


In solving this type of problem it may be seen that 
the solution will depend on the following six quanti- 
ties: y, pi, p, Eo, r and t. The first three quantities 
describe the gas, E. is the instantaneous energy release, 
and r and t are the independent variables in the space- 
time plane of interest. From dimensional analysis 
(using length, mass and time as the three basic quan- 
tities), it is known that five minus three or two non- 
dimensional (2) functions can be formed that will 
describe the subsequent motion (y is already dimen- 
sionless). In the spherical case, for example, these 
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assumed that for the strong blast wave p: > p: ~ 0, 
the non-dimensional time parameter .7 becomes zero 
and the motion is independent of time and is self- 
similar as noted above, and only two-dimensional con- 
stants E, and p: are essential for a solution of the 
problem. 

If we assume that the gas is perfect and inviscid, 
the equations of continuity, motion and energy (par- 
ticle isentropic) which are valid in the region bounded 
by the blast wave are given by: 


Op _ Apu) pu 

at ap TES DS «6 (1) 
Ou Ou 1090p 

— +u—+ =0 (2) 


Ot Or por 


od (p d (p 
3i(5) +05-(5) =0 (3) 


where, v = 1, 2 and 3 for planar, cylindrical and 
spherical flows, respectively. 

In general, this set of nonlinear, partial differential 
equations is too difficult to solve and recourse must 
be made to numerical techniques. However, in the 
case of intense blasts where only two dimensional 
quantities E. and p; affect the motion, the two inde- 
pendent variables r and t are so related that the equa- 
tions can be reduced to ordinary differential equations 
and exact self-similar solutions are possible’. The 
boundary conditions that have to be satisfied are that 
the particle velocity at the origin is zero and that the 


_ following Rankine-Hugoniot equations are satisfied 


across the shock wave at a radius equal to R: 
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For intense blasts, Ms is very large and the above re- 
lations reduce to 


2 

ke = y +1 S (4a) 
+1 

p2 = ee pi (Sa) 
2 

pb: = a Ss (6a) 


The solution of Eqs. (1) to (3), subject to the 
above boundary conditions, are given in detail in Re- 
ference 5 for planar, cylindrical and spherical flows, 
and the results are illustrated in Figure 2. It is seen 
that the pressure behind the shock wave is a maximum 
and falls off quite rapidly near the shock wave to level 
off to a nearly constant value for radial positions 
r/R<0.5. As expected, the fall-off in pressure is 
greatest for the spherical case, with its greatest free- 
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Self similar profiles for intense planar (—-—), 
cylindrical (——-——) and spherical blast 
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dom for expansion, and least for the planar case. The 
same tendency is even more accentuated in the density 
profiles, where it is seen that nearly all of the mass of 
gas engulfed by the blast is concentrated close to the 
shock front itself. The pressure and density profiles 
are reflected in the temperature profiles where it is 
seen that, subject to the assumption of an inviscid 
non-conducting gas, enormous temperatures are de- 
veloped towards the centre of the blast, as a result of 
the prevailing vanishingly small densities and finite 
pressures. This is particularly marked in the spherical 
case. The particle velocities decrease from their maxi- 
mum value behind the shock front to zero at the 
origin of the blast. It is seen that the curves differ only 
slightly in the range 0.5 <7r/R < 1, but here at a given 
r/R, the spherical flow velocity is the lowest and the 
planar velocity the largest. 

The motion of the shock wave is, of course, also 
found from the solution. However, an insight into this 
can be obtained from the dimensional considerations 
noted above. Since E. represents the total energy re- 
leased in the spherical case, per unit length in the 


TABLE 1 
SoME PROPERTIES OF INTENSE BLAST WAVES 


Geometry |R, (Eq. (8)) S = dR/dt 


E V6 1/6 1/2 
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cylindrical case and per unit area in the planar case, 
it has dimensions of ML’* T*, where M, L and T 
represent the usual basic quantities of mass, length and 
time. Consequently, the non-dimensional distance 
noted above can be expressed as 


r 


ao" oer (7) 
¢ ae f2/ (2+v) 
Pi 
where E = E./a(y) 


The function a(y) is a constant and is determined 
from the solution of the equations of motion. With 
the above notation the shock wave path can be ob- 
tained from Eq. (7) as 


1/(2+y) 
R= (=) f2/(2+r) (8) 
Pi 


where, for convenience, Y was set equal to unity. 
Using Eqs. (8) and (6a), the quantities in Table 1 
can be derived. From.the shock radius R, shock speed 
S, the pressure behind the shock p, or the overpressure 
AP., it can be seen that the spherical blast wave decays 
most rapidly with time or shock radius. 


It is worthwhile considering some actual applica- 
tions of the above equations. Eq. (8) was used by 
Taylor* to deduce the energy yield of the atomic ex- 
plosion in New Mexico. Measurements of shock radii 
with time from high speed photographs gave the 
points shown in Figure 3, which is a plot of Eq. (8) 
for the spherical case; that is 


5/2 log R(cm) — log t(sec) = 11.915 (9) 


The value 11.915 is the intercept of the 45° line with 
the ordinate axis at log t = 0, and equals % log E/p:. 
For p= 1.25 X 10° gm/cm"’, E = 6.76 X 10” X p; ergs 
= 8.45 X 10” ergs. For y = 1.4, the solution for the 
spherical blast yields a(y) = 0.851 or E. = 7.19 X 10” 
ergs, or, using a value of the energy content of 
TNT as 1000 cal/gm, this is equivalent to a yield of 





log t (sec) 


Figure 3 
Velocity of the spherical blast wave from 
an intense explosion‘ 
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Figure 4 
Shock path from an exploding wire™ 


17 kilotons of TNT. It should be noted that an up- 
ward translation of the straight line in Figure 3 would 
be indicative of a stronger or higher energy yield 
spherical blast and vice versa. The answer obtained in 
this manner is only approximate considering the sim- 
plifying assumptions of a point source, constant 
energy and a perfect inviscid gas. In an actual case 
a significant portion of the energy is dissipated 
through thermal radiation at later times®. Nevertheless, 
the method shows that a great deal of insight can be 
obtained from this type of analysis. 


An example of an intense cylindrical blast wave is 
the sudden explosion of a fine cylindrical metal wire, 
caused by the sudden input of a heavy electrical cur- 
rent pulse. The analysis of the motion of the shock 
wave is made more difficult by the complex manner 
in which energy is added, and the change of phase 
of the wire from a solid to a vapour that occurs at 
very small times. In addition, the imperfect equations 
of state for the metal and the ambient gas into which 
the blast propagates must be taken into account. In 
fact, the exploding wire becomes analogous to a hot, 
pressurized cylinder that is suddenly ruptured, rather 
than a line source. Consequently, imploding shock 
waves and wave interactions occur that invalidate the 
similarity assumptions. 

A recent survey of exploding wires is given in 
Reference 15, where Bennett’s experimental results in- 
dicate a variation favouring Rat’, rather than ¢** as 
given by Eq. (8). For this case a(y) =~ 1 (see Refer- 
ences 5 and 10). One of the best results for a 5 mil 





(a) 


Figure 5 
Blast wave produced in a T tube 
(a) schlieren record of the (r, t)-plane 
(b) explanatory sketch: air at 500 mm Hg 
and 300K°K, capacitors 2yf at 17,5 ky 
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exploding copper wire with a rated energy of 60 
joules/em (~ 15 milligram TNT/cm) is shown in 
Figure 4. Here, the shock wave is. advancing into air 
at 0.5 atm ata rate of R « t’’, However, in the range 
of % to 1 atm, a single plot of R « t** appeared to fit 
the remaining results even better. The fact that all of 
these experiments at different ambient densities gave 
an apparently constant E/p:, or that a decreasing p 
could apparently decrease the energy added to the 
copper wire is surprising. However, when this result 
is considered from the viewpoint of the cylindrical 
analogue to the shock tube noted above, it might be 
explained by the fact that a diaphragm pressure ratio 
that changes by a factor of 5 (as in this case) would 
hardly affect the strength of a strong shock wave™. 
Some critical comments regarding this problem by 
Bennett and Rouse can be found in Reference 15. 


(= 0.2 mm.Hg.) 
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Figure 6 
Shock wave paths in a T tube in air and argon 


A planar type of blast wave may be generated by 
the sudden addition of electrical energy to a gas by 
means of an electrode or electrodeless discharge. The 
so-called T tube, developed by Fowler and Kolb, 
is a typical example (see Reference 16 for a brief sur- 
vey). Such a device for adding energy to a gas in a 
plane (ideally), as investigated at UTIA, is shown 
schematically in an insert in Figure 5. Energy during 
a discharge is added to the gas across a gap from two 
1.5 microfarad, 30 kv, low inductance capacitors 
by means of ohmic heating and through the magneto- 
gasdynamic ponderomotive force applied as a conse- 
quence of the rigid back-strap. During an experiment 
the current, rate of current and voltage inputs are 
monitored across the gap by means of oscillograms. 
The energy input is therefore known as a function of 
time. The sudden addition of electrical energy coupled 
with the ponderomotive force generates a blast wave 
in the 22 mm pyrex tube, whose path is illustrated in 
Figure 5 in an (7, t)-plane schlieren record. 
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The initial conditions are noted in the figure, and 
from the oscillograms a computation has indicated that 
about 20, 8 and 5 joules/cm* were added to the gas in 
three damped half-cycles at about 3 usec intervals. 
The latter two probably did not contribute very much 
to driving the blast wave as there was no apparent 
strengthening of the shock wave after the first half- 
cycle. Figure 5 shows the hot, luminous expanding gas 
driving the shock wave in front of it. During the first 
25 usec it is difficult to separate the shock front from 
the luminous gas on the schlieren record. Some typical 
shock paths are shown in Figure 6 as solid lines and 
the dashed lines show the planar blast wave decay 
given by Eq. (8). In this instance a(y) ~ 1.08, where 
y = 1.4, and a(y) ~ 0.60, where y = 5/3”. Also in this 
case, similarity is not achieved during the early phase 
of the flow because the assumption of energy addition 
at a plane is not met. A finite source is probably closer 
to reality. At higher pressures (5 mm Hg) and after 
about 20 usec, R « t’* appears to be a reasonable ap- 
proximation for subsequent small time intervals. Be- 
low 20 usec the shock waves are weaker than pre- 
dicted by this power law. However, for the case of 
air at a pressure of 0.2 mm Hg, a power law profile 
appears to be entirely inadequate. The blast wave is 
considerably stronger at the low pressures, which 
would result from a higher effective diaphragm pres- 
sure ratio. Further details are given in Reference 17. 


It is worthwhile to note at this point that, in the 
analysis of an intense blast, at times when shock wave 
attenuation no longer permits the assumption that the 
pressure in front of the shock wave may be neglected 
by comparison with that behind it, the exact Rankine- 
Hugoniot shock conditions must be satisfied and that 
the solution is no longer self-similar. Eqs. (1) to (3) 
must now be solved numerically. This has been done 
by Brode’, who made use of the concept of an artificial 
viscosity in order to overcome the difficulty of flow 
discontinuities at the shock wave. The artificial vis- 
cosity spreads the shock front over about 6 mesh 
lengths and a continuous transition results, analogous 
to an actual shock transition. The results are shown 
in Figure 7 for air considered as a perfect gas and an 
imperfect gas. It is seen that for shock wave peak 
overpressures, in atmospheres, in the range AP, > 10 
atm, AP, « 1/K* (a slope of 72°) and at very low over- 
pressures, AP. — 1/)., the acoustic result for a spheri- 
cal wave (a slope of 45°). In the intermediate range 
AP, can be fitted with a polynomial in 1/A,”. 


The perfect gas solution lies above the imperfect 
gas case and is indicative of the energy that is bound 
in vibration, dissociation, electronic excitation and 
ionization, and is not available for driving the blast 
wave. Brode estimates that the imperfect point-source 
explosion is only about 50% as efficient as that in 
perfect air (y = 1.4)*. 

Brode has reported three explosions from very hot 
spheres of gas, where the interior density is assumed 
as that at standard conditions but at pressures of 122, 
2000 and 20,000 atm” *. In the latter case the required 
uniform temperature is about 386,000°K. The 122 atm 
sphere approaches the perfect gas point-source solu- 
tion at AP, ~ 20 and A, ~ 0.2, and the 2000 atm sphere 
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Figure 7 
Variation of peak overpressures (AP,) with shock 
radius (},) for various explosions. Point source, air, 





perfect (— — — —), imperfect (— - — - —); hot 
sphere 2 < 10‘ atm and 3.86 < 10° °K ( 3 
TNT (—--—-- »’ 


when AP. ~ 300 and A, ~ 0.8. The 20,000 atm sphere 
approaches the imperfect gas solution when AP, ~ 2500 
atm and A, ~ 0.035, as shown in Figure 7. The latter 
then diverges somewhat from the imperfect point- 
source curve and is brought closer to it by the over- 
taking of the second shock and finally the third shock. 
After an overpressure of 10 atm the hot-sphere solu- 
tion runs almost parallel to the imperfect gas point- 
source solution. The circled points are the data from 
a nuclear explosion, as reported in Reference 3, and 
are seen to agree very. well with Brode’s calculations, 
except at low overpressures when the explosion 
energy is small compared with the ambient energy of 
the air engulfed by the blast wave, and the accuracy 
in the calculation of maintaining a constant blast 
energy is reduced. 

The TNT curve on Figure 7 exhibits a decay in AP, 
~.1/hs initially”. Later, it cuts across the other three 
curves, and at low overpressures lies closest to the im- 
perfect gas solution. As Brode points out, every ex- 
plosion is affected by its initial history to some degree, 
and consequently they are difficult to compare or 
scale accurately. He indicates that a rough rule is that 
the peak overpressure in an explosion will decay in a 
manner similar to a point-source when the shock wave 
has engulfed a mass ten times the mass of the initial 
explosive. Furthermore, although the peak overpres- 
sure may be scalable the flow quantities within the 
blast interior may differ for each source, and identical 
scaling becomes very difficult”. 

Three profiles of overpressure to peak overpres- 
sure ratio AP/AP., density ratio p/p: and particle 
velocity u/u. from the blast centre to the shock front 
for the point-source, perfect gas explosion have 
been drawn on Figure 7, for + = 0.00147, 0.1436 and 
0.745. These have been included to illustrate how, in 
time, the blast wave profiles depart from the similar 
solution for the strong shock wave“, and develop 
Positive and negative phases when the overpressure 
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drops. A hot-spot remains near the origin even at late 
times, as indicated by the persistent near-zero densities 
and finite pressures. 

An important conclusion that one might draw is 
that a finite source explosion will eventually decay 
very much like a point-source blast. The effects of the 
wave interactions appear to vanish quickly if the 
source is hot (high sound speed). Since the cube of 
the inverse radius type of shock-decay relation applies 
to the spherical-shock overpressures long after the 


strong shock conditions really cease to be accurate, the 


corresponding attenuation relations would also be ex- 
pected to apply to the cylindrical and planar blasts. 
However, this is not entirely borne out by the ex- 
ploding wire or T tube experiments noted above. In 
these cases it is possible that a modified type of simi- 
larity solution based on a decreasing energy input with 
time might yield a more realistic type of power pro- 
file for the shock path at early times. Later on the 
decay may well be represented by that from a finite 
source explosion in a manner exhibited by the nuclear 
blast points in Figure 7. 


EXPLOSIONS AND IMPLOSIONS GENERATED 
FROM FINITE SOURCES 

The foregoing results indicate that a realistic 
approach to the propagation of blast waves would be 
to consider an energy source of finite extent for 
spherical, cylindrical and planar explosions. Once a 
finite source is assumed then the complementary prob- 
lem of implosions may be considered simultaneously. 
In order to simplify the analysis it is assumed that no 
energy is lost through radiation and that the flow is 
inviscid. The problem is then to determine the result- 
ing flow properties when a plane, cylindrical or 
spherical diaphragm, which separates a gas at higher 
or lower pressure (and temperature) from that of the 
external atmosphere, is suddenly ruptured generating 
an explosion or an implosion, respectively (Figure 8). 
In the case of a planar flow one assumes a very small 
length of chamber r. relative to the channel in the 
case of an explosion. On the other hand, in the case 
of an implosion the channel is made very short. 


For very small times, less than the time it takes for 
the head of the rarefaction wave or the shock wave 
to reach the origin in the case of an explosion and 
implosion, respectively, Eqs. (1) to (3) can be solved 
by using a power series solution or the method of 
characteristics” ”. The results are illustrated schema- 
tically in Figure 8. The planar case is the familiar 
shock tube case, which has been treated in great detail 
by many authors”. Here state (4) is at high pressure 
(and temperature) and state (1) at low pressure. Since 
the flow is symmetrical about the origin, consider that 
the right diaphragm is ruptured. Immediately a shock 
front S moves into region (1), compressing and 
heating the gas in region (2) in an irreversible manner. 
At the same instant a rarefaction wave R moves 
into state (4) and isentropically expands and ac- 
celerates the particles in the opposite direction to form 
state (3). This gives rise to a contact surface C, 
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Figure 8 
Spherical, cylindrical and planar explosions and implosions 


which separates the uniform state (2) existing at a 
high entropy level from that of state (3). 


For the cylindrical and spherical explosions the 
picture is changed. States (2) and (3) are no longer 
uniform. The shocks and contact surfaces decelerate. 
Consequently, state (2) is a region of decreasing 
entropy behind the shock path. The characteristic 
lines in the expansion fan accelerate towards the shock 
front”, and it will be shown that the existence of an 
imploding shock wave turns state (3) into a non- 
isentropic region. These effects are even more pro- 
nounced in the spherical case than in the cylindrical 
case (Figure 8). 

When an implosion occurs, the shock accelerates 
towards the origin and it is followed by an accelerat- 
ing contact surface. The characteristic lines in the 
expansion wave decelerate in a direction opposite to 
the shock front. State (2) is now a region of in- 
creasing entropy behind the shock path. The effects 
are again accentuated for the spherical flow. As ex- 
pected, for the same diaphragm pressure ratio Pu, 
a finite source explosion generates a stronger rare- 
faction wave for the spherical and cylindrical blasts 
than for the planar explosion, and vice versa for the 
implosion. Furthermore, it might have been antici- 
pated that, since in cylindrical and spherical flows the 
contact surface decelerates and it also behaves as a 
massless piston, it would send expansion pulses to 
decay the main shock wave and compression pulses to 
form a second shock at the tail of the rarefaction wave, 
in the case of cylindrical and spherical explosions. In 
implosions the reverse takes place. Here the compres- 
sion pulses reinforce the main shock and the expansion 
pulses reinforce the main rarefaction. wave. Conse- 
quently, a second shock wave does not occur in an 
implosion. This will be verified subsequently. 


Although the above approach for small times. is 
of interest, it does not solve the problem. A solution 
has to be obtained from a numerical integration of a 
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set of stable, finite difference equations that approxi- 
mate the nonlinear partial differential equations (Eqs. 
(1) to (3)). The boundary conditions of zero 
velocity at the origin and the Rankine-Hugoniot re- 
lations across shock waves must also be satisfied. It 
might be noted that although a solution is possible by 
using the method of characteristics for the planar case, 
it is complicated by the additional term in the con- 
tinuity equation for the cylindrical and spherical cases. 
This term leads to a singularity at the origin and causes 
the planar Riemann invariants [24/(y—1) + u] to vary 
along characteristic lines so that the method becomes 
complex, especially when unique values for the flow 
quantities no longer exist as a result of the formation 
of shock waves from steepening compression waves. 


The method of numerically integrating Eqs. (1) 
to (3) has been developed by Brode’ and successfully 
applied to a number of problems*” *. For conveni- 
ence, Eqs. (1) to (3) are written in Lagrangian form”. 
The independent variables are time t and the initial 
radial position of a gas element 7. The dependent 
thermodynamic and dynamic variables are made non- 
dimensional by using the ambient gas properties (for 
example, p = p/p:, “ = u/a, etc, and the bars can be 
omitted for simplicity). Similarily, the independent 
variables are made non-dimensional by using the pres- 
sure-reduced energy of the gas contained in the in- 
terior of the diaphragm. 


For example, 


VpsC,T. P. 
BL Oi Se at Pa (9) 
pi pr ¥o— 1 
where 
V = 4/3 7r,* spherical flows, vy = 3 
V = ar? per unit length, cylindrical flows, vy = 2 
V = 2r, per unit area, planar flows, vy = 1 


For spherical flows in particular, 


é es Py 
4/3mr2 Yys—1 





(10) 


From Eqs. (9) and (10) it is seen that the reduced 
energy per unit volume (reduced energy density) de- 
pends only on the diaphragm pressure ratio Pa and 
the specific heat ratio y. in all blast cases, and ap- 
proaches infinitely large values when Pa > © or 
y.: > 1. The usefulness of Eq. (10) will be considered 
subsequently. The following dimensionless quantities 
are now formed: 

A =r/e 


(11) 
T = at/e 


With the above notation, the Lagrangian equations of 
motion in non-dimensional form can be expressed for 
the spherically symmetric case, as follows’: 


ae en 12 
OX pd (12) 

Ou »* = 
ay (2 +9) (13) 


ie Se. 
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Figure 9 
Computed wave system of a blast from a 2 inch diameter 


sphere case Air/Air, p, = 22 atm, p, = 1 atm, T, = T, 
= 299°K (shock tube flow reference: S = 2110 fps, 
u, = 1250 fps, T, = 470°K, T, = 180°K) 


0 10 
oP am Pip + (y — 104) (14) 
Or por 
Or 
“% = ia (15) 


1 /r\3 
x =-(— (16) 
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The quantity X is a Lagrangian variable expressing 
the non-dimensional volume per steradian for a spheri- 
cal shell of radius 1. The quantity q is the dimension- 
less form of the artificial viscosity that provides the 
means of automatically spreading the shock wave into 
a thin, continuous transition”. Brode has found that 
the following form of q produces the required results 
without significantly distorting the rest of the flow 


field: 
ate ‘ Ou 
4 ie re) max ay | ae roe - 


Where, AX is the grid size (~ 1/100 distance traversed 
by the shock wave), N is the number of grid zones 
in the shock front (~ 6) and @u/@X is negative for 


shock waves so that ae( 5 =) for shocks and gq = 


for expansion waves or for contact surfaces. Eqs. i 
to (17) provide five equations for the five unknowns 
1, P, p, & and q in terms of X and +. A discussion of 
the finite difference equations that have been used by 


Brode to approximate the above partial differential | 


equations and the required stability conditions in the 
time interval Ar can be found in Reference 23. 


Spherical explosions and implosions 
Explosions have been produced by bursting 1 inch 
diameter to 5 inch diameter pressurized glass spheres. 
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Compressed air, helium, argon, sulfurhexafluoride and 
combustible mixtures of oxygen-hydrogen diluted 
with helium were used as driver gases. However, ex- 
tensive numerical solutions by Brode were done only 
for the following two cases: 2 inch diameter spheres 
containing air at 22 atm and helium at 18% atm and 
300°K. These cases were checked experimentally in 
some detail” *. It should perhaps be emphasized again 
that the shock-tube problem only has a simple analyti- 
cal solution as long as the diaphragm separates cham- 
ber and channel sections of sufficient length to elimi- 


_ mate immediate wave interactions. The amount of 


work required to obtain a solution for a single spheri- 
cal or cylindrical blast problem is rather extensive, 
and it is for this reason that the investigation was 
limited to the above two cases. Under the conditions 
noted above, the air and helium spheres contained an 
energy equivalent to 88 and 44 milligrams of TNT 
(1040 cal/gm), respectively. The nominal 50 mm 
diameter glass spheres were about 1 mm thick and 
weighed approximately 20 gm. The spheres were 
blown individually of pyrex or soda lime glass with 
stems about 65 mm long, 6 mm diameter and 1 mm 
thick. 

The original computations, as previously described, 
for the Air/Air and He/Air blasts were done in 
the non-dimensional radius-time (A, r)-plane and on 
several graphs have been converted for convenience 
to units used in the experimental work. The con- 
version factors for the (7,t)-plane are: 4 = 0.2 cor- 
responds to r = 1.1 inches, r = 0.2 corresponds to 
t= 90 usec for Air/Air, and 0.2 4 = 0.94 inches and 
0.2 r = 72 wsec for He/Air. A standard atmosphere 
of 14.7 psi and 0°C (subscript 0) was used by Brode 
as a reference atmosphere rather than the ambient 
conditions. 


The paths of the various wave elements are shown 
in the (7, t)-plane in Figure 9 for the case of Air/Air 
and in Figure 10 for He/Air. It is seen that when the 
sphere is ruptured, a decelerating shock wave S, 
moves into (1) followed by a decelerating contact 
front C. A rarefaction wave R_ propagates into 
(4) and a second shock wave S:, which forms at 
the tail of the rarefaction wave but faces in the op- 
posite direction to S:, is swept outward by the ex- 
panding high pressure gas. However, as the second 
shock wave gains strength, it overcomes this flow 
velocity and finally implodes on the origin. Refrac- 
tions at the contact front produce secondary waves 
(rarefaction and shock waves in the case of the air 
explosion, but only shocks for the case of the helium 
blast) and cause it to oscillate as indicated, and after 
some time several shock waves S,, S:, Ss move away 
from the explosion. Unlike similarly-facing shock 
waves in the plane case these waves do not overtake” 
at the times and distances considered here. For ex- 
ample, S. is caught in the negative-pressure phase of 
the primary shock wave where it encounters a low 
sound speed and opposing particle velocity that slows 
it up. In the case of the hot sphere, shown in Figure 7, 
S. does not encounter such conditions and does over- 
take to strengthen S,, and this is repeated again, by a 
third shock S,. Unlike the planar flow, regions (2) 
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Figure 10 
Computed wave system of a blast from a 2 inch diameter 
sphere case He/Air, p, = 18144 atm, p, = 1 atm, T, = T, 
= 301°K (shock tube flow reference: S = 2700 fps, 
u, = 1850 fps, T, = 600°K, T, = 197°K) 


and (3) are not uniform or isentropic owing to the 
existence of these shock waves of continuously vary- 
ing strength. 

It has been shown by Brode , that various types 
of blasts approach and have a decay rate like the limit- 
ing point-source solution for shock overpressure 
versus radius, after the primary shock S, has engulfed 
a mass of air of about tenfold that of the original 
explosive. For a hot or light gas this occurs much more 
rapidly than for a heavy gas for a given initial volume. 
Therefore, for the same size sphere r. and dia- 
phragm pressure ratio Pu, the combination He/Air 
gives initially a much stronger shock than say, SF./Air 
(shock tube effect)". However, SF. is a more efficient 
explosive when the point-source solution is finally 
reached and the scaling laws can be applied approxi- 
mately (energy loading effect — Eq. (10)). That is, 
the same shock overpressure would take place at a 
larger physical radius for the case of SFs. These effects 
are illustrated in Figure 11 for the cases of He/Air 
and Air/Air blasts used in the present experiments. It 
is seen that the He/Air case approaches the point- 
source solution faster, owing to the lower mass of 
helium, or since the sound speed in helium is much 
greater the effects produced by the wave interactions 
vanish sooner. Initially, the helium gives a larger peak 
overpressure AP, at a given radius R than the air 
case. However, for the same initial diaphragm pressure 
ratio P. the air sphere contains 5/3 as much energy 
as the helium sphere and, if both explosions would 
coincide with that of the point-source solution, AP, 
at a given R would be larger for the air blast than 
for the helium blast. From Figure 11 it is seen that the 
initial history of a given blast is not “forgotten” and 
no curve is coincident with any other. Consequently, 
identical scaling using different sources is not possible, 
although as noted above for many purposes the close 
proximity of the curves may be quite sufficient for 
scaling. 
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Figure 11 
Shock overpressure (AP,) for. various explosions. Point 
source, perfect gas (— — — —), imperfect gas (- ---); 
TNT, 1.5 gm/cm*, (—--—); He/Air, 18% atm, 
(— - —); Air/Air, 22 atm, ( ) 





It is worth noting from Eq. (10) that once a 
spherical flow solution has been found for a given 
P. and y. (this fixes the type of explosive), then it 
does not matter if r. or « (mass of explosive) is in- 
creased or decreased. The same solution applies, as 
this only increases or decreases 4 and the explosion 
is magnified or reduced in the physical (r, t)-plane. 

The Air/Air and He/Air explosions approach the 
point-source curve at positions when AP. decays in a 
more complex manner than AP, ~ | A," since the shock 
front is neither strong nor weak. The helium blast 
approximates the point-source, imperfect air attenua- 
tion when the shock is at A. ~ 0.5 or R ~ 2.4 inches, 
at a shock overpressure AP, ~ 1.8 atm. If one con- 
sidered a 20 kiloton point-source blast, then applying 
the energy scaling relation such an overpressure would 
exist at a distance R ~ (20 X 10° X 10° X 10°)” 
(44)°% (2.4/12) = 1540 ft or 0.29 miles. The Air/Air 
blast cuts the point-source imperfect air solution when 
A. ~ 1.1 or R ~ 6 inches and AP, ~ 0.35 atm, but it 
has a more gradual attenuation. The equivalent scaled 
distance for a 20 kiloton blast when R is at the same 





Figure 12 
Lines of constant particle velocity in the (), 7)-plane. 
Case He/Air (see Figure 10) 
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overpressure occurs at R ~ (2 X 10/88)" x 6/12 
= 3050 ft or 0.58 miles. Consequently, only shock 
wave effects at large distances from very intense ex- 
plosions with correspondingly small models could be 
simulated by using the present shock sphere apparatus. 

In addition to the wave traces in the (r,t)-plane 
it is also possible to plot all the thermodynamic and 
dynamic properties of the explosions (see Refer- 
ences 7, 8, 18, 21 and 24, for example) as functions 
of time or radial distance. Contours of constant par- 
ticle velocity and density are shown in the (A, r)- 
plane in Figures 12 and 13, respectively. The decrease 
in density behind the in¢ident and reflected rarefaction 
waves is quite apparent. The particle velocity however 
is increased through the incident rarefaction wave 
and decreased through the reflected rarefaction wave 
until it is zero at the origin. The density in front of 
the imploding shock wave near the origin is very 
small. The imploding shock also induces large negative 
velocities as it nears the origin. The effects are felt 





Figure 13 
Lines of constant density in the (), r)-plane. 


Case He/Air (see Figure 10) 


over the entire flow field and it is seen that a line of 
zero particle velocity exists behind the exploding 
primary shock wave. This is the beginning of the 
negative phase, and the pressure becomes subatmo- 
spheric to the left of this contour. It is of interest to 
note that the constant density lines in the (+, 4)-plane 
resemble their steady flow hypersonic counterpart in 
the (x, y)-plane, for a body whose axisymmetric shape 
is given approximately by the contact surface locus 
in the (t,7r)-plane (the blast wave analogy would 
really apply to the cylindrical explosion, since the 
spherical blast involves four dimensions and it has no 
analogue in steady flow). 


A second shock wave is also usually observed in 
the steady flow case. However, it always appears to 
be coupled with some overexpansion at the shoulder 
of a flat plate with a cylindrical nose, say, where the 
boundary layer thins and a recompression is neces- 
sary for the streamlines to match the body contour. 
Unlike an explosion, the second shock wave in super- 
sonic flow appears to be caused by viscous effects, 
but this has not been verified analytically. 
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Pressure profiles at various times after an Air/Air 
explosion (see Figure 9) 


As an example of a typical profile, the variation of 
pressure from the main shock front to the centre of 
the explosion is shown in Figure 14. (Any number of 
such profiles, or of density, temperature, particle 
velocity, dynamic pressure, and the static and dynamic 
pressure impulses, can be drawn from the detailed 
tabulated results.) The pressure profiles are shown for 
times of 130, 195, 310 and 420 msec. The contact 
surface is identified by a bar. At 130 usec the rare- 
faction wave has already reflected from the origin 
where the pressure is quite low. The second shock S; 
is imploding on the origin and the main shock S, has 
covered over 2 inches. As time goes on, the pressures 
behind the main shock and in front of the second 
shock are dropping, giving rise to a weaker first shock 
and a stronger second shock wave, respectively. At 
420 usec, S, has reflected from the origin and is very 
much stronger than S., which it is now following. 


The variation of temperature with time at the 
origin for the air blast is shown in Figure 15. It is seen 
that the temperature undergoes extreme changes from 
the reflected rarefaction wave to the reflected implod- 
ing shock wave (20°K to 3000°K). The origin is a 
singular point, and Guderley” has shown that the 
imploding shock wave in a perfect gas approaches 
infinite strength when it reaches the origin. For an 
imperfect gas (or perfect gas) considering viscosity 
and heat conduction this would hardly occur, further- 
more the continuum equations may not apply very 
close to the origin. However, this is not the reason 
why the temperature (or any other thermodynamic 
quantity) remains finite on the present graph. It is due 
to the artificial viscosity that spreads the shock wave 
over about 6 mesh points and consequently it is not 
possible to approach the origin to within a few 
hundredths of an inch. For all practical purposes this 
is quite sufficient. The experiments conducted at 
UTIA indicate a finite shock velocity at the origin, 
within the limits of resolution of the (r,t)-plane 
schlieren records for cylindrical and spherical flows. 

The question why the second or imploding shock 
wave occurs is often asked. It arises as a result of the 
spherical or cylindrical nature of the flow and can be 
understood by observing the path of the contact front, 
which indicates the manner in which the interior 


March, 1961 


high pressure gas expands to its nearly stationary value 
at a pressure somewhat below the ambient atmosphere. 
As this front decelerates, it sends out ahead of itself 
rarefaction pulses that overtake and decay the —_— 
shock. However, behind it compression pulses are 
sent out that steepen into the second shock wave. 
(The behaviour is qualitatively similar to the one- 
dimensional problem of impulsively accelerating a 
piston to a uniform velocity and then decelerating it 
instantaneously.) Such an effect does not occur in the 
planar case, where initially the contact front moves at 
a uniform velocity before it interacts with the shock 
or the rarefaction wave. 


It might be noted that in the case of an implosion, 
the motion of the contact front is such as to generate 
in front of it compression pulses or Mach waves to 
further increase the strength of the imploding shock 
wave. However, behind it rarefaction pulses are pro- 
duced. Consequently, a second shock cannot occur 
in an implosion. Furthermore, these rarefaction pulses 
make it possible to have a weaker initial rarefaction 
wave for a spherical implosion than in a shock tube 
flow. The inward motion of the contact front is 
stopped by the reflected imploding shock wave. 

Although explosion experiments have been made 
using spheres from 1 inch to 5 inch diameter, the 
2 inch diameter spheres have proven most satisfactory 
for the present 12 inch diameter field of view schlieren 
system. The 1 inch diameter sphere shrinks the phe- 
nomena at the origin while the 5 inch diameter magni- 
fies it too much in many cases. Consequently, most 
of the blast effects will be illustrated with the results 
from the 2 inch diameter spheres. 


Figure 16 shows a sequence of four shadowgrams 
that illustrate the bursting process of a 2 inch diameter 
glass sphere, pressurized at 300 psi, when it is hit with 
a Y% inch diameter mallet head and the subsequent 
development of the flow and the shock wave. The 
mallet head is attached to the slender rod which is part 
of a solenoid actuated breaker mechanism. At the to 
of the figures appears the rubber coupling that holds 
the stem through which the sphere is pressurized. The 
foil or wire at the right side forms part of the 
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Figure 15 
Temperature at the origin as a function of time after an 
Air/Air explosion (see Figure 9). 
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Figure 16 
Spark shadowgrams of the explosions generated from a 
22 inch diameter pressurized glass sphere. Case Air/Air 
Pp, = 300 psi, p, = 14.7 psi, T, = T, = 300°K, times after 
breaker struck the sphere, in y sec: (a) 200, (b) 255, 
(ce) 375, (d) 850 
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triggering circuit. It is seen that the gas rushes out 
through the fractured glass to generate compression 
waves which coalesce to form the shock wave. The 
contact region exhibits turbulence owing to its passage 
over the jagged glass fragments and its front is there- 
fore irregular. Although this has an effect on the 
shape of the shock front at small times, the wave soon 
stabilizes and becomes perfectly spherical as shown 
in Figure 16c. The mallet blow affects the collective 
shape of the glass fragments so that they assume an 
ellipsoidal form as in Figure 16d, where the appear- 
ance and diffraction of the second shock wave can 
also be seen” ”. The glass fragments are not small. It 
is only when they strike the wall of the steel sphere 
that they are pulverized. 


Figure 17 is a schlieren record of the bursting of 
a 2 inch diameter glass sphere taken with a high speed, 
multi-source spark camera™. In this case the sphere was 
ruptured through overpressure alone. A band of con- 
ducting paint circled the sphere and when the paint 
was severed it supplied the pulse for triggering the 
spark sources. The central frame, which is on the 
optical axis, is circular (1 inch diameter) and is used 
for alignment purposes only. The remaining 8 frames 
are on a 2% inch diameter pitch circle and undergo 
vignetting owing to the off-axis location of the spark 
light sources of the schlieren system”, whose 12 inch 
diameter beams also transluminate viewing windows 
of identical size. All the 9 frames are taken on standard 
4 inch X 5 inch Kodak Royal Pan cut film. The 
framing, sequence is from | to 8, at intervals noted in 
Figure 17. The growth of the explosion is illustrated 
very well in this series. Especially noteworthy is the 
second shock wave in frame 6, and its disappearance 
in frame 7. 


The same type of blast is shown in Figure 18, 
which is a schlieren photograph of the (7, t)-plane. 
The paths of the main shock wave, the second shock 
waves, the front of the contact region and the glass 
fragments are readily seen. From this and many other 
records it can be concluded that the imploding second 
shock and its reflection at the origin advance at finite 
speeds. The glass sphere breaks up in about 50 to 
60 mu sec and the fragments accelerate to velocities of 
over 200 ft/sec. As a result, the shock front does not 
form instantly but as a consequence of coalescing 
compression waves as noted above. This can be seen in 
Figure 18 at small times and is similar to the formation 
process of plane shock waves in a shock tube”. 


A comparison of theory and experiment for the 
case of the air blast is shown in Figure 19, Owing to 
the finite duration of the breaking process of the glass 
sphere, the actual primary shock takes some time to 
form and to achieve full strength, and consequently 
it falls behind the theoretical wave. The glass frag- 
ments are quickly accelerated to a uniform velocity 
(~ 250 ft/sec) and represents an equivalent loss in 
energy available to drive the explosion (~ 20%). 
These effects cause the second shock to implode about 
15 mwsec late. However, it races on and soon does 
match the theoretically predicted shock trace. The 
oscillations of the contact front are much milder than 
predicted by theory, and this no doubt arises from the 
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Figure 17 


High speed schlieren photographs of an explosion 
generated from a 2 inch diameter pressurized glass sphere. 
Case Air/Air, p, = 370 psi, p, = 14.5 psi, T, = 298°K, 

times in ,, sec from frame Nos. 1 to 8: 50, 100, 150, 
200, 250, 800, 900, 1000 





(a) (b) 


Figure 18 


Blast from a 2 inch diameter pressurized glass sphere 
(a) schlieren record of the (r, t)-plane, 
(b) explanatory sketch: G = glass sphere, 
F = glass fragments, S, = main shock, 
S, = second shock, C = contact surface. 
Case Air/Air, p, = 326 psi, p, = 14.4 psi, 
T, = T, = 297°K 
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Figure 19 


Comparison of theory and experiment of the wave system 

following the explosion from a pressurized glass sphere. 

Case Air/Air, p, = 326 psi, p, = 15 psi, T, = T, = 300°K. 

(Note that only the boundaries of the experimental results 
have been given for convenience) 


fact that the actual contact front becomes quite tur- 
bulent as it passes over the glass fragments. Conse- 
quently, the expected refraction process is more 
gradual. The large scale mixing at the contact front 
also causes the cold gas to occupy a larger (and 
asymmetrical) volume than predicted by theory. Some 
detailed discussions of the blast generated by an ex- 
ploding glass sphere can be found in Reference 24, 
and more recently in Reference 25. 


Two (r,t)-schlieren records for a He/Air ex- 
plosion appear in Figure 20, where S., S. and S, are 
clearly seen. The implosion of S. is masked by the 
glass fragments since it occurs sooner than in air. The 
excellent reproducibility for the two runs for the 
6 inch and 12 inch radii are shown in Figures 20a and 
20b and is quite typical of all runs.’ The experimental 
results for the helium blast appear in Figure 21, and 
are seen to be qualitatively similar to those described 
for the air case. Owing to the initial formation delay, 
the second experimental shock wave appears where the 
calculated third shock occurs. 


It is worth noting that the final results shown in 
Figures 19 and 21 were obtained in a very much im- 
proved apparatus”, as illustrated in Figure 22, whereas 
the original explosion runs were made in a simple 
wooden box”. The glass spheres and cylinders were 
housed in a spherical shell of 3 ft diameter that could 
be pressurized or evacuated for explosions, implosions 
or wave interaction studies. 


Some experiments were also conducted at the same 
diaphragm pressure ratio Pa. but at half the absolute 
pressure level in order to verify Eq. (10), and to simu- 
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(a) (b) 


Figure 20 
Schlieren records of the (r, t)-plane showing the blast from 
a 2 inch diameter pressurized glass sphere. Case He/Air, 
Pp, = 826 pei, p, = 14.7 pei, T, = T, = 300°K 
(a) symmetrical view; (b) view for maximum shock path 


late an altitude explosion where the explosive energy 
density has been reduced in the same ratio as the am- 
bient atmosphere. This is simple to do in the present 
case but is not so easily done with TNT, for example. 
The runs were in good agreement in every detail with 
those taken at ambient conditions (at sea level). How- 
ever, this should not be misconstrued as a verification 
of scaling at high altitude, where sizeable excitation of 
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Figure 21 
Comparison of theory and experiment of the wave system 
following the explosion from a pressurized glass sphere. 
Case He/Air, p, = 264 psi, p, = 14.7 psi, T, = T, = 300°K. 
(Note that only the boundaries of the experimental 
results have been given for convenience) 
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internal modes can take place which are both density 
and temperature dependent. Under such conditions it 
is not possible to use the blast data obtained at sea 
level in order to extend it to high altitude by means 
of the energy parameter « in a simple manner. It 
should be noted that TNT cannot be scaled identically 
at altitude, because by analogy with Eq. (10) the 
diaphragm pressure ratio Pa is changed as a result 
of the reduced ambient pressure p:. To reproduce 
identical scaling, a solid explosive of lower energy 
density would have to be used to match the new 
pressure p. (see Reference 18). 





Figure 22 
View of steel sphere, schlieren system, 
and ancillary equipment 


Piezo-pressure records taken at given radii as a 
function of time for these types of explosions are dis- 
cussed in some detail in Reference 25. As these gauges 
effectively measured the unsteady pitot pressure, they 
were useful for determining the instantaneous over-all 
Rankine-Hugoniot pressure jump across the incident 
and reflected shock waves. Miniature blast type gauges 
for measuring flow-pressure profiles have as yet not 
been developed. Two representative pressure records 
for Air/Air and He/Air explosions are shown in 
Figure 23. Good agreement was obtained between the 
pressure records and those of the (1, t)-plane schlieren 
photographs regarding shock arrival and shock 
strength. . 


It is worth noting that the head of the spherical 
rarefaction wave can also be seen when a gas with 
a high index of refraction, such as sulfurhexafluoride, 
is used in a large glass sphere which also helps to de- 
crease the opacity, as shown in Figure 24. The im- 
ploding head moves at a constant sound speed of 
about 450 ft/sec until it reflects from the origin and 
then the speed decreases through the rarefaction wave. 
The second shock is also seen to implode and reflect 
from the origin about 1.3 millisec after the sphere 
breaks, and illustrates that the wave interaction effects 
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persist for a relatively long time in a low sound 
speed gas. 

The effect of increasing the charge size or sphere 
diameter is well illustrated in Figure 25a when it is 
compared with Figure 18. For example, the cor- 
responding distance between S, and S, is reduced by 
half in Figure 25a. Cranz" discovered that times and 
distances in spherical blasts of a given explosive, for 
example, scale directly with the radius of the charge. 
Consequently, if the charge radius is doubled or the 
mass increases eightfold, a given peak overpressure, 
say, will occur at twice the distance and time than 
those for the original charge. This property is already 
inherent in the present numerical solution becatse of 
the definition of + and A, which were made non- 
dimensional by the total explosion energy. 
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(b) He/Air, p, = 326 psi, p, = 14.7 psi, T, = T, = 300°K 


Figure 23 


Oscilloscope records of the response of an SLM piezo- 
pressure gauge to the blast from a 2 inch diameter 
glass sphere at a radius of 5.7 inches 


The effect of a high diaphragm pressure ratio Pu 
is illustrated in Figure 25b, where the shock is propa- 
gating into the 10 mm Hg atmosphere at 2600 ft/sec 
with almost negligible attenuation. Consequently, very 
strong spherical shocks can be produced by using a 
hot, high pressure driver and a low ambient atmo- 
sphere. This result might have been expected since 
at the instant of rupture, conditions at a spherical or 
cylindrical diaphragm are as in an ordinary shock 
tube. That is, the planar solution applies at the instant 
of rupture but immediately after the cylindrical or 
spherical geometric decay takes effect. However, with 
combustion runs significant improvements in the pro- 
duction of strong shock waves should be possible if 
they are coupled with a low pressure atmosphere. 
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(a) (b) 


Figure 24 


Blast from a pressurized glass sphere containing 
sulfurhexafluoride 
(a) schlieren record of (r, t)-plane, 
(b) explanatory sketch: S, = main shock, 

S, = second shock, R = head of rarefaction wave, 
= 5 inch diameter intact glass sphere, 
F = fragments of shattered sphere. 

Case SF,/Air, p, = 70 psi, p, = 14.7 psi, T, = T, = 300°K 





(b) 


Figure 25 


Scaling properties of blast waves 
(a) mass scaling effects of a blast from a 1.0 inch 
diameter pressurized glass sphere. Case Air/Air, 
P, = 326 psi, p, = 14.7 psi, T, = T, = 297°K, 
(b) scaling properties at low pressure atmospheres for a 
2 inch diameter glass sphere, p, = 326 psi, 
P, = 0.013 psi, T, = T, = 295°K, 
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Figure 26 
Combustion driven explosions 
(a) spherical combustion in a 4 inch diameter < 4 inch 
long glass cylinder, p, = 40 psi (0.2 H, + 0.1 O, + 0.7 He) 
Pp, = 14.7 psi (Air), T, = 300°K 
(b) 2 inch diameter glass sphere, p, = 200 psi 
(0.2 H, + 0.1 O, + 0.7 He), p, = 14.5 psi, T, = 303°K 


Experiments using combustion drivers for 2 inch 
diameter to 5 inch diameter spheres have recently 
been performed at UTIA, using a stoichiometric mix- 
ture of oxygen and hydrogen diluted with 70% 
helium (partial pressures)”. The combustible gas was 
initially premixed in two stages, the hydrogen and 
helium from one tank was admitted first and oxygen 
from another was used to complete the charge. The 
three gases are now premixed (strict safety precau- 
tions must be taken) in a single high pressure cylinder 
to improve combustion through the homogeneity of 
the mixture. The mixtures are ignited at the centre of 
the sphere from a small spark obtained from a pulse 
transformer, using a 300 volt primary and 6000 volt 
secondary. The combustion process is illustrated in 
Figure 26a, where, for convenience, it is taking place 
in a 4 inch diameter cylinder, 4 inches long, so that 
it may be observed through the optically transparent 
end plates. The moment the spark occurred is shown 
by the horizontal light time. This timing marker was 
put on using another spark source. It is seen that the 
flame builds up in velocity during the first 200 usec 
and then maintains an almost constant speed of 50 ft/ 
sec up to 1700 usec when it starts to decelerate”, as 
a result of reflected compression waves generated by 
the piston-like flame front, and at 2800 usec the 
cylinder bursts. From these and other results it can 
be concluded that only a slow moving deflagration 
front occurs since the spark is very weak. The full 
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pressure and temperature (theoretically nearly an 
order of magnitude greater than at ambient condi- 
tions)” is not realized in any of the present spheres 
(or cylinders) as they will normally not withstand 
more than about 600 psi. As a result, the spheres burst 
at each particular maximum strength and the com- 
bustion process is quickly extinguished by the incident 
and reflected rarefaction waves. Consequently, al- 
though a noticeable improvement in shock velocity is 
obtained the full benefit of combustion is not realized. 
Figure 26b illustrates this point. It is seen that the 
second shock wave reflects from the origin very early 
(less than 100 uw sec after rupture) owing to the pre- 
dominantly hot helium in the driver gas mixture. An 
acceleration of the main shock wave occurs after 
about 50 uw sec and its maximum speed is only 2000 ft/ 
sec. Therefore, the most significant improvement in 
shock strength would be obtained if the combustion 
drives were coupled with stronger spheres bursting in 
a low pressure atmosphere. Further details will be 
found in Reference 32. 

Some earlier attempts to generate spherical im- 
plosions may be found in Reference 25. In the present 
experiments, the 3 ft diameter steel sphere was pres- 
surized up to 125 psi, which is the limiting pressure 
of the existing windows in the sphere. The glass 
sphere, up to 5 inch diameter, was either evacuated 
and broken with the solenoid-operated breaker or the 
sphere failed automatically at its weakest section, 
owing to the crushing effect of the ambient atmo- 
sphere. Two representative results are shown in 
Figures 27 and 28. 


Figure 27 is a high speed schlieren record taken at 
50 u sec intervals of an implosion from a 5 inch dia- 
meter glass sphere. It is seen that, when the mallet hit 
the sphere, the fragments that formed the lower cap 
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Figure 27 


High speed schlieren record of an imploding 5 inch 
diameter sphere. Case Air/Air, p, = 15 psi, p, = 100 psi, 
T, = T, = 300°K. Times from 1 to 8 in y sec, 900, 1100, 
1300, 1600, 1900, 2200, 2500, 2800: mechanical break 
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: Figure 28 
High speed schlieren record of an imploding 5 inch 
diameter glass sphere. Case Air/Air, p, = 75 psi, 
Pp, = 125 psi, T, = T, = 300°K. Times from 1 to 6 in 
u sec, 600, 800, 1000, 1200, 1800: auto break 


were carried inside the sphere by the inrush of high- 
pressure air. Unlike the explosion cases, the glass 
sphere did not fragment instantly, illustrating that a 
compressive load tends to keep it intact. The very 
interesting fringe pattern that covers the field consists 
of Mach waves that have been generated by the 
spherical rarefaction wave moving into the gas at rest. 
(The knife edge in the schlieren system was vertical 
and to the right looking from the camera. Therefore, 
light areas on the left of any frame are indicative of 
expansions and dark areas of compressions and vice 
versa on the right hand side.) The low pressure region 
around the periphery of the broken sphere caused by 
the inflow is clearly seen in frames 1 and 2 in the 
region of the lower cap. The shower of glass frag- 
ments soon impinged on the remainder of the sphere 
and broke it up, and then emerged through the upper 
cap of the sphere. The resulting flow was completely 
non-uniform. 

Figure 28 illustrates a similar implosion run, but 
this time the sphere was broken by the crushing load 
of the external pressure. The painted conducting ring, 
which is clearly seen, supplied the trigger pulse for 
the camera. The right hand side of the sphere failed 
first as shown in frame 1. This time it is evident that 
the entire sphere has cracked since there is an efflux 
of gas from the sphere cracks on the left hand side. It 
must be concluded that compression or shock waves 
were formed by the piston action of the inrushing 
air, such that the pressure and temperature inside the 
broken: glass sphere locally exceeded the ambient 
pressure. The density of the gas however is still lower 
than ambient (white). The glass fragments formed a 
“shaped charge” that completely penetrated and obli- 
terated the left side of the sphere, as seen in the 
remaining frames. 

A similar type of implosion is shown in Figure 29, 
which is a schlieren photograph of the (r, t)-plane. 
A 5 inch diameter sphere containing air at 300 mm Hg 


- was used under an external pressure of 65 psi of 30% 


SF. and 70% air, in order to get a high index of 
refraction for the schlieren flow visualization. Rare- 
faction wave pulses are seen for about 800 u sec to 
the left of the record (black lines) and when the glass 
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fragments start moving towards the origin, the ex- 
pansion appears stronger. About 400 wu sec later the 
fragments have reached their minimum collective ap- 
parent diameter and are then swept outward by the 
compression wave (broad black band on the right). 
The inner compressed air also is seen to be moving 
out at the left as a lower density (black) flow. Quali- 
tatively therefore, the correct type of implosion re- 
cord is obtained except that the shattering process 
does not take place instantly, and the waves form 
gradually. The interior flow is completely hidden by 
the opaque glass fragments so that it is not possible 
to. tell whether any sharp imploding shock wave 
exists. 


In view of the high speed schlieren records of the 
(x, y)-plane and those of the (7, t)-plane it must be 
concluded that this method of generating spherical im- 
plosions will not prove successful until a type of glass 
“diaphragm” is found that will shatter uniformly and 
very rapidly, and even then interior observations may 
prove impossible. The present results substantiate 
those obtained previously in Reference 25. 


In summary, one can state that the use of pres- 
surized glass spheres for generating spherical explo- 
sions has proved quite successful. The agreement with 
theory is satisfactory, considering the limitations that 
are imposed when a physical diaphragm, which has 
mass and takes a finite time to break up, is introduced. 
The motion of the second shock wave in the interior 
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(a) (b) 


Figure 29 


Spherical implosion 
(a) schlieren record of the (r,t)-plane showing an 
imploding 5 inch diameter glass sphere : 
(b) explanatory sketch: G = glass sphere still intact, 
F = fragment of shattered sphere, p, = 65 psi 
(0.3 SF, + 0.7 Air), p, = 5.8 psi (Air), T, = T, = 298°K 
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Figure 30 


Types of glass spheres and cylinders used for explosion 
and implosion experiments 


(a) 125 mm diameter < 1 mm wall, 60 gm, 
(b) 50 mm diameter < 1 mm wall, 20 gm, in holder, 


(c) 50 mm diameter < 65 mm long X 65 mm long 
< 1% mm wall, 40 gm, open end, 


(d) 50 mm diameter <X 50 mm long < 1% mm yall, 
50 gm, closed end (3.2 mm wall), 


(e) 100 mm diameter < 50 mm long X< 3 mm yall, 
200 gm, closed end (3.2 mm wall), 


(f) 100 mm diameter < 100 mm long X< 3 mm wall, 
300 gm, closed end (3.2 mm wall) 


and exterior gas has been established. There is little 
doubt from the schlieren records that the second 
shock wave does not reach infinite strength as it 
implodes and reflects at the origin as expected for a 
viscous, heat-conducting gas. Attempts were made to 
photograph the luminosity that might have been pro- 
duced at the origin in argon and air but without 
success. This substantiates that even if strong implod- 
ing shocks were generated, the amount of gas affected 
by excitation at the origin is negligibly small. 


Cylindrical explosions and implosions 

During the past few months glass cylinders have 
been used to generate explosions and implosions™ “. 
Open and sealed end cylinders have been used for this 
purpose (Figure 30). The sealed end cylinders were 
made by the Corning Glass Co., with glass end plates 
fused to the cylindrical body with a minimum of 
optical distortion. Both types of cylinders were held 
quite rigidly in a frame between two glass plates, ap- 
proximately 12 inch X 12 inch X % inch, in order to 
assure that a cylindrical flow geometry was main- 
tained. 


The cylinders were either pressurized and ruptured, 
using a solenoid-operated mallet, or they were rup- 
tured by overpressure, when filled with a gas or com- 
bustible mixture of oxygen-hydrogen and helium (as 
for spheres), and an exploding wire was used for 
ignition in case of open end cylinders. Rubber gaskets 
and O-rings held the cylinders against the glass plates, 
which were drilled and fitted with two % inch dia- 
meter, copper electrodes. A 0.0025 inch diameter 
copper wire was stretched across the electrodes in the 
centre of the cylinder, and it, was exploded by dis- 
charging an energy of 50 joules stored in a 25 uf 
capacitor at 2000 volts. The wire ignited the gas mix- 
ture and at a later time the cylinder ruptured and 
generated an explosion. 


A blast from a 2 inch diameter closed end cylinder 
appears in Figure 31. It is seen that the wave system 
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is very similar to the spherical explosion. The head 
of the rarefaction wave H is easily visible now 
through the flat end plates. The first shock S, and 
the imploding and exploding motion of the second 
shock S., as well as the path of the contact front C, 
can readily be distinguished. 


Figure 32 shows a much stronger explosion from 
an open end cylinder. The rarefaction wave and the 
main shock wave are very clear in this record. The 
second shock wave is obscured by the internal flow 
and reflected shock waves. Sound pulses are seen in 
the interior of the pressurized cylinder (50 mu sec) 
prior to its complete collapse. The pulses which arise 
from the initial cracks in the cylinder are not seen in 
the exterior atmosphere, owing to its low density 
perhaps or the asymmetry of the break. 


Figure 33 shows a cylindrical blast, using a stoi- 
chiometric mixture of oxygen and hydrogen diluted 
with 70% helium at a total pressure of 88 psi, which 
was initiated by an exploding wire. The moment the 
wire exploded is indicated by the luminous streak 
across the schlieren record. The shock wave from the 
exploding wire and its attached rarefaction wave move 
back and forth within the cylinder and distort the 
motion of the flame front. The burning mixture ap- 
pears as a white region, and is finally extinguished by 
the expansion waves about 1100 mu sec after initiation. 
The black line through the centre is the position of 
the electrode. It can be seen that the shock wave was 
accelerating and finally was moving at 1700 ft/sec 
at a radius of about 5 inches from the centre of the 
explosion. Since every cylinder has its own bursting 
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(a) (b) 


Figure 31 


Blast from a 2 inch diameter pressurized glass cylinder 
(a) schlieren record of the (r, t)-plane, 
(b) explanatory sketch: G = sealed end glass cylinder 
2 inch diameter < 2% inch long: mechanical break 
F = glass fragments, S, = main shock, S, = second shock, 
H = head of rarefaction wave, C = contact surface. Case 
Air/Air, p, = 115 psi, p, = 14.5 psi, T, = T, = 297°K 
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(a) (b) 


Figure 32 
Cylindrical explosion from a pressurized open end glass 
cylinder 2 inch diameter X 2% inch long 
(a) schlieren record of the (r, t)-plane, 


(b) explanatory sketch: G = open end glass cylinder, 
mechanical break, F = glass fragments, 
H = head of rarefaction wave, case Air/Air, 
P, = 350 psi, p, = 14.7 psi, T, = T, = 300°K 


pressure, which is below the calculated constant 
volume pressure conditions, combustion will continue 
for a short time after rupture and will help to ac- 
celerate the shock front, but because the breaking 
strength of a cylinder is relatively low the full benefit 
of the high pressure and temperature due to com- 
bustion is not attained. The motion of the second 
shock is rather obscured in the combustion runs and 
it is difficult to distinguish. 


From the above it can be concluded that it is quite 
feasible to generate cylindrical explosions using glass 
cylinders by any of the three methods noted. The 
combustion runs are more desirable when strong shock 
waves are required and for their reliability in obtain- 
ing a successful break. However, their initial condi- 
tions are not as well known as in explosions initiated 
with a mechanical breaker. Comparing the cylindrical 
and spherical explosions it can be stated that the 
spherical bursts are much more reliable and sym- 
metrical. They do not suffer from gasket leaks and are 
much simpler to set up. In addition, the cylindrical 
explosions are not free from side-wall boundary layers. 


For the analysis of the cylindrical explosions, 


' Eqs. (12) to (17) may again be solved, employing the 


method developed by Brode. This has not yet been 
done, since it was felt necessary to establish first a 
successful method of generating such explosions prior 
to undertaking extensive calculations. 
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It was noted previously that a real difficulty with 
spherical implosions from ruptured glass spheres is 
that the interior flow is obscured by the opacity of 
the fragmented glass. To overcome this difficulty some 
cylindrical implosions were initiated using sealed end 
and open end cylinders“. The sealed end cylinders 
were not found to be useful for this work, as the end 


; plates shattered and the net effect was an obstructed 


interior view nearly similar to the spherical implosions 
seen above. Large open end cylinders appear to offer 
some promise and their use is presently being investi- 
gated. Analyses also in this case will not be attempted 
until a feasible experimental method has been estab- 
lished. It should be noted that Payne” has numerically 
solved a number of cases for the cylindrical implosion. 
No second shock wave was found for the reasons 
given previously. Payne used a form for the artificial 
viscosity that did not prove as useful as the one em- 
ployed by Brode (Eq. (17)), since it appeared to 
affect the entire flow field to some extent. It would be 
quite worthwhile to repeat this work using an 
equation of the form given by Eq. (17) for the arti- 
ficial viscosity. Such calculations would be initiated 
at UTIA if a successful method of generating cylin- 
drical implosions were developed. 


Planar analogues of explosions and implosions 
If a very short chamber relative to the channel 
is used to establish a shock tube flow, then the rare- 
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(a) (b) 
Figure 33 
Combustion driven cylindrical explosion 
(a) schlieren record of the (r, t)-plane, 


(b) explanatory sketch: G = open end cylinder 2 inches 
< 2% inches; auto break, F = glass fragments, 
E = exploding wire, W = wire generated shock, 
B = burning flame, p,= 88 psi (0.2 H, + 0.1 O, 
+ 0.7 He), p, = 14.7 psi, T, = 299°K 
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faction wave will quickly reflect from the closed end, 
interact with the contact region, overtake the shock 
wave and establish a process of shock decay. A peaked 
type shock wave profile is then formed which re- 
sembles a blast wave. If the initial shock is very strong 
then the shock motion may perhaps have a type of 
strong planar blast wave decay similar to Eq. (8), 
after it had engulfed a certain mass of air correspond- 
ing to the spherical cases shown in Figures 7 and 11. 
In the planar case, a graphical or numerical solution 
using the method of characteristics would be quite 
feasible, especially since many of the interaction 
problems can be solved in closed form”. 


In a similar manner, when the channel of a shock 
tube is made very small compared with the chamber, 
then a flow analogous to an implosion can be gen- 
erated. The shock wave “implodes” on the origin with 
uniform strength (end of channel), reflects from the 
closed end with an increased uniform strength, inter- 
acts with the contact surface and then with the rare- 
faction wave where it is decayed. Of course, it is only 
the attenuation of the shock wave that bears an 
analogy with the cylindrical and spherical blast waves 
in the explosion or implosion cases. Before that time 
the waves and states are uniform in the planar case 
and non-uniform in the cylindrical and spherical cases. 
It would be of some interest to investigate these two 
problems analytically and experimentally in a shock 
tube such as the Wave Interaction Tube at UTIA, 
which is well suited for this purpose. 


BLAST WAVE SIMULATORS 


The effects of blast waves can be simulated by 
different methods and to different degrees. Each 
method has its own inherent advantages and draw- 
backs. The following methods have been used: 


(1) full scale tests’, 

(2) chemical explosions in open ranges” ”, 

(3) chemical explosions in closed ranges”, 

(4) shock tubes, 

(5) conical sector tubes”, and 

(6) explosively driven (Primacord) pressure vessels”. 


A full discussion of these facilities is beyond the scope 
of the present paper. However, it is worth noting that 
complete simulation can only be achieved in a full- 
scale test but with the inherent disadvantages of 
hazards, expense, and the lack of laboratory type of 
control and repeatability. The other facilities may 
simulate various properties of explosions and some of 
these better than others. For example, the blast from 
pressurized glass spheres or cylinders is most useful 
for the investigation of the shock wave and the flow 
properties, since the overpressures and their durations 
are comparatively small for significant size model tests. 
Shock tubes are very useful for the determination of 
aerodynamic loads on rigid models in a blast“ and for 
proving structural components. Typical schlieren 
photographs of the flow over slender“ and gable 
models with relatively short and long shock-transit 
times are shown in Figures 34a and 34b. In both cases 
the shock diffraction over the models gives rise to 
Mach-type shock reflection, which includes a vortex 
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(a) diffraction of a plane shock wave over a slender 
model. Model thickness 4% inch, p, = 185 mm Hg, 
P,, = 5.20, M, = 1.04, 2 inch X 7 inch shock tube, 
primary shock moves from left to right 





(b) diffraction of a plane shock wave over a gable model. 
Primary shock wave moves from right to left and has a 
subsonic flow behind it 


Figure 34 


and a slipstream (contact surface). In Figure 34a the 
flow behind the shock wave is just supersonic and in 
Figure 34b it is subsonic and contains a series of Mach 
waves generated in the flow. For a slender solid model 
(Figure 34a), with a short shock-transit time or rapid 
pressure equalization, aerodynamic drag may be pre- 
dominant. However, for a sizeable panelled structure 
(Figure 34b) with a long transit time, then peak pres- 
sures and the pressure impulse may cause the most 
severe damage. 

For dynamic testing of sizeable structural com- 
ponents, a shock tube facility would become enor- 
mously large and expensive. For such tests high over- 
pressures and long flow durations (of the order of the 
response time of the structure) are required, and 
explosively driven simulators, such as specially de- 
signed high pressure vessels” or sector tubes”, appear 
to be very useful. It is worth noting that since all of 
the blast energy in a sector tube can be concentrated 
over a small solid angle rather than over 4 steradians, 
very large amplification factors are theoretically pos- 
sible for very small cone angles [(4/d6)’*, or 5.25 X 10* 
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for a 1° cone angle]. Consequently, relatively small 
charges can be used to produce effectively large ex- 
plosions. In practice, such high magnification factors 
are not obtained owing to losses that arise from elastic 
and viscous effects. However, the experimentally ob- 
tained gain is still vey impressive™. 

A blast from relatively small hemispherical TNT 
charges exploded in a steel hemisphere of about 50 
to 100 ft in diameter, for example, would only have a 
magnification factor of two for the same charge 
weight. However, it would be more versatile than a 
narrow sector tube with regard to the number of 
models that might be tested per run at different radii, 
and it would have the advantage of simulation of 
terrain and water surfaces as well as the possibility of 
testing at low pressure (high altitude) testing. With 
a 7 lb hemisphere of TNT (1060 cal/gm) it would 
be possible to obtain an overpressure of 25 psi at 
11 ft from the origin, about 3.4 millisec after the 
initiation of the blast. The positive phase would last 
about 2.2 millisec and the negative phase about 5.6 
millisec. At a 50 ft radius, the overpressure behind the 
blast wave is down to 1.5 psi and should not present 
any structural difficulty. On the other hand, a pressure 
of 25 psi is reached at about 1560 ft in a 20 kiloton 
TNT blast, giving pressure profile durations of about 
140 times as great as for the 7 lb hemisphere. Conse- 
quently, for dynamic structural tests, flow times of 
several hundred milliseconds would be required. 


SOME RESEARCH APPLICATIONS OF EXPLOSIONS 
GENERATED FROM PRESSURIZED GLASS SPHERES 
AND CYLINDERS 

It was noted previously that successful explosions 
can be generated from pressurized glass spheres and 
cylinders. Blasts produced in this manner are relatively 
safe in comparison with chemical sources, and the 
initial conditions are well known. The spherical ex- 
plosions have the added advantage of the absence of 
any viscous boundary layer effects, and are even 
simpler to generate than cylindrical blasts. Conse- 
quently, spherical explosions produced in this manner 
can be applied to a variety of interesting research 
problems, such as underwater explosions, blasts in a 
supersonic flow”, shock wave collisions“, shock wave 
diffractions over various models, shock wave refrac- 
tion at gaseous and liquid interfaces and shock waves 
at extremely low densities without boundary layer 
interference. Properties of cylindrical and spherical 
deflagration and detonation waves may well be investi- 
gated using the same techniques. 


In the following subsections a few of the problems 
noted previously will be considered. As soon as the 
flow loses complete spherical symmetry, the analysis 
becomes very complex. It is necessary to introduce 
an additional space variable and more difficult bound- 
ary conditions. Consequently, problems on spherical 
shock wave collision, diffraction and refraction have 
as yet not been solved numerically, using the methods 
developed by Brode, as they tax the capabilities of 
present day high speed computers. It is hoped that 
such solutions will become possible in the near future. 
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In the meantime experimental feasibility studies have 
been made of some of the problems outlined pre- 
viously. 


Underwater explosions 


It is quite feasible to study underwater blasts gen- 
erated from pressurized glass spheres, using compressed 
gases or combustible mixtures as drivers. Some pre- 
liminary runs have just been conducted in a commer- 
cial aquarium (8 inches X 16 inches X 10 inches deep) 
filled to a depth of about 8% inches of water. Distilled 
water is most satisfactory as it reduces light absorption 
and scattering. Pressurized glass spheres 2 inch dia- 
meter containing air over 300 psi were used to pro- 
duce the blasts, which were recorded with the high 
speed schlieren camera. One run containing 200 psi of 
a combustible mixture (0.2 H: + 0.1 O. + 0.7 He) 
was also photographed and a schlieren record of the 
(7, t)-plane was obtained. All the runs were self-breaks 
due to overpressurization. A trigger pulse to the spark 
light sources was provided by a painted circuit 
circling the glass sphere which was severed when the 
sphere burst. 


High speed schlieren photographs of an under- 
water explosion generated from a 2 inch diameter 
glass sphere at 310 psi are shown in Figure 35a. In 
this experiment the schlieren system was not perfectly 
adjusted for the water, but rather for the air as the 
working medium. The record also suffers from a slight 
accidental double exposure which accounts for the 
flake-like background in some frames. Nevertheless, 
it does illustrate the blast phenomena quite well. In 
frame 1, the air-water contact surface appears as a 
black band and above it is the wider band of the 
aquarium metal framework. The sphere and its metal 
pipe holder are also seen. Tap water was used in the 
tank for this run. Frame 1 was taken 15 wu sec after 
the sphere ruptured. The hydrodynamic blast wave 
appears well formed but somewhat displaced below 
the ruptured sphere, which still appears intact. Unlike 
an air burst the glass fragments are constrained by 
the water from flying off. Taking half of the differ- 
ence between the shock and glass sphere diameters, for 
a distance base line, gives the shock wave an average 
velocity of 1.84 inches in 15 mu sec or 12,300 ft/sec. 
Using a sound speed in water of 4800 ft/sec results 
in an apparent shock Mach number Ms ~ 2.12. Em- 
ploying the calculated shock wave properties in water 
given in Reference 42, a corresponding pressure of 
nearly 35,000 atm would exist behind a blast wave at 
this shock velocity. Consequently, it must be con- 
cluded that the electronic timing, any optical re- 
fraction effects that might produce an enlarged shock 
wave image, and the encasement effect of the glass 
sphere will have to be carefully investigated before 
the actual velocities of the shock wave are accurately 
established. 


Since vertical knife edges (to the right, lgoking 
from the camera) were used, the upper and lower 
portion of the shock do not appear so well defined. 
The left portion of the shock also suffers from the 
removal of light from an already weakly illuminated 


129 








3 


(a) High speed schlieren photograph of an underwater 
explosion generated from a pressurized glass sphere. 
Aquarium 8 inches X< 16 inches < 10 inches, water depth 
8% inches, 2 inch diameter glass sphere pressurized at 
310 psi, surrounding conditions atmospheric, framing rate 
from 1 to 8 in ,, sec, 15, 60, 100, 200, 350, 700, 1300, 2000 





(b) High speed schlieren records of underwater explosions 
generated from pressurized glass spheres. Aquariums: 
8 inches < 16 inches < 10 inches, water depth 8% inches. 
Row No. 1: 2.1 inch diameter sphere, bursting pressure 
300 psi, framing rate from 1 to 4 in y sec, 11, 30%, 50, 
70. Row No. 2: 2.1 inch diameter sphere, bursting pressure 
365 psi, framing rate as above. Row No. 3: 2 inch diameter 
sphere, bursting pressure 330 psi, framing rate 201% and 
50 uw sec. Ambient, temperature 300°K at 1 atm 


Figure 35 


background. Many overtaking compression waves are 
seen behind the shock wave in the white area (com- 
pression) to the right of the frame. Unlike an air 
burst, the characteristic change from a white to a 
black area, as the gradients from the rarefaction zone 
of the peaked profile take effect behind the shock 
wave, is not seen here. The reflected rarefaction pulses, 
generated at the air-water contact surface as the main 
blast wave and overtaking compression or Mach waves 
refract at the contact surface, are seen in frame 2. 
These waves do travel with the sound speed in water 
at about 4800 fps, as will be noted later. 


The compressed air bubble is seen to emerge in 
frame 4 and does not appear to have changed in 
frame 5. In frame 6 it has grown to a 3.0 inch dia- 
meter, in frame 7 to a 3.9 inch diameter and in 
frame 8 to a 4.6 inch diameter, giving an average 
velocity of about 60 ft/sec and 40 ft/sec, respectively. 
The air bubble appears effectively opaque, as a result 
of the spherical lens action of the surrounding water, 
which deflects the light rays very much like the glass 
sphere itself. In frame 7 the expanding bubble is now 
raising the water surface and in frame 8 this effect is 
still more apparent. Even in this last frame, after a 
flow duration of 2 millisec, the glass sphere still ap- 
pears intact and its diameter is almost unchanged. 


Figure 35b contains frames from three additional 
underwater explosions taken in similar aquariums (the 
glass Side plates invariably broke during a run) all 
filled with distilled water. The schlieren system was 
adjusted for the water as the working medium. The 
improved quality of the schlieren photographs is quite 
noticeable. The advance of the blast wave and its re- 
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fraction at the contact surface are well illustrated in 
row No. 1. The apparent average shock speed from 
frame 1 is 2 inches in 11 uw sec or 15,200 ft/sec; from 
frame 2, 3.1 inches in 30.4 sec or 8500 ft/sec. The 
difference between frames 2 and 3 gives 1.16 inches in 
19.4 4 sec or a rapid decay to nearly 5000 ft/sec. The 
advance of the reflected rarefaction wave head in 
frames 2, 3 and 4 gives it a speed very close to the 
sound velocity of 4800 ft/sec. 


In row No. 2, the blast wave appears to have 
travelled 2.5 inches in 11 mw sec or an apparent speed 
of 19,200 ft/sec. It is worth noting that the bursting 
pressure was also higher here. The reflected rare- 
faction wave head again has a velocity of nearly 
4800 ft/sec, but somewhat slower. (A measurement 
of distance does not have the same accuracy as that of 
time. However, the velocities quoted are not meant 
to be definitive as the experiments are only of a pre- 
liminary nature.) In row No. 3, an average shock 
velocity of 7800 ft/sec is obtained for the first 
20% m sec and is consistent with the previous results, 
as the shock decay is rather’rapid. From frame 2, the 
average speed over the first 50 mu sec is 6100 ft/sec. 
The difference between frames 1 and 2 again gives an 
average velocity of 5000 ft/sec, which is approaching 
the acoustic value and is in agreement with the pre- 
vious result. 


The above preliminary data have shown that this 
method of generating underwater blasts offers some 
very worthwhile possibilities, and a systematic re- 
search program is presently underway. The apparent 
values quoted for the shock speeds were given for 
illustration purposes only and should not be miscon- 
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strued as being actual speeds. Further details will be 
found in future UTIA publications as they become 
available. 

It is perhaps worthwhile to include the (7, t)-plane 
schlieren record of an underwater explosion, shown 
in Figure 36. A combustible mixture (0.2 H: + 0.1 Os 
+ 0.7 He) at 200 psi was used as the driver gas. The 
instant of ignition is shown by the luminous horizontal 
line and the first shock wave appears about 1400 u sec 
later. Poor quality schlieren (due to the flexing of the 
thin glass plates perhaps) make the passage of the 
shock wave almost invisible except at the very edges 
of the film strip. In this case, the hot air bubble in 
contact with the water gives rise to some peculiar 
spreading and motions of the contact region. "The 
record was included for its qualitative value only but 
it does illustrate that good quantitative records should 
be possible. Schlieren photographs of the (7, t)-plane 
would be much more definitive in accurately com- 
puting the shock and air bubble velocities from these 
continuous-time records. 


Head-on collision of spherical shock waves 

An investigation of this type of interaction, using 
blasts from 2 inch diameter. glass spheres, has proven 
quite successful. Since numerical solutions existed for 
Air/Air and He/Air explosions for spheres at 22 and 
18 atm, respectively, these were chosen for the initial 
studies. Consequently, by the time the spherical shock 
waves collided they were relatively weak. Much 
stronger collisions can be generated by using low am- 
bient pressures (Figure 25b) coupled with combustible 
driver gases. The boundary condition at the point of 
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(a) (b) 


Figure 36 

Underwater explosion generated from a pressurized glass 

sphere (a) schlieren record of the (r, t)-plane, 

(b) explanatory sketch: G = glass sphere, 2 inch diameter, 
filled with 200 psi (p,) of a combustible mixture 
(0.2 H, + 0.1 O, + 0.7 He); ambient conditions, 
aquarium 8 inches < 16 inches < 10 inches, water 
depth 8 inches, surrounding atmosphere 14.7 psi at 
300°K, S = main shock wave, C = constant region 
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Head-on collision of planar shock waves in the (u, p)- and 
(x, t)-planes. Conditions in a shock tube before collision 
(t = t,) and after collision (t = t,) are indicated 


intersection of the spherical shock waves can be deter- 
mined from planar wave theory”. This type of inter- 
action is illustrated in Figure 37 in the (x,t)- and 
(u, p)-planes. State (1) is a rest state. States (2) and 
(5) are the states behind the backward and forward 
facing approaching shock waves. It is seen that after 
the collision two shock waves that are separated by a 
contact surface recede from the point of interaction. 
The new states (3) and (4) have a higher pressure 
than the previous states and their particle velocity 
will follow the stronger shock wave. Once the initial 
shock strengths P., and Ps: and state (1) are known, 
the final wave strengths and states (3) and (4) can 
be found“. 


After the instant of interaction, new boundary 
conditions must be satisfied that involve both regular 
and Mach-type reflections. Consequently, the flow 
becomes asymmetrical and very complex. A solution 
would now contain three independent variables (two 
spatial and one temporal) and the reflection boundary 
conditions noted above. Such a solution based on the 
method developed by Brode has not been obtained 
to date. 


However, the collision problem has been investi- 
gated using optical techniques, and a comparison be- 
tween the planar theory and the experimental results 
at the point of intersection was possible. A view of 
the two 2 inch diameter glass spheres, with their 
centres 9 inches apart, and of the solenoid-operated, 
spring-loaded breaker mechanism, is shown in 
Figure 38. The same technique was employed in gen- 
erating a double blast as for a single explosion, with 
the exception that the plunger now had two mallet 
heads to break the two spheres at approximately the 
same instant. Although there was some variation in 
sphere geometry and strength of the glass this was 
usually achieved. 


Figure 39a shows a shadowgram of approaching 
shock waves from He/Air and Air/Air explosions. 
Some diffracted waves about the breaker can be seen, 
otherwise the blasts are as in Figure 16. In Figure 39b 
the shock waves have collided and penetrated. Only 
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Figure 38 


Breaker mechanism used for head-on collision of 
spherical shock waves 


a small change in shock wave curvature is noted owing 
to the low strength of the colliding shock waves, and 
for the same reason no contact front can be observed 
in the penetration region between states (3) and (4). 

A schlieren record of the (7, t)-plane of the same 
type of interaction is shown in Figure 40. The left- 
hand sphere contains helium at 500 psi and the right- 
hand sphere air at 250 psi. The approximate location 
of the contact surface is shown in the explanatory 
sketch. The absolute velocities of the shock waves, 
measured in laboratory coordinates, divided by the 
sound speed in state (1), are given in the caption as 
measured absolute shock Mach numbers. These are in 
fair agreement with theory. There are some difficulties 
involved in accurately measuring shock speeds from 
curved paths, and further details are given in Re- 
ference 14. 

The second shock waves, in both blasts, are clearly 
visible. The one from the helium explosion appears 
sooner owing to the higher sound speed. The station- 
ary character of the contact surfaces, after about 
300 usec from the start of the blast, is clearly illus- 
trated and is in keeping with the model noted pre- 
viously for the generation of the second shock wave. 
In this record the glass fragments attain speeds of 
210 and 240 ft/sec for the air and helium blasts, re- 
spectively. The fragments collide about 1500 mu sec 
after the initiation of the explosions. At approximately 
this time the shocks reflected from the steel walls of 
the 3 ft diameter sphere are also’seen, and they implode 
and reflect at the origin. This process is repeated 
manifold. 


Figure 41 shows the details of the interaction 
region in the (r,t)-plane for the case of two 2 inch 
diameter glass spheres under identical conditions, 
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He/Air, 326 psi/14.7 psi, T = 295°K, except that the 
left sphere breaks 160 u sec earlier. As a result a col- 
lision of unequal shock waves takes place as shown. 
Figure 42 shows a typical schlieren record of the 
head-on collision process, taken at nominal time in- 
tervals of 50 u sec, or at a framing rate of 20,000 pic- 
tures per second. In this case, both spheres contained 
helium at an initial pressure of 326 psi. The series of 
photographs start at frame 1, where it is seen that the 
left sphere breaks first and as a result the shocks 
when they collide are of unequal strength. The de- 
velopment of the wave systems appears in frames 1 to 
8, inclusive, and covers an interval from 100 to 
500 « sec. The shocks approach the collision point in 
frame 5 and have collided in frame 6. In frames 7 
and 8 the forward facing shock is interacting with 





(a) left sphere, Air/Air, 400 psi/14.7 psi. Right sphere, 
He/Air, 326 psi/14.7 psi, T = 298°K, 2 inch diameter 
spheres, time delay 297 sec 





(b) left sphere, Air/Air, 400 psi/14.7 psi. Right sphere, 
He/ Air, 326 psi/14.7 psi, T = 295°K, 
time delay = 373 y sec 


Figure 39 


Shadowgrams of the head-on collision of spherical 
shock waves 
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Figure 40 
Schlieren record of the (r,t)-plane showing the head-on 
collision of spherical shock waves. Left sphere, He/Air, 
500 psi/14.7 psi, right sphere, Air/Air, 250 psi/14.7 psi, 
T = 296°K, S,, = 1.04, S,, = 1.24, S,, = 1.07, S, = 0.91 


the contact region from the right sphere. The station- 
ary character of the contact regions beyond frame 4 
is quite noticeable. 


Figure 43 shows a similar result at a slightly later 
time of 252 to 605 u sec. Of particular interest in this 
series is the second shock, which has reflected at the 
origin and is seen to pass through the contact region 
in frame 4, and appears prominently in frame 6. The 
second shock wave from the right-hand sphere is 
more difficult to detect, since the schlieren knife edge 
is on the right and the shock tends to blend into the 
dark background on the photograph. Further ex- 
amples and interesting details may be found in Re- 
ference 14. The above work has indicated that it is 
possible to generate head-on collision interactions. The 
actual experiments were conducted at low shock 
strength because of the availability of numerical solu- 
tions only for the primary blast cases. However, the 
same methods can be used to generate stronger shock 
collisions by simply using low pressure ambient atmo- 
spheres. A complete theoretical-experimental com- 
parison was not possible owing to the lack of numeri- 
cal solution for an asymmetric problem of this type. 
Where it was possible to apply planar wave theory, 
such as to the point of interaction, fair to good agree- 
ment was obtained. 


CONCLUSION 


Analytical solutions exist both for explosions and 
implosions when they are very intense and for ex- 
plosions when they are very weak. For example, 
Taylor’s self-similar solution for an intense point- 
source blast has been substantiated from shock speed 
measurements of an atomic explosion at small times, 
when the assumptions of constant energy and a per- 
fect gas do not appear to seriously affect the solution. 
At later times when the strong shock conditions cease 
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Figure 41 
Detailed view of the collision of two unequal spherical 
shock waves. Both spheres 2 inch diameter, He/Air, 
326 psi/14.7 psi, T = 295°K, time interval 
between breaks 160 y, sec 





Figure 42 
High speed schlieren photographs of the collision of two 
spherical shock waves. Both spheres 1.9 inch diameter, 
He/Air, 326 psi/14.5 psi, T = 296°K, time delays in y see 
for frames 1 to 8, 100, 150, 250, 300, 350, 400, 450, 500 


to apply and considerable energy has been lost 
through radiation, the experimental results are in good 
agreement with the numerical solution for the blast 
generated from a high pressure sphere of very hot 
gas. The analytical extensions to intense cylindrical 
and planar blast waves do not appear to be well repre- 
sented by exploding wires or magnetogasdynamically 
driven shock waves in a T tube. In these cases it may 
be necessary to treat the explosion as a finite source 
with a more complex type of relation for the addition 
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Figure 43 


High speed schlieren photographs of the collision of 
spherical shock waves. Both spheres, 1.9 inch diameter, 
= 295°K, time delays in y, sec 
for frames 1 to 8, in 50 y sec intervals, 250 to 600 


He/ Air, 326 psi/14.5 psi, 


of the blast energy. However, the use of blast wave 
theory for the analysis of the related problems in 
steady, hypersonic, two-dimensional and axisymmetric 
flows have been very successful. 


The theoretically predicted wave systems from 
finite source explosions generated by pressurized glass 
spheres, including the motion of the second shock 
wave, have been satisfactorily verified from schlieren 
records of the (7, t)-plane, considering that limitations 
are imposed by an actual diaphragm of finite strength 
and mass during the early development of the blast. 
Ample scope for the measurement of pressure, tem- 
perature, density and velocity profiles still remains in 
order to fully substantiate the numerical predictions. 
Various driver gases and combustible gas mixtures 
have been successfully used to produce explosions of 
varying initial energy and intensity. As for intensity, 
a low pressure atmosphere provides a simple means of 
producing initially strong blast waves. 


The technique of using pressurized glass spheres 
has been successfully extended to the production of 
cylindrical explosions by using pressurized open or 
sealed end glass cylinders. Compressed gases and com- 
bustible gas mixtures ignited by exploding wires have 
been used as drivers. However, the cylindrical ex- 
plosions are not as consistently symmetrical or as 
easily generated as spherical blasts. The same methods 
when applied to spherical and cylindrical implosions 
have to date not yielded the desired results, owing to 
the severe asymmetrical bursting properties of glass 
spheres and cylinders under compressive loading. 
More sophisticated materials or rupturing methods 
may be required. 


The attenuation of spherical shock waves from 
various finite source explosions appear to have a decay 
rate similar to that of a point-source explosion, after 
the blast wave has engulfed a mass of air equal to 
approximately ten times the initial mass of the ex- 
plosive. Consequently, it might appear that explosions 
from different sources would be completely scaleable. 
However, as Brode points out, the initial history of 
an explosion persists and identical scaling may not be 
possible even if the decay rate is similar. However, for 


134 


the purposes of some aerodynamic 
simulation this may not be a seri- 
ous obstacle, and for structural 
simulation even cruder types of 
explosion generators may be used 
as long as they supply a blast 
pressure profile of sufficient dura- 
tion compared with the response 
time of the structure. 


Use can be made of pressurized 
glass spheres to generate more 
complex types of blast phenomena 
and wave interactions. Some pre- 
liminary results of underwater ex- 
plosions have already substan- 
tiated the worth of this method 
and a detailed study of such blasts 
is now underway. The head-on 
collision of spherical shock waves 
has been successfully produced 
using this method for weak shock interactions, since 
for this case only were numerical solutions of the 
initial explosions available. The experimental work 
could readily be extended to stronger shock waves, but 
a comparison with theory would not be possible be- 
cause numerical solutions of the interacting flows are 
too difficult to handle on present day, high speed com- 
puters. Such a comparison was only possible near the 
point of the shock intersection, where only the ap- 
plication of planar shock wave interaction theory gave 
fair to good agreement with the shock wave paths 
in the (r, t)-plane. 

Similar remarks regarding numerical analyses 
would also apply to studies of spherical and cylindrical 
shock diffraction and refraction in gases or in water, 
or to explosions from pressurized glass spheres and 
cylinders in supersonic flows. Pressurized glass cylin- 
ders and spheres also appear to have useful applications 
in combustion investigations and for spherical shock 
wave transition studies in a rarefied atmosphere. The 
latter would have the advantage of a complete absence 
of interference from the viscous boundary layer, 
which is always present in shock tube flows. 
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DESIGN OF AN INLET DUCT 


FOR THE PROPELLER-TURBINE INSTALLATION 


ON A STOL AIRCRAFTi 


by Dr. D. H. Henshaw,* A.F.C.A.L 


The De Havilland Aircraft of Canada Limited — 





Figure 1 
STOL Caribou 


_—— De Havilland Aircraft of Canada Limited is 
currently designing a propeller-turbine installation 
for its STOL Caribou. The Caribou is a twin-engined 
utility aircraft capable of performing exceptionally 
short takeoffs and landings from rough, unprepared 
fields. It is currently in production and powered by 
the reliable Pratt and Whitney R2000 reciprocating 
engine (Figure 1). 

With STOL aircraft, the engine installation must 
be such as to provide the maximum possible power 
for the takeoff manoeuvre. Special design problems 
are encountered with the Caribou because of the low 
takeoff speeds permitted by the exceptionally effective 
flaps, which operate in the strong flow of air from 
the propeller. The installation must also be suited to 
the high speed cruise condition. 

The turbine engine that is being fitted to the 
Caribou is manufactured by the General Electric 
Company under the designation T-64. This unit has 
a remote propeller gear box driven by a shaft protrud- 
ing forward along the axis of the compressor. The 
general arrangement of the engine may be seen in 
Figure 2. 

With the T-64 type of engine, the airframe manu- 
facturer assumes responsibility for the air supply to 





+Paper read at the Joint I.A.S./C.A.I. Meeting in Montreal on 
the 18th October, 1960. 
*Project Engineer 
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the compressor. On the Caribou, the air is con- 
veniently supplied by means of an S shaped duct with 
entry above the propeller spinner. The flow is then 
over the gear box, downwards past the power shaft 
and thence into the engine intake. Figure 3 illustrates 
the T-64 engine and intake installation as adopted for 
the DHC-4 turbine application. 

The design of the intake is of interest because of 
its influence on engine performance. Engine power 
can be increased, fuel consumption decreased and en- 
gine operation improved by the use of a duct that 
supplies the engine with a uniform flow of air at the 
highest possible pressure. Performance calculations 
have shown that appreciable increases of payload can 
be achieved through the careful refinement of the 
nacelle and duct details. 

The inlet arrangement chosen for the T-64 is in 
the form of 4 scoop above the propeller spinner. The 
inlet is canted and raked to provide the optimum ar- 
rangement for pressure recovery and drag at both low 
and high speeds. A high critical Mach number has 
been assured by an appropriate choice of NACA 
Series 1 Cowl lines’. A boundary layer by-pass has 
been fitted to prevent ingestion of the disturbed air- 
flow from the propeller root sections. 

The difficult factors in the inlet design were the 
wide ranges of angle of attack and inlet velocity ratio. 

High angles of attack tend to give external flow 
separation over the upper intake lip, particularly for 
the “dead” engine in single-engine climbs. Highly 





Figure 2 
Intake duct on General Electric T-64 
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INTAKE DUCT 
Tee ENGINE 





ENGINE WACELLE 


Figure 3 
Caribou installation of T-64 turbine 


cambered cowl lines were adopted to minimize drag 
under these conditions. In addition, a slight downward 
inclination of the intake axis reduced the effective 
angle of attack. The rake of the inlet reduces the 
tendency to external flow separation on the upper lip 
at high angles of attack. 

For the high inlet velocity ratios occurring in 
takeoff, some internal flow separation will occur on 
the inner lips of the intake duct. According to 
Seddon’, the duct immediately downstream of the 
inlet should be so proportioned that flow reattachment 
can take place. This generally requires a constant area 
duct for a length of about one inlet radius immediately 
downstream of the inlet. Inlet losses have also been 
reduced by the use of as large an inlet area and in- 
ternal lip radii as are compatible with the external 
flow requirements in the high speed cruise condition. 

In the formulation of the internal design of the 
duct, calculations were carried out to determine the 
nature of the potential flow solutions. The vectorial 
equations describing the flow in the duct are 


=> 
pu, = — pgrad@ 
s > 
div (pu,) = 0 


To simplify the equations, the flow was assumed 
to be incompressible and so the above equations re- 
duce to the familiar La Place equation y*¢ = 0. 

To avoid the formidable difficulties encountered in 
analytical or digital calculations, it was decided to use 
analogue methods. 

The analogue considered in the present connection 
was one relying on the properties of an isotropic 





Figure 4 
Electrical analogue for duct flow 
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conducting medium. The current flow is described as 

— 

1 =ograde 

= —_— 
div i =0 
—, 

where i is the current vector, e is the voltage and o 
is the conductivity. When o is independent of posi- 
tion, the equations reduce to the La Place equation 
V7e=0 

The velocity potential for the fluid flow and the 
electric potential both satisfy the La Place equation. 
Experimental solutions for the electric potential may, 
therefore, be interpreted as solutions for the velocity 
potential ¢. 

One of the more prominent names in the literature 
of electrical analogues is Malavard. He and his co- 
workers have carried out an extensive development of 
analogue techniques for the solution of many varied 
problems. His early studies were based to some extent 


on the pioneering efforts of E. F. Relf, 1924, and G. I. 
Taylor, 1928. 
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FIGRE GLASS MOULDED TANK 
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Figure 5 
Electric analogy tank for duct 


Studies of three-dimensional flow problems have 
been rather limited, because of difficulties encountered 
in the representation of lifting surfaces. The general 
literature refers almost exclusively to the two-dimen- 
sional flow problems associated with plane or axially 
symmetric flow. In fact, Dr. J. H. T. Wade* has re- 
cently discussed the application of analogue tech- 
niques to the design of an axially symmetric air intake 
for a small gas turbine. Popularity of the two-dimen- 
sional analogues is perhaps limited by the availability 
of alternative numerical methods. For the case of two- 
dimensional plane flow, a rapid approximate analysis 
has been developed by Valentine’. In addition, the 
relaxation methods of Southwell may be applied to 
both plane and axially symmetric flows. 

For three-dimensional flows, the alternatives to the 
analogue techniques are generally unattractive because 
of the effort required for even approximate solutions. 
An invesigation was therefore made to ascertain the 
suitability of analogue techniques for the study of 
flow in a three-dimensional duct. As a result of dis- 
cussions with R. J. Templin of the National Research 
Council, Ottawa, it was decided to proceed with an 
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electrical analogy using a conducting fluid such as tap 
water. 

Following some exploratory measurements with a 
two-dimensional analogy, a three-dimensional model 
of the duct was molded to exact contour using a 
fiberglass material. The electric potential of 110 volts, 
60 cycles, was applied through copper plates at each 
end of the duct. The electric potential distribution 
was measured with a voltmeter as shown in Figure 4. 
Note that the symmetry of the duct has permitted the 
adoption of the half model technique (Figure 5). 


Many of the measurements were concerned with 
the determination of the flow velocities (Figure 6). 
For this purpose, a dual probe was used; this had two 
electrodes at a fixed distance apart. The voltage differ- 
ence was proportional to the component of velocity 
in the direction of the line joining the electrodes. 


Velocity measurements were made at the engine 
inlet. These were useful for the development of a 
design capable of meeting the engine manufacturer’s 
specification for flow uniformity. In this investigation, 
an effort was made to assess the sensitivity of the inlet 
flow velocities to the method of simulation of the 
compressor inlet. It was anticipated that the com- 
pressor would tend to smooth out the distribution of 
flow into the engine. This smoothing effect was simu- 
lated by the introduction of a perforated sheet of 
plexiglass at the compressor face position. The rela- 
tively high voltage drop across the plexiglass resulted 
in a more uniform flow of electric current at the 
compressor face. The measurements showed that com- 
pressor restraint effects were small compared with the 
turning effects. 


Another calculation was performed to determine 
the flow velocities in the neighbourhood of the power 
shaft. This shaft was circular in cross section but was 
to be fitted with a streamline fairing. The beneficial 
effects of a favourable pressure gradient on the fairing 
induced flow losses have already been noted by 
Seddon and Spence’. 


Calculations showed that the original fairing pro- 
posal was situated in an adverse pressure gradient and 
so could be expected to generate relatively large flow 
losses. It was determined that a somewhat longer fair- 
ing would position the trailing edge in a favourable 
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Figure 6 
Electrical analogue flow velocity on plane of symmetry 


gradient. The resultant beneficial effects on the flow 
in the wake were expected to reduce the flow losses 
to a low level. 


In the design of the Caribou propeller-turbine in- 
stallation, flow calculations carried out by the electri- 
cal analogue method have contributed to the rapid 
and economic development of a high efficiency airflow 
system. It is hoped that this brief presentation will 
contribute to a more widespread application of the 
analogue method to the design of high efficiency flow 
systems. 
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SECRETARY’S LETTER 


Mr annual tour of the western Branches has always 

been a most enjoyable interlude and since I have 
been making these tours with the President they have 
become even more fun. This year the tour took place 
from the 16th to the 20th January; the weather was 
kind to us — it fact it was bright and clear all the time 
— and all our journeyings worked out as planned. 


VANCOUVER 


The President and I joined forces in Vancouver 
on the Monday. He had come up from a business trip 
to the States and I had come by DC-8 (what an aero- 
plane!) across Canada. On my arrival and before | 
met the President, Mr. Hartley, the Branch Chairman, 
took me around the Canadian Pacific Air Lines shops. 
Spending, as I do, so much time in Headquarters, these 
rare glimpses of reality are most refreshing. 


In the afternoon we had a meeting with the 
Branch Executive Committee and in the evening we 
attended the January meeting of the Branch, which 
had been rearranged from the 18th to enable us to be 
present. In Vancouver, as in most of our Branches, 
I recognize so many faces nowadays — and forget so 
many names, alas — that much as I enjoy these meet- 
ings I always leave them feeling a little guilty at not 
having greeted personally some of my old friends. | 
hope that they will understand how bewildering this 
sudden reunion with fifty or sixty people can be. 


CALGARY 


Next morning bright and early we set out for 
Calgary, where we were met by Mr. Fenby and 
and F/O Flaherty, the Chairman and Secretary of the 
Branch, respectively. The President went on into town 
but the rest of us waited at the airport to meet Dr. 
Whyte of the NRC, who was coming from Ottawa 
to speak to the Branch. Dr. Whyte and I were taken 
to lunch at Canadian Pacific Air Lines (Repairs) at 
Lincoln Park. I had not been there before and was 
interested to see this modern version of the “Mod 
Centre”, which I knew so well in the States during 
the war. 
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Later in the afternoon we had an informal meeting 
with the Branch Executive Committee in the Presi- 
dent’s room in the hotel, a room I might add with a 
magnificent view of the foothills country, rimmed by 
the white peaks of the Rockies across the entire 
horizon. The Branch held a dinner meeting in the 
evening. I think that the highlights were an unex- 
pected burst of French from the President in opening 
his excellent (I nearly said “otherwise excellent”) ad- 
dress — this was brought on by a few pointed remarks 
about la Belle Province by Mr. Fenby in introducing 
him — Dr. Whyte’s paper, and the customary western 
hospitality and renewal of acquaintances. It was a 
most enjoyable visit. 


COLD LAKE 


On the Wednesday we were due at Cold Lake and 
the President, Dr. Whyte and I flew up to Edmonton 
in the morning. I had spoken to F/L Smith, the Chair- 
man of the Cold Lake Branch, on the telephone from 
Calgary and as a result it was decided that the Presi- 
dent and Dr. Whyte should remain in Edmonton; the 
membership of the Cold Lake Branch has fallen con- 
siderably during the last year, due to transfers and for 
other reasons, and it is becoming an increasingly diffi- 
cult job to get there, we felt that it would be rather 
an imposition on the President and the Speaker to ask 
them to tackle the journey for perhaps a rather small 
meeting. So I went alone; and, had a wonderful time. 


A seven hour bus ride took me there. This in itself 
was a worthwhile experience, running through that 
remote country with brief stops at half a dozen iso- 
lated communities, with their Ukranian churches and 
grain elevators along the railway. I was met by F/O 
Bailey, the Branch Secretary, who kindly entertained 
me to dinner at his home before we went on to the 
meeting. It was a small meeting, but very pleasant and 
chatty; S/L Fryers stepped into the breach as speaker, 
and discussed his recent paper on Pressure Altitude 
Altimetry, which was published in the December issue 
of the Journal. We also discussed some of the current 
difficulties of the Branch. . 
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Cold Lake: Mr. C. B. Jeffery, F/L N. H. Smith and 
the Taylorcraft — running 


I had intended to return to Edmonton by train 
but Mr. Jeffery offered to fly me back in his rather 
elderly Taylorcraft and I readily accepted. Accord- 
ingly Mr. Jeffery, F/L Smith and I duly assembled 
next morning on the frozen surface of Cold Lake 
itself (this is a little way from the RCAF Station) 
where some half dozen small aircraft, some on skis 
and some on wheels, were soaking in the near-zero 
temperature. Mr. Jeffery’s Taylorcraft is an admirable 
aeroplane but easy starting when cold is not one of its 
virtues. When we had got it started, it suffered rather 
severely from mag drop and we had to take all the 
plugs out. All this was a lot of fun, if a little chilling; 
another return to reality. I hadn’t swung a prop for 
about 25 years, and then on wet English grass rather 
than Northern Alberta ice. However we duly arrived 
at Edmonton, where I found the President wondering 
whether he should begin looking for a new Secretary 
or whether it was worth alerting Search and Rescue 
to look for the old one. 


EDMONTON 

Late that afternoon Mr. van Horne, Chairman of 
the Edmonton Branch, and Mr. Quick, the Chairman 
of his Programmes Committee, met us in the Presi- 
dent’s hotel room and we discussed Branch and Insti- 
tute affairs. After dinner together, with Dr. Whyte, 
we proceeded to the Branch meeting, where again we 
met many old friends; both the President and I had, 
of course, attended the Mid-season Meeting in Ed- 
monton last year and had very pleasant memories of 
Edmonton hospitality. 

This was another good meeting and fortunately 
the President was able to stay to the end, though he 
had to leave that night on his return journey to the 
east. 
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WINNIPEG 


Because the President could not be with me on the 
Friday and because Friday is a bad night for meetings 
and because we would all be back in Winnipeg again 
for the Mid-season Meeting, no Branch meeting had 
been arranged in Winnipeg when I arrived on the last 
visit of my tour. However I reached Winnipeg by 
mid-day and was taken by Mr. Baker, the Branch 
Secretary, to see the site of the forthcoming Mid- 
season Meeting and to talk over our plans with the 
hotel. This done, we went to the Flying Club where 
I conferred with the Committee which had been 
working on the Mid-season Meeting plans. At the time 
of writing, I think that all that can be done to make 
this Meeting a success has been done and the Branch 
is to be congratulated on the preparation it has made. 


I spent the evening with the Branch Executive 
Committee at the home of W/C Evans, their Chair- 
man, and we talked about all manner of administrative 
and policy matters, with a few digressions into such 
topics as man-powered flight and snowboats (which, 
for the uninitiated as I was, are in effect aero-engine- 
powered sleds). 


And so home by the late night plane. It was a good 
tour. Somewhere along the way someone mentioned 
that my correspondence lately had been so full of the 
Institute’s woes that they had felt it was high time 
for me to come and see them and be cheered up. 
Cheered up I certainly was; it is most encouraging to 
witness the work that the Institute is doing across the 
country and the momentum that the western Branches 
have built up. I am sure that the President was simi- 
larly impressed. 


DUES DUE 


March is the end of our fiscal year and April the 
beginning of the next. During March we send out next 
year’s dues bills and I would suggest that they should 
be paid promptly. Certainly the Bylaws provide for 
a period of grace — technically a member does not 
become “in arrears” until October — but if payment 
is put off till then, we have all sorts of complications 
about loss of privileges and reminder letters, which 
we should avoid if we can. And we can, if members 
would pay when they get their bills and before they 


forget all about them. 
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BRANCHES 


Vancouver 


Reported by M. G. Brechin 


January Meeting 
The monthly meeting of the Branch 


was held on January 16th, 1961, at 


the RCAF Officers’ Mess Sea Island, 
and was attended by 46 members and 
4 guests. 

Mr. F. L. Hartley, Branch Chair- 
man, called the meeting to order and 
introduced Mr. David Boyd, Presi- 
dent of the CAI, and Mr. H. C. Lutt- 
man, Secretary, to the members and 
guests present. 


Mr. Boyd addressed the meeting 
with a few informal words, saying 
how pleased he was to be able to 
meet the members of the Vancouver 
Branch. He went on to say that the 
CAI was going through a few diffi- 
cult years due to the “cut-back” in 
aircraft manufacturing in Canada, and 
“where we were approaching future 
plans on horseback we are now walk- 
ing”. However, membership in the 
Institute is growing steadily in spite 
of the difficulties. He urged more 
membership from the hangar people 
such as Mechanics, Inspectors and Air 
Engineers — people who know their 
way around aircraft. These people are 
as interested as the “brass hats and 
slide rule boys”. Our future lies in 
development — not in size but in how 
effective we are. Aeronautics is the 
people, knowledge and education in 
aviation. 





Vancouver Branch Executive Committee and guests 
(lto r) Mr. M. G. Brechin, S/L A. E. Falls, Mr. David Boyd, Mr. F. L. Hartley, 
Mr. H. C. Luttman, Mr. R. J. McWilliams, Mr. A. G. Reddy 
and Mr. G. W. T. Roper. 


Mr. Luttman, in his address to the 
membership, said how pleased he was 
to be back visiting the Vancouver 
Branch. He laid emphasis on partici- 
pation for membership from the 
hangar and the shortage of papers for 
the Journal, especially those papers 
which lend themselves to the interest 
of people on the hangar floor. 

Presentation of membership certi- 
ficate was made to Mr. C. F. Tinley 
by the President. 

The Programme Chairman, Mr. 
A. G. Reddy, presented a short re- 
port advising that the hopes of the 





Presentation of membership certificate to Mr. C. F. Tinley by the President 
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Vancouver Branch in obtaining Sir 
Frank Whittle to address the joint 
CAI/SAE meeting in February were 
still high. 

The programme for the evening 
consisted of a panel discussion on the 
subject, “Constant Frequency vs Fre- 
quency Wild Generating Systems for 
Aircraft”. The panel members were 
Mr. F. P. Schregardus from Boeing 
Aircraft and Mr. L. C. Bryan, being 
the proponents of constant frequency, 
and Mr. J. A. Love and Mr. R. L. 
Jubinville for frequency wild systems. 
The moderator, Mr. G. A. Worden, 
was hard pressed at times to keep the 
discussion at a constant frequency, as 
the lively discussion which followed 
could be described as frequency wild. 

The panel members were thanked 
by Mr. W. E. Regan for their par- 
ticipation and for the lively interest 
aroused. 


Calgary 
Reported by F/O L. A. Flaherty 
January Meeting 

The January meeting of the Cal- 
gary Branch was a dinner meeting 
held at the Al San Club on 17th Janu- 
ary, 1961, at 7.30 pm. The meeting 
was presided over by the Chairman, 
Mr. G. H. Fenby. There were 16 
members and 10 guests present. 
Among the guests, we were honoured 
in having our President, Mr. David 
Boyd, and the National Secretary, 
Mr. H. C. Luttman, present. 
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Before the introduction of the 
guest speaker, Mr. Luttman spoke 
briefly to the members on the general 
policies affecting the CAI. He gave 
a short résumé of the SAE Interna- 
tional Congress and Exposition which 
he attended in Detroit. He also ex- 
tended an invitation for the submis- 
sion of individual papers for the 
Journal, and explained how the new 
layout of the Journal would permit 
individual binding of papers. 

Mr. Boyd then addressed the mem- 
bers, beginning his address in French, 
but switching quickly to English, 
when he saw the consternation on 
the faces of some of the members. 
He suggested later, that we should 
all endeavour to learn French. Mr. 
Boyd expressed the hope that the 
CAI collectively and individually, by 
Branches and members, would en- 
deavour to take the lead in putting 
the spirit of romance back into avia- 
tion. He stressed the need for en- 
couraging students to be interested 
in the aircraft business. He urged the 
members to make an effort to in- 
crease the membership, and he touch- 
ed briefly on the financial situation of 
the CAI. Upon the conclusion of his 
address, he was thanked by Mr. 
Fenby. 


Mr. W. E. Jamison introduced the 
guest speaker, Dr. R. B. Whyte, 
whose subject was “Qualification 
Testing of Fuels and Lubricants at 
the NRC”. 


Dr. Whyte explained the methods 
and purposes for qualification testing, 
and supported his talk with slides 
showing the actual equipment in use 
at the National Research Laboratories 
in Ottawa. It was interesting to note 
that these facilities are now available 
to any oil company wishing to be- 
come a supplier to the armed forces. 
Through this facility, the company 
can obtain qualification testing and, if 
acceptable, be certified by that branch 
of the armed forces concerned. Prior 
to this facility being made available, 


CAI 


25th-26th May, 1961—Annual Gen- 
eral Meeting, Toronto, Ont. 


BRANCHES 


Vancouver 


22nd April — Canyon Garpens, An- 
nual General Meeting and Dance. 
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oil companies in Canada were at a 
considerable disadvantage. 

Dr. Whyte’s talk was followed by 
a question and answer period, at the 
conclusion of which S/L R. S. Crosby 
thanked Dr. Whyte on behalf of the 
Calgary Branch. 

The Chairman, Mr- Fenby, then 
extended the thanks of the Calgary 
Branch to the President and the Na- 
tional Secretary for visiting us and 
expressed the hope of continued 
liaison and good will. 


Halifax-Dartmouth 
Reported by Lt A. M. Percy 


January Meeting 

The regular meeting was held in 
the cinema of the CPOs’ Mess, HMCS 
Shearwater, on Wednesday, 18th 
January, 1961. A total of 45 members 
and guests were present. The meeting 
was presided over by LCDR G. M. 
Cummings, the Branch Chairman. 

The Chairman stated that CPO 
G. E. Blackwell had offered his ser- 
vices as Secretary since the former 
Secretary, CPO A. C. Green, was 
now serving at sea, and that Lt A. M. 
Percy had been elected as Vice 
Chairman of Publicity. 

The Chairman then invited Mr. 
W. G. Stewart to introduce the guest 
speaker, Mr. L. A. Dickinson. Mr. 
Dickinson, who is Head of the 
Rocket Engine Development Section, 
Propulsion Wing, CARDE, gave a 
very interesting and informative lec- 
ture entitled “Rocket Engine De- 
velopment in Canada”. 

The speaker presented his lecture 
in two sections. The first, which was 
illustrated by a number of slides, was 
devoted to describing the factors in- 
fluencing the design of a rocket motor 
and the choice of a suitable propel- 
lant to satisfy the many and varied 
requirements. The second section was 
a film showing stages in the produc- 
tion and testing of a rocket propul- 
sion vehicle of the Black Brant type 


COMING EVENTS 


Tour Speakers 


Advanced Systems Planning Looks at 
VTOL in Peace and War, 
W. S. Geppes. 


20th March — Vancouver 
21st March — Edmonton 
22nd March — Calgary 
23rd March — Winnipeg 


culminating in the successtul launch- 
ing of such a vehicle. 

The lecture was followed by a 
lively discussion period during which 
Mr. Dickinson mentioned that Cana- 
dian research was now taking place 
to investigate the cheniistry of the 
upper atmosphere, with a view to 
improving long range weather fore- 
casting and long range radio com- 
munications. 

Mr. R. Wallworth thanked the 
speaker for the interesting presenta- 
tion of his lecture and for the ready 
and informative answers to the varied 
questions raised during the discussion 
period. 


Cold Lake 
Reported by F/O D. N. Bailey 


January Meeting 

The January meeting of the Branch 
was held in the GIS Theatre on the 
18th January, 1961. 7 members and 
2 guests attended. 

After welcoming those present, the 
Branch Chairman, F/L N. H. Smith, 
introduced the Secretary of the CAI, 
Mr. H. C. Luttman, who gave a short 
talk on the purpose and aims of the 
CAI. This was followed by a short 
discussion of Branch problems, with 
special reference to the future of the 
Cold Lake Branch. 

The speaker for the evening, S/L 
W. R. Fryers, was introduced by Mr. 
C. B. Jeffery. The speaker’s subject 
“Pressure Method of Measuring Al- 
titude”, covered the research and 
planning he had done on the measur- 
ing of aircraft altitude with reference 
to pressure levels rather than to feet 
or kilometers. He had obviously gone 
into the problem thoroughly and had 
solutions to the problems of aircraft 
flight level separation, a standard 
flight level indication and had a ten- 
tative design for method of presenta- 
tion to the pilot. 

After an informal discussion period, 
the speaker was thanked on behalf of 
those present by F/O D. N. Bailey. 


A Review of Technical Requirements 
for the Operation of Aircraft at Sea, 
CDR J. F. Frank. 

22nd March — Montreal 

23rd March — Quebec 


Solutions to Aircraft Store Separation 
Problems, T. B. FessenpEN. 


19th April — Halifax-Dartmouth 
20th April — Quebec 
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BOOKS 


Transport Phenomena. By R. B. Biro, 
W. E. Srewart anv E. N. LicHtroot. 
John Wiley and Sons, New York, 1960. 


‘780 pages. Illus. $13.75. 


This book, written by three mem- 
bers of the Chemical Engineering 
Department at the University of 
Wisconsin, is an important addition 
to the literature of fundamental en- 
gineering science. It deals with those 
physical phenomena in which the 
transport in fluid or solid media of 
momentum, energy or mass are the 
central features. Thus the topics dealt 
with include, among others, viscous 
flows (momentum transport), heat 
transfer (energy transport) and dif- 
fusion (mass transport). The close 
correspondence and mathematical 
analogies among these transport pro- 
cesses are emphasized. The material 
is cleverly organized on a “two-co- 
ordinate” basis; that is, each of the 
22 chapters is defined by (a) the 
type of transport and (b) the entity 
being transported. The entities are 
those already mentioned (momen- 
tum, energy and mass) and the types 
of transport include transport by 
molecular motion, by laminar and 
turbulent flow, transport across fluid 
and solid interfaces, and transport by 
radiation. These and other associated 
phenomena are basic to many engi- 
neering processes and design prob- 
lems. They have widespread and im- 
portant applications in aeronautics 
and astronautics, of which one note- 
worthy example is the re-entry heat- 
ing problem. In the case of ablation, 
all three processes take place simul- 
taneously and are strongly coupled. 


Since the book is intended for use 
as a text book, many excellent prob- 
lems are included at the end of each 
chapter to be worked by the student. 
These should prove to be of great 
help to teachers. 


The level of scholarship is gratify- 
ingly high. The authors exhibit a 
comprehensive and fundamental grasp 
of an extremely wide variety of sub- 
ject matter, and have effectively put 
this across in the text. However, it 
must be said that students will not 
find this book ‘easy going’. The de- 
tails of the mathematical arguments 
are rather condensed, and the level 
of mathematics employed is quite ad- 
vanced (but not unnecessarily so). 


March, 1961 


This book may be said in a sense 
to be three in one — Viscous Flow, 
Heat Transfer and Mass Transfer. 
The case for a unified treatment in 
a single volume has been well put by 
the authors, and this reviewer has no 
doubt of the value of this approach 
in a work of reference. He has some 
reservations as to its desirability from 
the pedagogical standpoint, however. 
The volume of material included in 
this way may be too much for one 
course. Experience with classroom 
use will show how correctly the 
authors have anticipated future de- 
velopments in engineering curricula. 
Indeed this important book may exert 
a substantial influence on them. 

B. Erkin 


An Introduction to the Theory of Air- 
eraft Structures. By Dr. D. WituiaMs. 
MacMillan Co. Canada, 1960. 448 pages. 
Illus. $10.25. 


The first sentence of the preface 
sums up this book admirably: “The 
aim of this book is to provide a back- 
ground for the solution of aircraft 
structural problems as distinct from 
providing a collection of solutions to 
individual problems”. 


Dr. Williams is Deputy Chief 
Scientific Officer in the Structures 
Department of the Royal Aircraft 
Establishment, Farnborough. He has 
had many papers published in various 
technical journals, plus many RAE 
reports, to establish him as an author- 
ity on the subject of aircraft struc- 
tures. 


Chapter 1 presents a sufficient 
amount of the theory of elasticity 
that is required later in the book. It 
is well presented and easy to follow. 
The two chapters on strain energy 
methods and strain energy examples 
are a welcome addition in books of 
this type. Too often strain energy 
theories are derived and it is up to 
the student to apply them. Chapters 4 
and 5 deal with “individual displace- 
ments” and application and removal 
of constraints. This latter chapter 
only touches on the basic elements of 
the subject, with the author himself 
suggesting reference to Southwell’s 
work on relaxation methods for a 
more comprehensive background. 
Chapter 6 deals with torsion of shells 





and box beams and considers axial 
warping of box beams. Chapter 7 
deals with bending shear stresses in 
thin walled beams, both open and 
closed sections. Both these chapters 
contain much that can be found in 
appropriate textbooks. Chapter 8 deals 
with shear lag and is presented in an 
excellent manner. Rectangular box 
beams and beams representative of 
fuselages are considered. An excellent 
collection of references on shear lag 
are added at the end of the chapter. 
Chapter 9 deals with pressure cabin 
problems and includes such items as 
neutral holes, reinforcing around 
holes and effect of cabin frames and 
bulkheads. An interesting example of 
a typical fuselage shows that conven- 
tional frames at 20 inch pitch can at 
best be only 38% effective in helping 
to reduce the maximum hoop stress 
in the skin. It was disappointing to 
see only about four pages on the sub- 
ject of fatigue and fail-safe aspects 
of pressure cabins. Dr. Williams could 
have reduced the contents of the last 
four chapters in order to discuss this 
subject in more detail. The last four 
chapters deal with the dynamics of 
the aeroplane. They will be especially 
useful to the aircraft stressmen in 
providing a good background on the 
problems his neighbour engineer in 
the dynamics department is confront- 
ed with. More could have been in- 
cluded on power spectrum analysis 
of gusts since this “new” technique is 
now an accepted procedure of 
analysis. Throughout the book, Dr. 
Williams has emphasized the import- 
ance of stiffness as well as strength. 
In present day design, for instance, 
wing flexibility is one of the most im- 
portant problems, and therefore the 
addition of these four chapters is not 
out of place in a structures book of 
this kind. 


The mathematics involved are of 
such a level that the graduate engi- 
neer should not have too much 
trouble in reading this book. 


The book is one which can be re- 
commended, to the post-graduate en- 
gineer and the aircraft stressman — 
the former for studying, the latter as 
a reference book. 


E. AUBREY 
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Avionics Research: Satellites and Pro- 
blems of Long Range Detection and 
Tracking. Edited by E. V. D. Grazirr, 
E. Ricurer anp J. Voce. Pergamon Press, 
1960. 257 pages. Illus. $10.00. 


The 19 papers published in this 
book were presented at the AGARD 
Symposium held in Copenhagen in 
October, 1958, under NATO aus- 
pices. 

Immediately after the unannounced 
launching of Sputnik I, radio astro- 
nomy equipments were converted to 
radar sets and the Minitrack stations 
modified to operate on 20 mc/sec and 
40 mc/sec. Eight of the papers dis- 
cuss the results obtained. Radar de- 
tection was not reliable; its failures 
being attributed to polarisation rota- 
tions in the ionosphere and to satel- 
lite motions causing the effective 
echoing area to fluctuate over wide 
limits. Such motions could also be 
inferred from the amplitude varia- 
tions of the signals from the satellite 
radio transmitters. 


The facilities of the Journal are offered 
Gee of ete So See mentee Oe 


notice of positions 4 
published for two consecutive months and 
will thereafter be oe unless their 
reinstatement is ee © 


Box No., to which cxpuinios 
dressed (c/o The Secretary), 33 


- 
to each notice submitted Z an 
individual. 


The Institute reserves the right to decline 
onsidered unsuitable for this 


tion if circumstances so 


POSITION WANTED 


Box 111 Technician: Aircraft Main- 
tenance Technology graduate from 
the Southern Alberta Institute of 
Technology and Art desires position 
in industry. Will appreciate all offers. 
Excellent references available. 


SUMMER EMPLOYMENT 


Box 112 Student: Second year engi- 
neering student at Ecole Polytech- 
nique de Montréal seeks position for 
the period 15th May to 15th Septem- 
ber. Bilingual with technical drawing 
and surveying experience. 
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Detection of large ground targets 
by ionospheric radar is discussed and 
mention is made of the probable de- 
tection of a satellite using this tech- 
nique. Some observers, however, feel 
that the echo was returned from the 
ionization trail behind the satellite, a 
view which could be supported by 
the data on auroral echoes and the 
responses from meteor trails. 

Other papers deal with the electro- 
hydrodynamic effects of an ionized 
medium on the motion of the satellite. 
The connection between this and the 
detection problem seems somewhat 
tenuous; and a similar criticism could 
be levelled at the description of the 
re-entry of an IRBM where initial 
detection was by a chance visual ob- 
servation. 

A paper on echoing area charac- 
teristics deals with simple shapes, 
meteor trails and the size of the 
moon. Unfortunately the continua- 
tion of the discussion to artificial 
satellites was not attempted. 


APPOINTMENT NOTICES 


Box 113 Student: First year student 
in Aircraft Maintenance Technology 
at the Southern Alberta Institute of 
Technology and Art seeks position, 
in western Canada, for the period 
May 15th to September Ist. 


Box 114 Student: First year student 
in Aircraft Maintenance Technology 
at the Southern Alberta Institute of 
Technology and Art seeks position, 
preferably in Ontario, for the period 
May 15th to September 15th. 


Box 115 Student: Student having 
completed his first year in the Air- 
craft Maintenance Technology course 
at the Southern Alberta Institute of 
Technology and Art seeks position in 
the Edmonton or Calgary areas. 


Box 116 Student: Second year post- 
graduate student in Mechanical Engi- 
neering at Nova Scotia Technical Col- 
lege seeks position from the ist May 
in the propulsion field, with an out- 
look for permanent employment after 
graduation. No geographical prefer- 
ence. Experienced in supervision. 


The theory and philosophy of 
radar detection receive cursory treat- 
ment in two other papers. One in- 
teresting paper describes the use of 
noise radiated by the sun as a remote 
signal source in the measurement of 
antenna patterns, particularly when 


ground reflection effects are encoun- 
tered. 

The “principal purpose of the 
Avionics Panel of AGARD is to 
highlight profitable areas of avionics 
research and development”. This aim 
is not realized in these papers, other 
than to point out the incomplete 
knowledge, at the time, of what 
might be loosely called “propagation 
phenomena”. The book is not likely 
to provide much background data for 
aeronautical engineers unless many of 
the references are consulted; for the 
electronics specialist recent literature 
supplants much of the information 
presented. 


J. CuisHoLm 


Box 117 Student: First year student 
in Aircraft Maintenance Technology 
at the Southern Alberta Institute of 
Technology and Art seeks position, 
preferably in western Canada, for the 
period March 20th to September Ist. 
Hard worker and willing to learn. 


Box 118 Student: Second year stu- 
dent in Aeronautical Engineering at 
the Southern Alberta Institute of 
Technology and Art seeks position, 
preferably in western Canada, for the 
period May 15th to September 15th. 
Experience includes two years Mech- 
anical, structural and topographical 
drafting. 


Box 119 Student: First year student 
at the Southern Institute of Tech- 
nology and Art seeks position any- 
where in western Canada. Experience 
includes a Private Pilot’s Licence and 
two months service in the RCAF 
Reserve Tradesmen Training Plan to- 
wards an Airframe Technician Group 
I rating. 

Box 102 Student: Fourth year Mech- 
anical Sciences student at McGill 
University seeks position in plant 
engineering. Experience includes four 
years workshop and drafting prior to 
entering university. 
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BL-1705 Mounting System — shown 
protecting pressure transducer. 


in control for aircraft instruments 


e’s a new Lord mounting system for small equipment or 
instruments on jet aircraft that offers big advantages! 






i 2 Ot 2 Sek ea a mee = = oo tt =e eae oo oe ae ' 
It’s versatile. Basic design is adaptable to a variety of : : 5 
units: transducers, indicators, electronic tubes, gages, warning 2 typical speci fications s 
and timing devices, baroswitches. Mounting i este : f f 
arrangements permit use in instrument panel, nacelle : Application: pressure transducer. Weight: equipment = § 
or ates aietrese tocations. i = 1.2 Ibs., base — .38 Ibs. Shock/vibration protec: f 
® tion: all-attitude. Isolators: four Special BTR Multiplane 2 
It's high-performance. Lightweight system provides » Mountings. System natural frequency: 45 cps. Operat- * 
excellent all-attitude control of high-frequency excitations plus , ing temperature range: —65° to +300°F. Environ- i 
attenuation of shock and structure-borne noise. ; mental resistance: unaffected by fungi, dust, sand, 
BTR® Elastomeric Mountings possess extreme environmental 5 salt atmosphere, oil, Saome. Construction: centet-of- 
resistance, excellent damping, superior endurance, J gravity suspension, fail-safe. Mounting arrangemers: 
; > * clearance holes or clinch nuts for base, bulkhead or @ 
consistent performance over —65° to +300°F. range. é overhead attachment. ' 
It’s economical. Simplified design gives you Lord quality i aac meaeer i 
at a lower price than competitive bases. Long é al | | piil | Ls a 
service life means your maintenance costs will be lower, too. L > VERTICAL & | INPUT: 5 IN./SEC.| | i 
It's proved. The advanced design and performance ' = 200 , 1 
of this Lord system have been thoroughly proved in actual 3 a 3 
service on today’s jet airliners. j : 100 i 
w 
This mounting system is an example of Lord ingénuity. i < 0 t 
It indicates why you can continue to expect Lord to produce the a 2 3040 6 8 ie | 
best in vibration/shock/noise control for the aerospace ' i 
environment. To put this ingenuity to work on your project, j j 
or to get further data on the BL-1705 base, contact —— a 


the nearest office of Railway & Power Engineering Corporation, 
Limited. 
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